Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,539 T→C 11.3% I474T (ATC→GCC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,5390T→C11.3% 58.2 ‑1.2 ‑1.1296I474T (ATC→ACC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (10/16):  ref base (118/114):  total (145/151)
Fisher's exact test strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.29e-02

Only 100 of 316 total aligned reads displayed.

TCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTC  >  NC_014034/1397439‑1397639
                                                                                                    |                                                                                                       
tCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGAt                                                                                                         >  2:1136457/1‑101 (MQ=255)
 ccccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATc                                                                                                        <  2:3923732/101‑1 (MQ=255)
   ccccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCAt                                                                                                      <  1:1994561/101‑1 (MQ=255)
    cccGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                                                                                                     <  2:3734429/101‑1 (MQ=255)
       gcgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c                                                                                                 <  2:666211/101‑1 (MQ=255)
        cgGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGc.c..c                                                                                              <  1:3032414/101‑1 (MQ=255)
         gggggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..Cg                                                                                             >  2:4089654/1‑101 (MQ=255)
           gggCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGcc                                                                                           <  2:173657/101‑1 (MQ=255)
             gCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTc                                                                                         >  2:1900110/1‑101 (MQ=255)
                tGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgc                                                                                      <  2:4141995/101‑1 (MQ=255)
                tGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgc                                                                                      >  1:762079/1‑101 (MQ=255)
                 ggATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAgcg                                                                                     >  1:579157/1‑101 (MQ=255)
                      gAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCcgcg                                                                                >  1:2818027/1‑101 (MQ=255)
                      gAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCcgcg                                                                                >  1:2193962/1‑101 (MQ=255)
                        aaCACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCt                                                                              >  1:1868301/1‑101 (MQ=255)
                         acacCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTg                                                                             <  2:3436610/101‑1 (MQ=255)
                           cccGTCATCACCACGCAACCCCCCTGTACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGAc                                                                           <  1:3737700/100‑1 (MQ=255)
                           acCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGAc                                                                           <  1:889485/101‑1 (MQ=255)
                           acCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGAc                                                                           >  2:2120712/1‑101 (MQ=255)
                             cGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGt                                                                         >  1:165736/1‑101 (MQ=255)
                              gTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTc                                                                        >  1:504215/1‑101 (MQ=255)
                                 atcaCCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCt                                                                     >  1:1196816/1‑101 (MQ=255)
                                   caccacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcc                                                                   <  1:1001709/101‑1 (MQ=255)
                                    accacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTccc                                                                  >  1:692614/1‑101 (MQ=255)
                                      cacGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTccccc                                                                <  1:131714/101‑1 (MQ=255)
                                       ccGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcccccc                                                               <  1:1448125/100‑1 (MQ=255)
                                       acGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTcccccc                                                               >  2:4002436/1‑101 (MQ=255)
                                        cGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcg                                                              <  2:1955200/101‑1 (MQ=255)
                                         gCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcgc                                                             >  1:2925137/1‑101 (MQ=255)
                                          cAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcgcg                                                            <  1:2593696/101‑1 (MQ=255)
                                          cAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcgcg                                                            >  1:3726973/1‑101 (MQ=255)
                                             ccGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCGTGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgggg                                                         <  1:3932774/101‑1 (MQ=255)
                                             ccGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgggg                                                         <  1:1563843/101‑1 (MQ=255)
                                              cgctcggccccaccGGGCGCTCCCCCCGCGGGCGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTCACGTCGCTCCCCCCGCggggg                                                        <  1:2545808/92‑1 (MQ=37)
                                              cGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGTGCCTCGCTGACGTCGCTCCCCCCGCggggg                                                        >  2:1468803/1‑101 (MQ=255)
                                                 ccTGCCCCCCCGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGt                                                     <  1:67039/101‑1 (MQ=255)
                                                 ccTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGt                                                     >  2:3495317/1‑101 (MQ=255)
                                                  cTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTg                                                    >  2:3158082/1‑101 (MQ=255)
                                                  cTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTg                                                    >  2:1306909/1‑101 (MQ=255)
                                                   tGCCCCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTgg                                                   <  1:3099114/101‑1 (MQ=255)
                                                   tGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTgg                                                   >  1:318928/1‑101 (MQ=255)
                                                   tGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCCCGGGGGCGTgg                                                   >  1:1708627/1‑101 (MQ=255)
                                                     caccacCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGAt                                                 <  1:1468908/101‑1 (MQ=255)
                                                       ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg                                               <  2:2282511/101‑1 (MQ=255)
                                                        cgagcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                               <  1:1361849/96‑1 (MQ=0)
                                                        aggatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCCCGGGGGCGTGGAtgaca                                             >  2:3612589/6‑97 (MQ=0)
                                                          gatgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGGGGGTCGGAAc                                           >  2:2431248/4‑101 (MQ=0)
                                                           cGGTCGCTCCACCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                           >  1:3849619/1‑101 (MQ=255)
                                                            ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                          <  1:647964/101‑1 (MQ=255)
                                                             gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                          >  2:2944194/1‑101 (MQ=2)
                                                             gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                         >  1:4216445/1‑101 (MQ=255)
                                                              tcccTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                        <  2:2036093/98‑1 (MQ=255)
                                                              tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                        >  2:1209741/1‑101 (MQ=255)
                                                              tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                        >  2:4167502/1‑101 (MQ=255)
                                                               ctcTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcc                                       <  2:1838111/99‑1 (MQ=255)
                                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCCTGGATCCAAAcggggt                                      >  1:2113117/1‑96 (MQ=2)
                                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                       >  2:4238696/1‑100 (MQ=2)
                                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcg                                       <  2:2836700/101‑2 (MQ=2)
                                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcc                                       <  1:1926276/101‑1 (MQ=255)
                                                               cGCCCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaa                                        <  1:3568048/101‑3 (MQ=0)
                                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACaaaa                                       >  2:1595413/1‑98 (MQ=2)
                                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcg                                      >  1:2392472/1‑101 (MQ=255)
                                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcg                                      >  2:2977803/1‑101 (MQ=255)
                                                                ccTCCCCCCGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcg                                     <  2:578656/100‑7 (MQ=2)
                                                                 ccccccccGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaa                                      <  2:4411418/99‑5 (MQ=0)
                                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaa                                      >  2:3279913/1‑97 (MQ=2)
                                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaa                                      >  2:3063328/1‑97 (MQ=2)
                                                                  tCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatg                                    >  1:4346928/1‑97 (MQ=2)
                                                                  tCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAATcgatg                                    >  2:2410476/1‑97 (MQ=2)
                                                                  tCCCCCCGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaag                                     >  2:879322/1‑96 (MQ=0)
                                                                  tCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCcgcg                                    >  2:2156565/1‑101 (MQ=255)
                                                                  cccccccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatg                                    <  1:1595519/100‑5 (MQ=2)
                                                                   ccccccGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaagc                                    <  1:3140497/101‑7 (MQ=0)
                                                                   ccccccGCGGGGGCGTGGATCGAAACACCAATGCCGAAGC.CAAGGCCGCAGCAGCG.GCTGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACaaaaagc                                    >  1:3327065/1‑95 (MQ=0)
                                                                    cccccGCGGGGGCGTGGATCGAAACTATGCTGGC.TCGA.C..CGAGTTCACGGCGGGCTATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgggtcgaaga                                 <  2:3531133/101‑11 (MQ=0)
                                                                     ccccGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGAt                                 <  1:4136213/101‑1 (MQ=255)
                                                                      cccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatca                                >  1:2726014/1‑93 (MQ=0)
                                                                      cccGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG..CGGC.GAACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATcgatgatca                                <  2:3311248/101‑9 (MQ=0)
                                                                          cgGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCg                            >  2:2883607/1‑101 (MQ=255)
                                                                          cgGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCg                            >  1:1776753/1‑101 (MQ=255)
                                                                           gggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGt                           >  1:2130288/1‑101 (MQ=255)
                                                                           gggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGt                           >  1:99937/1‑101 (MQ=255)
                                                                            tctgCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTg                          <  1:235923/98‑1 (MQ=255)
                                                                            ggggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTg                          <  1:1134708/101‑1 (MQ=255)
                                                                             gggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGa                         <  2:793693/101‑1 (MQ=255)
                                                                             gggCGTGGATCGAAACAACCCGATCATGGC.C..CCCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGGGGATCGAAAAATCCGCGGATCGGCGTgc                         >  1:1087594/1‑100 (MQ=255)
                                                                              ggCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAc                        >  1:4282269/1‑101 (MQ=255)
                                                                                cGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGGGAccc                      >  2:1122285/1‑101 (MQ=255)
                                                                                 gTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccg                     <  2:2700299/101‑1 (MQ=255)
                                                                                     aTCGAAACACCCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc                 <  1:2268745/101‑1 (MQ=255)
                                                                                     aTCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgcccc                 <  1:3510730/101‑1 (MQ=255)
                                                                                      tCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgccccc                >  2:4369634/1‑100 (MQ=255)
                                                                                       cGAAACACCCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgccccgc               <  2:2193823/101‑1 (MQ=255)
                                                                                       cGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGAcccgccccgc               <  2:547999/101‑1 (MQ=255)
                                                                                        gAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACccgccccgcc              >  1:2682028/1‑101 (MQ=255)
                                                                                         aaaCAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCcgccccgccc             <  2:2391187/101‑1 (MQ=255)
                                                                                            caacCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATc          <  2:1676376/101‑1 (MQ=255)
                                                                                             aacCCGATCATGCC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCg         <  2:694784/101‑1 (MQ=255)
                                                                                               cccGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGt       <  1:4015748/101‑1 (MQ=255)
                                                                                                 cGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCCTCGGTCg     >  1:3533363/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
TCCCCCCGCGGGGGCGTGGATCGAAACACCGTCATCACCACGCAACCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGC.C..CGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCGGATCGGCGTGACCCGCCCCGCCCATCGGTCGCTC  >  NC_014034/1397439‑1397639

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG