Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,918,291 C→T 15.4% intergenic (‑26/+25) 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,2910CT15.4% 429.1 / 55.3 174intergenic (‑26/+25)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 C (75/72);  new base T (15/12);  total (90/84)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 182 total aligned reads displayed.

TCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACC  >  NC_000913/3918192‑3918390
                                                                                                   |                                                                                                   
tcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGcc                                                                                                     <  1:2539554/100‑1 (MQ=255)
tcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGcc                                                                                                     >  1:1671887/1‑100 (MQ=255)
ctcATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGcc                                                                                                     <  1:1730952/98‑1 (MQ=255)
  cATCTCTTTAGTGACTGTGTGCGTGTTCTGGTCGCCTGCTATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTTAATGCCTTGCGGCCTg                                                                                                   <  2:1567561/100‑1 (MQ=255)
  cATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTg                                                                                                   >  1:798866/1‑100 (MQ=255)
  cATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCACCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTg                                                                                                   >  1:3473704/1‑100 (MQ=255)
   aTCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGc                                                                                                  <  1:1618685/100‑1 (MQ=255)
    tCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGcc                                                                                                 >  1:1357478/1‑100 (MQ=255)
     cGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGcc                                                                                                 <  1:508410/99‑1 (MQ=255)
       cTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTa                                                                                              <  2:3554115/100‑1 (MQ=255)
       cTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTa                                                                                              <  1:481946/100‑1 (MQ=255)
         ttAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAg                                                                                            <  1:267729/100‑1 (MQ=255)
            gTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCa                                                                                         >  1:3442058/1‑100 (MQ=255)
            gTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCa                                                                                         <  2:2828522/100‑1 (MQ=255)
            gTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCa                                                                                         <  2:2968678/100‑1 (MQ=255)
             tGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCaa                                                                                        >  1:3378594/1‑100 (MQ=255)
             tGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCaa                                                                                        >  1:2702475/1‑100 (MQ=255)
              gATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAg                                                                                       <  2:2267034/100‑1 (MQ=255)
              gATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAg                                                                                       >  1:941888/1‑100 (MQ=255)
              gATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAg                                                                                       >  2:1714699/1‑100 (MQ=255)
                tCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAgcc                                                                                     <  2:1621692/100‑1 (MQ=255)
                tCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAgcc                                                                                     >  2:418194/1‑100 (MQ=255)
                 cTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGgccg                                                                                    <  1:2766324/100‑1 (MQ=255)
                 cTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAgcc                                                                                     >  2:2880216/1‑99 (MQ=255)
                  ttttGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGgccgc                                                                                   <  2:1560310/100‑1 (MQ=255)
                  ttttGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAgccgc                                                                                   <  1:733913/100‑1 (MQ=255)
                     tGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAg                                                                                <  2:1869214/100‑1 (MQ=38)
                     tGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAg                                                                                <  1:2358135/100‑1 (MQ=38)
                     tGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTCTTACGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAg                                                                                <  2:2562704/100‑1 (MQ=37)
                      gCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAgg                                                                               <  1:2602092/100‑2 (MQ=38)
                      gCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGa                                                                               >  2:3546227/1‑100 (MQ=255)
                      gCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGa                                                                               >  1:2577454/1‑100 (MQ=255)
                      gCGTGTTCTGGACGCTGGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAgg                                                                               <  1:893520/100‑2 (MQ=37)
                        gtgtTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACg                                                                             >  1:1152723/1‑100 (MQ=255)
                        gtgtTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACg                                                                             >  1:3604936/1‑100 (MQ=255)
                        gtgtTCTGGACGCTTGCGATCTTACGACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAg                                                                                <  2:3400613/97‑1 (MQ=37)
                         tgtTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGgccgcaaggcat                                                                            <  2:3602508/100‑8 (MQ=37)
                          gtTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGtt                                                                           >  1:929799/1‑100 (MQ=255)
                             cTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTAcc                                                                        >  2:726065/1‑100 (MQ=255)
                              tGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCa                                                                       >  2:492972/1‑100 (MQ=255)
                              tGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCAGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTAcc                                                                        >  1:2500345/1‑99 (MQ=255)
                                  cGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGAt                                                                   >  2:2665251/1‑100 (MQ=255)
                                     ttGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTgg                                                                >  1:3104715/1‑100 (MQ=255)
                                      tGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGTGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTggg                                                               <  1:2469555/100‑1 (MQ=255)
                                      tGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTggg                                                               >  1:3459588/1‑100 (MQ=255)
                                      ngCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTggg                                                               >  1:636097/2‑100 (MQ=255)
                                       gCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGGCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGt                                                              <  2:603162/100‑1 (MQ=255)
                                       gCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGt                                                              <  2:145905/100‑1 (MQ=255)
                                            cTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCtt                                                         <  1:699695/100‑1 (MQ=255)
                                                 tacGTATCGCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACTTTACCAGGATTGGGTTGCTTTGAAg                                                    <  2:1319853/100‑1 (MQ=255)
                                                  acGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGa                                                  >  2:2551136/1‑100 (MQ=255)
                                                   cGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGa                                                  <  2:947357/100‑1 (MQ=255)
                                                     tATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAAt                                                >  2:1104036/1‑100 (MQ=255)
                                                     tATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAAt                                                <  1:544394/100‑1 (MQ=255)
                                                      aTCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTATCAGGATTGGGTTGCTTTGAAGGAATc                                               <  1:2216420/100‑1 (MQ=255)
                                                        ctctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGa                                             >  2:2937514/1‑100 (MQ=255)
                                                          ctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAgg                                           >  2:1234055/1‑100 (MQ=255)
                                                          ctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAgg                                           >  1:769812/1‑100 (MQ=255)
                                                          ctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAgg                                           <  1:3521680/100‑1 (MQ=255)
                                                            tttGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGAt                                         <  1:2481616/100‑1 (MQ=255)
                                                                cgcCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCt                                     >  2:2370971/1‑100 (MQ=255)
                                                                cgcCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCt                                     >  2:3032197/1‑100 (MQ=255)
                                                                 gcCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTGACAGGATTGGGTTGCTTTGAAGGAATAGAGGCTg                                        >  2:3215260/1‑96 (MQ=255)
                                                                   cGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTc                                  <  1:1615365/100‑1 (MQ=255)
                                                                    ggCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCa                                 <  1:1085593/100‑1 (MQ=255)
                                                                      ccATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGc                               <  1:314612/100‑1 (MQ=255)
                                                                      ccATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGc                               >  2:876329/1‑100 (MQ=255)
                                                                         tGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTg                            <  1:1642373/100‑1 (MQ=255)
                                                                         tGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTg                            <  1:3091167/100‑1 (MQ=255)
                                                                           aGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGcc                          >  1:427006/1‑100 (MQ=255)
                                                                            tcTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCt                         >  2:211134/2‑100 (MQ=255)
                                                                             cTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCtt                        >  2:3337242/1‑100 (MQ=255)
                                                                                ctcctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGa                     <  2:2034344/100‑1 (MQ=255)
                                                                                 tcctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAt                    <  1:2828714/100‑1 (MQ=255)
                                                                                  cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAtt                   >  1:500084/1‑100 (MQ=255)
                                                                                  cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAtt                   <  2:3287072/100‑1 (MQ=255)
                                                                                   ctcAATGCCTTGCGGCCTGGCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAttt                  <  2:1197540/100‑1 (MQ=255)
                                                                                   ctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAttt                  <  1:3128723/100‑1 (MQ=255)
                                                                                     cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcg                >  1:2955731/1‑100 (MQ=255)
                                                                                     cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcg                <  1:1106334/100‑1 (MQ=255)
                                                                                     cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcg                <  2:473984/100‑1 (MQ=255)
                                                                                      aaTGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcg                >  2:1684623/1‑99 (MQ=255)
                                                                                      aaTGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtc                 <  1:421647/98‑1 (MQ=255)
                                                                                        tGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcgtcg             <  2:2597140/100‑1 (MQ=255)
                                                                                        tGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcgtcg             <  2:24014/100‑1 (MQ=255)
                                                                                        cgtctggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcgtcg             >  2:563958/7‑100 (MQ=37)
                                                                                          tctggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGtt           >  1:2390542/5‑100 (MQ=38)
                                                                                          ccTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGtt           >  2:3141898/1‑100 (MQ=255)
                                                                                           ctggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTg          >  2:2116646/4‑100 (MQ=38)
                                                                                             agCGGCTTACCTTAGGTCAAGCCGCCAGACATTACCAGGAGTGGATTGCCTTGAAGGAATCGAGGATGCAGCTCACCTTGCGATCACTTTCCTAGTTGTa        >  2:3175076/2‑100 (MQ=255)
                                                                                              gCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGGAGGAATCGGGGATGCCTTTCAGCGTGCCTTCGATTTCGTCGTTGTa        >  2:2445454/1‑99 (MQ=255)
                                                                                              gCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAg       >  1:2054800/1‑100 (MQ=255)
                                                                                              gCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAg       >  1:2238961/1‑100 (MQ=255)
                                                                                               cGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGc      >  1:156475/1‑100 (MQ=255)
                                                                                               cGGCTTGCCTTAGGGCAGGCAGGCAGACGTTACTAAGATTGGGTTGCTTTGAAGGAACCGAGGATGCCTTTCACCTTGCCTTCGATTTCGCCGTTGTa        >  2:2987574/1‑98 (MQ=255)
                                                                                               cGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGc      >  2:788408/1‑100 (MQ=255)
                                                                                                ggCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGcc     >  1:3592817/1‑100 (MQ=255)
                                                                                                ggCCCGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGcc     >  2:1403723/1‑100 (MQ=255)
                                                                                                 gCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTACCTTCGATTTCGTCGTTGTAGCCa    >  2:1812421/1‑100 (MQ=255)
                                                                                                   cTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAcc  <  1:2439643/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACC  >  NC_000913/3918192‑3918390

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