Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,917,784 T→A 6.0% D161V (GAC→GTC)  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,7840TA 6.0% 291.6 / 4.4 93D161V (GAC→GTC) atpGATP synthase F1 complex subunit gamma
Reads supporting (aligned to +/- strand):  ref base T (49/38);  new base A (4/2);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGCATTGA  >  NC_000913/3917637‑3917933
                                                                                                                                                   |                                                                                                                                                     
ggTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGTCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTtcgtcg                                                                                                                                                     <  1:1045469/150‑1 (MQ=255)
 gTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTcgtcgt                                                                                                                                                    <  2:721877/150‑1 (MQ=255)
 gTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTcgtcgt                                                                                                                                                    <  2:483722/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                                                                                                                     >  2:242030/1‑150 (MQ=255)
                                cAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGAt                                                                                                                     >  1:980506/1‑150 (MQ=255)
                                 aGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATc                                                                                                                    >  2:1501555/1‑150 (MQ=255)
                                 aGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATc                                                                                                                    >  1:1238840/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                                                                                            <  1:1010376/150‑1 (MQ=255)
                                                         aaCGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATcc                                                                                            >  2:701292/1‑150 (MQ=255)
                                                         aaCGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATcc                                                                                            >  2:924991/1‑150 (MQ=255)
                                                           cGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGGCTGCAACATCACTTTTACCGGACCGATCAGTTCGTACAGGGAAGTTTTATcccc                                                                                          <  2:1073105/150‑1 (MQ=255)
                                                           cGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATcccc                                                                                          >  1:593181/1‑150 (MQ=255)
                                                                gcagcTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATccc                                                                                           <  2:1187464/144‑1 (MQ=255)
                                                                gcagcTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATccc                                                                                           >  1:1187464/1‑144 (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:1345240/1‑150 (MQ=255)
                                                                                cGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGcaac                                                                     >  1:1348910/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:1450012/150‑1 (MQ=255)
                                                                                                                               gCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGa                      <  1:1444217/150‑1 (MQ=255)
                                                                                                                               gCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGa                      >  1:955501/1‑150 (MQ=255)
                                                                                                                                  tGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGCGCCGATCa                   >  1:1323276/1‑150 (MQ=255)
                                                                                                                                      ttgcaggCCTACGACGAAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGc               >  2:990549/6‑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                                                                          >  2:662738/4‑89 (MQ=255)
                                                                                                                                        gcaggCCTACGACGAAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTgg                                                                          <  1:662738/86‑1 (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)
                                                                                                                                                   |                                                                                                                                                     
GGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGCATTGA  >  NC_000913/3917637‑3917933

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: