Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,917,781 T→A 5.3% E162V (GAA→GTA)  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,7810TA 5.3% 307.5 / 5.0 95E162V (GAA→GTA) atpGATP synthase F1 complex subunit gamma
Reads supporting (aligned to +/- strand):  ref base T (50/40);  new base A (3/2);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

GATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGCAT  >  NC_000913/3917632‑3917930
                                                                                                                                                     |                                                                                                                                                     
gATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCtt                                                                                                                                                       <  1:1237909/150‑1 (MQ=255)
 aTCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTtc                                                                                                                                                      <  1:635018/150‑1 (MQ=255)
     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                                                                                                                        <  2:270751/150‑1 (MQ=255)
                               aTGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCgg                                                                                                                        <  1:50033/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:1290599/1‑150 (MQ=255)
                                                      aTGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAg                                                                                                 >  2:1284700/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:924991/1‑150 (MQ=255)
                                                              aaCGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATcc                                                                                         <  1:1010376/150‑1 (MQ=255)
                                                              aaCGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATcc                                                                                         >  2:701292/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                                                                               <  1:994886/150‑1 (MQ=255)
                                                                        gctggctgATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCgg                                                                               >  2:258755/1‑150 (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:509154/150‑1 (MQ=255)
                                                                                                                                  aaGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGcc                     <  2:221009/150‑1 (MQ=255)
                                                                                                                                    gCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGa                   >  1:955501/1‑150 (MQ=255)
                                                                                                                                    gCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGa                   <  1:1444217/150‑1 (MQ=255)
                                                                                                                                    gCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGa                   <  1:1450012/150‑1 (MQ=255)
                                                                                                                                       tGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGCGCCGATCa                >  1:1323276/1‑150 (MQ=255)
                                                                                                                                           ttgcaggCCTACGACGAAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGc            >  2:990549/6‑150 (MQ=255)
                                                                                                                                           cAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGc            <  1:762103/150‑1 (MQ=255)
                                                                                                                                           cAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGc            >  1:379724/1‑150 (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:529061/150‑1 (MQ=255)
                                                                                                                                               cGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAgg        >  1:1471216/1‑150 (MQ=255)
                                                                                                                                                   cTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGc    <  1:941277/150‑1 (MQ=255)
                                                                                                                                                   cTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGc    >  2:221380/1‑150 (MQ=255)
                                                                                                                                                     tcgtcgTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGCAt  <  2:1034053/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATCGGGTTCGTACAGGTAATCCCAGGATTTATGTTTCAGATCATCATCATCTGATGCCGGTAACGGCAGCAGCTGGCTGATGGTCGGAACCTGAGACATGGTGTTAATAAATTTGTTGCTGACAATGTAAAGCTTGTCCAGACGGCCTTCGTCGTAGGCCTGCAACATCACTTTTACCGGACCGATCAGTTCGGACAGGGAAGGGTTATCCCCCATGCCGGTGACCTGGGCAACAACATTGCCGCCCACGGAGTTGAAGAACGACACGCCTTTCGAGCCGATCATTGCGAGGTCGCAT  >  NC_000913/3917632‑3917930

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: