Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,917,787 T→A 5.3% Y160F (TAC→TTC)  atpG ← ATP synthase F1 complex subunit gamma

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,917,7870TA 5.3% 308.4 / 3.8 94Y160F (TAC→TTC) atpGATP synthase F1 complex subunit gamma
Reads supporting (aligned to +/- strand):  ref base T (50/39);  new base A (3/2);  total (53/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGCATTGAACG  >  NC_000913/3917638‑3917936
                                                                                                                                                     |                                                                                                                                                     
gTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTcgtcgt                                                                                                                                                       <  2:483722/150‑1 (MQ=255)
gTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTcgtcgt                                                                                                                                                       <  2:721877/150‑1 (MQ=255)
  tCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAg                                                                                                                                                     <  2:71998/150‑1 (MQ=255)
     tACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGcc                                                                                                                                                  <  2:1070203/150‑1 (MQ=255)
      aCAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTACGACGAAGGccg                                                                                                                                                 <  1:1356855/150‑2 (MQ=255)
        aGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGc                                                                                                                                               <  2:776040/150‑1 (MQ=255)
         ggTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCa                                                                                                                                              >  1:203636/1‑150 (MQ=255)
         ggTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCa                                                                                                                                              <  1:654003/150‑1 (MQ=255)
           tAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAAc                                                                                                                                            <  1:360381/150‑1 (MQ=255)
            aaTCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACa                                                                                                                                           >  1:509154/1‑150 (MQ=255)
              tCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCAGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATc                                                                                                                                         >  2:529061/1‑150 (MQ=255)
                ccAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCAc                                                                                                                                       >  1:1073105/1‑150 (MQ=255)
                 cAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACt                                                                                                                                      <  1:256767/150‑1 (MQ=255)
                    gATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCAGCAACATCACtttt                                                                                                                                   >  1:385015/1‑150 (MQ=255)
                         aTGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCgg                                                                                                                              <  1:50033/150‑1 (MQ=255)
                         aTGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCgg                                                                                                                              <  2:270751/150‑1 (MQ=255)
                            tttCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGAcc                                                                                                                           <  1:1616030/150‑1 (MQ=255)
                             ttCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCg                                                                                                                          <  2:1362290/150‑1 (MQ=255)
                             ttCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCAGCAACATCACTTTTACCGGACCg                                                                                                                          <  2:385015/150‑1 (MQ=255)
                              tCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGAGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGa                                                                                                                         <  2:121822/150‑1 (MQ=255)
                               cAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGAt                                                                                                                        >  1:980506/1‑150 (MQ=255)
                               cAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGAt                                                                                                                        >  2:242030/1‑150 (MQ=255)
                                aGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATc                                                                                                                       >  1:1238840/1‑150 (MQ=255)
                                aGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATc                                                                                                                       >  2:1501555/1‑150 (MQ=255)
                                         tcatcTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGAc                                                                                                              <  2:395540/150‑1 (MQ=255)
                                              tGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGa                                                                                                         <  2:1478340/150‑1 (MQ=255)
                                                aTGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAg                                                                                                       >  2:1284700/1‑150 (MQ=255)
                                                aTGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAg                                                                                                       >  2:1290599/1‑150 (MQ=255)
                                                aTGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCGTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAg                                                                                                       <  2:1368700/150‑1 (MQ=255)
                                                 tGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAgg                                                                                                      <  1:192045/150‑1 (MQ=255)
                                                 tGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAgg                                                                                                      <  2:1053896/150‑1 (MQ=255)
                                                  gCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAggg                                                                                                     <  1:293485/150‑1 (MQ=255)
                                                   ccGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGt                                                                                                    <  2:932638/150‑1 (MQ=255)
                                                    cGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGtt                                                                                                   >  1:1148983/1‑150 (MQ=255)
                                                        aaCGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATcc                                                                                               >  2:701292/1‑150 (MQ=255)
                                                        aaCGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATcc                                                                                               >  2:924991/1‑150 (MQ=255)
                                                        aaCGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATcc                                                                                               <  1:1010376/150‑1 (MQ=255)
                                                          cGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGGCTGCAACATCACTTTTACCGGACCGATCAGTTCGTACAGGGAAGTTTTATcccc                                                                                             <  2:1073105/150‑1 (MQ=255)
                                                          cGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATcccc                                                                                             >  1:593181/1‑150 (MQ=255)
                                                               gcagcTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATccc                                                                                              >  1:1187464/1‑144 (MQ=255)
                                                               gcagcTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATccc                                                                                              <  2:1187464/144‑1 (MQ=255)
                                                                  gctggctgATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCgg                                                                                     >  2:258755/1‑150 (MQ=255)
                                                                  gctggctgATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCgg                                                                                     <  1:994886/150‑1 (MQ=255)
                                                                   ctggctgATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGt                                                                                    >  2:1363156/1‑150 (MQ=255)
                                                                     ggctgATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGa                                                                                  >  1:696191/1‑150 (MQ=255)
                                                                        tgATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCt                                                                               >  2:762103/1‑150 (MQ=255)
                                                                        tgATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCt                                                                               >  1:581130/1‑150 (MQ=255)
                                                                         gATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTg                                                                              <  1:662938/150‑1 (MQ=255)
                                                                          aTGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTgg                                                                             <  1:324182/150‑1 (MQ=255)
                                                                            ggTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGc                                                                           <  2:980506/150‑1 (MQ=255)
                                                                             gTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGca                                                                          >  2:1450012/1‑150 (MQ=255)
                                                                             gTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGca                                                                          >  2:1444217/1‑150 (MQ=255)
                                                                              tCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGcaa                                                                         >  2:255330/1‑150 (MQ=255)
                                                                               cGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGcaac                                                                        >  1:1348910/1‑150 (MQ=255)
                                                                               cGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGcaac                                                                        >  1:1345240/1‑150 (MQ=255)
                                                                                ggAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGcaaca                                                                       >  2:105744/1‑150 (MQ=255)
                                                                                        agaCATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTgcc                                                               >  1:577181/1‑150 (MQ=255)
                                                                                             tGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCa                                                          <  1:1075764/150‑1 (MQ=255)
                                                                                                tgtTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACgg                                                       >  1:365033/1‑150 (MQ=255)
                                                                                                 gtTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGa                                                      >  1:1034053/1‑150 (MQ=255)
                                                                                                   taataaATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCACCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGt                                                    >  1:1006937/1‑150 (MQ=255)
                                                                                                     ataaATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTg                                                  <  1:701292/150‑1 (MQ=255)
                                                                                                        aaTTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGTAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTgaag                                               <  2:1563857/150‑1 (MQ=255)
                                                                                                        aaTTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTgaag                                               >  2:668161/1‑150 (MQ=255)
                                                                                                        aaTTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTCCCCTGGGCAACAACATTGCCGCCCACGGAGTTgaag                                               >  2:401777/1‑150 (MQ=255)
                                                                                                         aTTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTgaaga                                              >  1:337776/1‑150 (MQ=255)
                                                                                                             gttgCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGa                                          >  2:1493578/1‑150 (MQ=255)
                                                                                                              ttgCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACggt                                         <  1:668161/150‑3 (MQ=255)
                                                                                                               tgCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGaca                                        >  1:73709/1‑150 (MQ=255)
                                                                                                                 cTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACg                                      >  1:616203/1‑150 (MQ=255)
                                                                                                                  tGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGc                                     >  1:144627/1‑150 (MQ=255)
                                                                                                                     cAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCtt                                  >  2:1502235/1‑150 (MQ=255)
                                                                                                                     cAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCtt                                  >  2:204558/1‑150 (MQ=255)
                                                                                                                        tGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCg                               >  1:1169213/1‑150 (MQ=255)
                                                                                                                         gTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGa                              >  1:1004158/1‑150 (MQ=255)
                                                                                                                          tAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCTCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAg                             <  1:1501555/150‑1 (MQ=255)
                                                                                                                            aaGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGcc                           <  2:221009/150‑1 (MQ=255)
                                                                                                                            aaGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGcc                           <  2:509154/150‑1 (MQ=255)
                                                                                                                              gCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGa                         >  1:955501/1‑150 (MQ=255)
                                                                                                                              gCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGa                         <  1:1450012/150‑1 (MQ=255)
                                                                                                                              gCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGa                         <  1:1444217/150‑1 (MQ=255)
                                                                                                                                 tGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGCGCCGATCa                      >  1:1323276/1‑150 (MQ=255)
                                                                                                                                     cAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGc                  >  1:379724/1‑150 (MQ=255)
                                                                                                                                     cAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGc                  <  1:762103/150‑1 (MQ=255)
                                                                                                                                     cAGACGGCCTTCGACGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCTGACAGGGAAGCGTTATCCCACATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGc                  >  2:588640/1‑150 (MQ=255)
                                                                                                                                      aGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCg                 >  2:429387/1‑150 (MQ=255)
                                                                                                                                       gcaggCCTACGACGAAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTgg                                                                             <  1:662738/86‑1 (MQ=255)
                                                                                                                                       gcaggCCTACGACGAAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTgg                                                                             >  2:662738/4‑89 (MQ=255)
                                                                                                                                       gACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGa                >  2:588807/1‑150 (MQ=255)
                                                                                                                                        caggCCTACGACGAAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGCCACGCCTTTCGAGCCGCTCATTGCGAg               >  2:1547378/3‑150 (MQ=255)
                                                                                                                                         cGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAgg              >  1:1471216/1‑150 (MQ=255)
                                                                                                                                         cGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAgg              <  1:529061/150‑1 (MQ=255)
                                                                                                                                             cTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGc          >  2:221380/1‑150 (MQ=255)
                                                                                                                                             cTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGc          <  1:941277/150‑1 (MQ=255)
                                                                                                                                               tcgtcgTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGCAt        <  2:1034053/150‑1 (MQ=255)
                                                                                                                                                 gtcgtAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGCATTg      <  1:1397647/150‑1 (MQ=255)
                                                                                                                                                  tcgtAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGCATTGa     >  2:1329571/1‑150 (MQ=255)
                                                                                                                                                   cgtAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGCATTGaa    <  1:1363156/150‑1 (MQ=255)
                                                                                                                                                     tAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGCATTGAACg  <  2:577181/150‑1 (MQ=255)
                                                                                                                                                     tAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAACGGTTTTCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGCCACGCCTTTCGAGCCGATCATTGCGAGGTCGCATGGCACg  >  1:1468602/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGCATTGAACG  >  NC_000913/3917638‑3917936

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