Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,918,304 A→G 16.6% intergenic (‑39/+12) atpG ← / ← atpA F1 sector of membrane‑bound ATP synthase, gamma subunit/F1 sector of membrane‑bound ATP synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,918,3040AG16.6% 372.1 / 65.8 168intergenic (‑39/+12)atpG/atpAF1 sector of membrane‑bound ATP synthase, gamma subunit/F1 sector of membrane‑bound ATP synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (66/74);  new base G (14/14);  total (80/88)
Fisher's exact test for biased strand distribution p-value = 8.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.55e-01

Only 100 of 170 total aligned reads displayed.

GATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATC  >  NC_000913/3918206‑3918403
                                                                                                  |                                                                                                   
gATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAg                                                                                                    <  2:1522863/100‑1 (MQ=255)
    ttttGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGgccgc                                                                                                <  1:3635965/100‑1 (MQ=255)
     tttGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGgccgc                                                                                                <  1:3852793/99‑1 (MQ=255)
     ctcttccgatCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGgccgca                                                                                               <  1:3732342/91‑2 (MQ=255)
       tGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAg                                                                                             <  2:3780289/100‑1 (MQ=255)
       tGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAg                                                                                             <  2:270538/100‑1 (MQ=255)
       tGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAg                                                                                             >  2:2091615/1‑100 (MQ=255)
        gCGTGTTCTGGACGCTTGCGATCTTACTTCGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTGGCGGCTTGCCTTAGGGCAGGCCGCAAgg                                                                                            <  1:2260970/100‑2 (MQ=37)
        gCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAgg                                                                                            <  2:502412/100‑2 (MQ=38)
        gCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAgg                                                                                            <  1:3318653/100‑2 (MQ=38)
        gCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAgg                                                                                            <  2:3224610/100‑2 (MQ=38)
        gCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAgg                                                                                            <  2:1767230/100‑2 (MQ=38)
        gCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAgg                                                                                            <  2:2421725/100‑2 (MQ=38)
        gCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGa                                                                                            <  1:3019690/100‑1 (MQ=255)
         cGTGTTCTGGACGCTTGCGATCTTACGACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGAc                                                                                           <  2:2007658/100‑1 (MQ=255)
          gtgtTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACg                                                                                          >  1:1067791/1‑100 (MQ=255)
          gtgtTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACg                                                                                          <  2:2050244/100‑1 (MQ=255)
            gtTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGtt                                                                                        >  2:1087390/1‑100 (MQ=255)
               cTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTAcc                                                                                     >  1:257759/1‑100 (MQ=255)
                 ggACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAg                                                                                   >  1:791317/1‑100 (MQ=255)
                  gACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAgg                                                                                  >  1:1801550/1‑100 (MQ=255)
                      cTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGAtt                                                                               >  2:3138148/1‑99 (MQ=255)
                       ttGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCTCCAGACGTTACCAGGATTgg                                                                             >  1:1317229/1‑100 (MQ=255)
                        tGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTgg                                                                             >  1:60832/1‑99 (MQ=255)
                        tGCGATCTTACTACGTATCTATTTTGCGCCGGCCATGAGCTTCTCCTCAATGCATTGCGGCCTGCCCTAAGGCAAGCCGCCAGCCGTTACCAGGCTTggg                                                                            >  1:536611/1‑100 (MQ=255)
                         gCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGt                                                                           <  2:1180283/100‑1 (MQ=255)
                          cGATCTTACTACGTATCTCTTTTTCTTTTTTCATTTGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGttt                                                                          >  2:2492755/1‑100 (MQ=255)
                          cGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGtt                                                                          <  1:866186/100‑1 (MQ=255)
                          cGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGtt                                                                          >  1:1953670/1‑100 (MQ=255)
                            aTCTTACTACGTATTTCTTTTGCGCCGGCCATGGGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGc                                                                        <  2:1406466/100‑1 (MQ=255)
                            aTCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGc                                                                        >  1:79099/1‑100 (MQ=255)
                            aTCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGc                                                                        >  1:550286/1‑100 (MQ=255)
                             tCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCt                                                                       >  1:1470661/1‑100 (MQ=255)
                              cTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCtt                                                                      <  1:313108/100‑1 (MQ=255)
                              cTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCt                                                                       >  2:2495797/1‑99 (MQ=255)
                                tactacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTg                                                                    >  2:563197/1‑100 (MQ=255)
                                 actacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGAAAGCCGCCAGACGTTACCAGGATTGGGTTGCttt                                                                     >  2:3323846/1‑98 (MQ=255)
                                  ctacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGaa                                                                  <  1:1087397/100‑1 (MQ=255)
                                  ctacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGaa                                                                  >  1:1910616/1‑100 (MQ=255)
                                  ctacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGaa                                                                  >  1:2665963/1‑100 (MQ=255)
                                   tacGTTTCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTCCCCTAAGGCAAGCCGCCAGACGTTACCAGGATGGGGTTGCTTTGaa                                                                  >  2:3359233/1‑99 (MQ=255)
                                   tacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAg                                                                 <  1:226447/100‑1 (MQ=255)
                                    acGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAgg                                                                <  2:1794985/100‑1 (MQ=255)
                                     cGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGa                                                               >  1:54657/1‑100 (MQ=255)
                                     cGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGa                                                               >  2:1243641/1‑100 (MQ=255)
                                       tATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAAt                                                             >  1:2827575/1‑100 (MQ=255)
                                        aTCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATc                                                            <  1:2304052/100‑1 (MQ=255)
                                          ntctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGa                                                          >  2:2928546/2‑100 (MQ=255)
                                           tctTTTGCGCCGGCCCTGAGCTTCTCCTCACTGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGCTTGCTTTGAAGGAAGCGAg                                                         >  2:3730200/1‑100 (MQ=255)
                                           tctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAg                                                         <  1:3544042/100‑1 (MQ=255)
                                              tttGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGCAGGAATCGAggtt                                                      <  1:2673371/100‑3 (MQ=255)
                                                  cgcCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGTCGCCAGACGTTACTAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCt                                                  >  2:587865/1‑100 (MQ=255)
                                                      ggCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCa                                              <  1:687751/100‑1 (MQ=255)
                                                      ggCCATGAGCCTCTCCTCAATGCCTTGCCGCCTGCCCTAAGTCAAGCCGCCAGACGCTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCa                                              <  2:899077/100‑1 (MQ=255)
                                                        ccATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGc                                            <  1:2120339/100‑1 (MQ=255)
                                                             aGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGcc                                       <  2:504255/100‑1 (MQ=255)
                                                              gCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCt                                      <  1:202619/100‑1 (MQ=255)
                                                                ttCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTc                                    <  2:1446019/100‑1 (MQ=255)
                                                                   tcctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAt                                 <  2:1404007/100‑1 (MQ=255)
                                                                   tcctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAt                                 <  1:3012980/100‑1 (MQ=255)
                                                                    cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAtt                                >  1:1639322/1‑100 (MQ=255)
                                                                    cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAtt                                >  1:120788/1‑100 (MQ=255)
                                                                     ctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAttt                               <  2:1669569/100‑1 (MQ=255)
                                                                     ctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAttt                               >  2:307309/1‑100 (MQ=255)
                                                                      tcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtc                              >  1:3274331/1‑100 (MQ=255)
                                                                       cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcg                             <  2:3601559/100‑1 (MQ=255)
                                                                       cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcg                             >  2:2780225/1‑100 (MQ=255)
                                                                       cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcg                             <  2:2232404/100‑1 (MQ=255)
                                                                          tGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcgtcg                          <  1:475615/100‑1 (MQ=255)
                                                                              ttGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGt                      >  2:231493/1‑100 (MQ=255)
                                                                               tGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTa                     >  2:1417601/1‑100 (MQ=255)
                                                                               ggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGCTTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTa                     >  2:840111/2‑100 (MQ=255)
                                                                               ggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTa                     >  1:740661/2‑100 (MQ=255)
                                                                               ggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTa                     >  1:197767/2‑100 (MQ=255)
                                                                               ggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTa                     >  2:127801/2‑100 (MQ=255)
                                                                               ggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTa                     >  2:217902/2‑100 (MQ=255)
                                                                                gCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGACGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGt                      >  2:2391499/1‑98 (MQ=255)
                                                                                gCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAg                    >  1:1897108/1‑100 (MQ=255)
                                                                                gCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCCTGCCTTCGATTTCGTCGTTGTAg                    >  1:2552048/1‑100 (MQ=255)
                                                                                 cGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTGGc                   >  2:1689873/1‑100 (MQ=255)
                                                                                  ggCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGcc                  >  2:2748873/1‑100 (MQ=255)
                                                                                  ggCCTGCCCTAAGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAGTCGAGGCTGCCTTTCAGCTTGCCTTCGATTCCGTCGTTGTCGcc                  >  2:1446461/1‑100 (MQ=255)
                                                                                  ggCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGcc                  <  2:56065/100‑1 (MQ=255)
                                                                                    ccTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAc                <  2:2322937/100‑1 (MQ=255)
                                                                                     cTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAcc               <  2:797996/100‑1 (MQ=255)
                                                                                      tGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCg              >  2:3829167/1‑100 (MQ=255)
                                                                                       gCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCTGCTTGCCTTCGATTTCGTCGTTGTAGCCACCgg             <  2:2602667/100‑1 (MQ=255)
                                                                                       gCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCgg             >  1:3783136/1‑100 (MQ=255)
                                                                                       gCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCgg             <  2:2997258/100‑1 (MQ=255)
                                                                                         ccTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTc           >  1:3403903/1‑100 (MQ=255)
                                                                                          cTAAGGTAAGACGCCAGACGTTACCAGTATTGGTTTGCTTTGGATGTATCGAGCATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCATCGGTCt          <  2:2187650/100‑1 (MQ=255)
                                                                                           tAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTg         <  1:1723062/100‑1 (MQ=255)
                                                                                           tAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTg         >  1:2640773/1‑100 (MQ=255)
                                                                                           tAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTg         <  1:1541561/100‑1 (MQ=255)
                                                                                              ggCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGtt      <  1:309506/100‑1 (MQ=255)
                                                                                               gCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTg     <  1:241396/100‑1 (MQ=255)
                                                                                                cAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTg     >  2:536475/1‑99 (MQ=255)
                                                                                                cAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGa    >  1:1332475/1‑100 (MQ=255)
                                                                                                cAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGa    >  1:3940924/1‑100 (MQ=255)
                                                                                                cAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGa    >  2:929522/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATC  >  NC_000913/3918206‑3918403

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: