Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,479,733 G→T 100% D806Y (GAT→TAT)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,479,7330GT98.1% 660.6 / 3.7 208D806Y (GAT→TAT) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base T (98/106);  total (99/109)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 208 total aligned reads displayed.

GGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTGGCATC  >  CP009273/1479593‑1479877
                                                                                                                                            |                                                                                                                                                
ggCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAAt                                                                                                                                           <  2:338414/148‑1 (MQ=255)
 gCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAAtc                                                                                                                                          >  1:2038281/1‑148 (MQ=255)
    gACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTcc                                                                                                                                       >  2:2828523/1‑148 (MQ=255)
            aCGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCa                                                                                                                               <  2:1634831/148‑1 (MQ=255)
            aCGCATTCTTCCGTGAAAACATGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCa                                                                                                                               >  1:3805561/1‑148 (MQ=255)
              tcATTCTTCCGTGAAAACCTGAAACTACGGGCGGAATTAGAAGAGCAGGACCACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCATCCACTATGTAATCTCCGCTCGCCACt                                                                                                                             <  2:1998556/147‑1 (MQ=255)
              gCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACt                                                                                                                             <  2:232774/148‑1 (MQ=255)
               cATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACtt                                                                                                                            >  2:1544616/1‑148 (MQ=255)
                   cttcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGAc                                                                                                                        <  2:1493700/148‑1 (MQ=255)
                     tcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAg                                                                                                                      <  1:1727522/148‑1 (MQ=255)
                     tcCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAg                                                                                                                      <  1:3196384/148‑1 (MQ=255)
                      ccGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGc                                                                                                                     <  2:1981228/148‑1 (MQ=255)
                      ccGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATAACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGc                                                                                                                     <  1:544570/148‑1 (MQ=255)
                            aaaCCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCtggtgg                                                                                                               >  1:3420613/1‑148 (MQ=255)
                              aCCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaa                                                                                                             >  2:745134/1‑148 (MQ=255)
                               ccTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaaa                                                                                                            >  1:1362375/1‑148 (MQ=255)
                                cTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaaaa                                                                                                           >  1:1960097/1‑148 (MQ=255)
                                 tGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGaaaaa                                                                                                          <  1:1421882/148‑1 (MQ=255)
                                    aaCTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGt                                                                                                       <  2:466303/148‑1 (MQ=255)
                                    aaCTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGt                                                                                                       <  2:2282448/148‑1 (MQ=255)
                                     aCTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTc                                                                                                      <  2:3626223/148‑1 (MQ=255)
                                      cTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCa                                                                                                     >  2:1874256/1‑148 (MQ=255)
                                               gAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGaaa                                                                                            <  2:2791295/148‑1 (MQ=255)
                                                  gTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGc                                                                                         <  1:2287312/148‑1 (MQ=255)
                                                  gTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGc                                                                                         <  2:3581191/148‑1 (MQ=255)
                                                   tAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGcc                                                                                        <  2:2439923/148‑1 (MQ=255)
                                                     gaagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTg                                                                                      >  2:1930285/1‑148 (MQ=255)
                                                     gaagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTg                                                                                      <  1:1610963/148‑1 (MQ=255)
                                                      aagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGa                                                                                     >  1:2779153/1‑148 (MQ=255)
                                                      aagaGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGa                                                                                     <  1:3585610/148‑1 (MQ=255)
                                                       agagCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAt                                                                                    <  1:7442/148‑1 (MQ=255)
                                                       agagCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAt                                                                                    >  1:600243/1‑148 (MQ=255)
                                                         agCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAttt                                                                                  >  1:1899929/1‑148 (MQ=255)
                                                         agCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAttt                                                                                  >  1:1155102/1‑148 (MQ=255)
                                                         agCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGACAATTCGTCAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGAttt                                                                                  >  2:1607265/1‑148 (MQ=255)
                                                          gCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTg                                                                                 <  2:2570728/148‑1 (MQ=255)
                                                          gCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTg                                                                                 <  1:108276/148‑1 (MQ=255)
                                                          gCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTg                                                                                 <  1:471470/148‑1 (MQ=255)
                                                          gCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTg                                                                                 <  1:108131/148‑1 (MQ=255)
                                                           cTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGc                                                                                >  2:2167120/1‑148 (MQ=255)
                                                           cTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGc                                                                                <  2:795486/148‑1 (MQ=255)
                                                              gAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCa                                                                             >  2:1761764/1‑148 (MQ=255)
                                                                acacAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAAc                                                                           >  1:674483/1‑148 (MQ=255)
                                                                  acaAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACtt                                                                         >  2:3487110/1‑148 (MQ=255)
                                                                  acaAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACtt                                                                         <  2:3669287/148‑1 (MQ=255)
                                                                  acaAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACtt                                                                         >  2:3343012/1‑148 (MQ=255)
                                                                      aTCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATNTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGaa                                                                     >  1:1937822/1‑148 (MQ=255)
                                                                       tCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGaaa                                                                    >  1:2782817/1‑148 (MQ=255)
                                                                        cACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGaaaa                                                                   >  2:2784032/1‑148 (MQ=255)
                                                                         aCGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGaaaaa                                                                  <  1:2435868/148‑1 (MQ=255)
                                                                         aCGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGaaaaa                                                                  <  1:515218/148‑1 (MQ=255)
                                                                               ccgcGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATg                                                            <  2:1789741/148‑1 (MQ=255)
                                                                                cgcgATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGt                                                           <  2:2987812/148‑1 (MQ=255)
                                                                                cgcgATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGt                                                           <  2:3290421/148‑1 (MQ=255)
                                                                                   gATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGa                                                        <  2:549267/148‑1 (MQ=255)
                                                                                    atatTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGAt                                                       <  1:1624593/148‑1 (MQ=255)
                                                                                     tatTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATc                                                      >  1:3614568/1‑148 (MQ=255)
                                                                                     tatTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTAAGTAATCTCCGCTCGCCACAACGACAGCTGCTGGAAAAAAGTCAGCCGCGAAACGCCTGATATGATCAACTTTGAAAAAAGCATGTTGATc                                                      >  2:284419/1‑148 (MQ=255)
                                                                                      atTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCa                                                     >  2:3317811/1‑148 (MQ=255)
                                                                                       ttCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCaa                                                    <  2:168904/148‑1 (MQ=255)
                                                                                        tCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCaaa                                                   >  1:3266718/1‑148 (MQ=255)
                                                                                        tCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCaaa                                                   >  1:1231772/1‑148 (MQ=255)
                                                                                         cTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAag                                                  >  2:2788906/1‑148 (MQ=255)
                                                                                           ggTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAagag                                                <  2:667891/148‑1 (MQ=255)
                                                                                            gTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAgg                                               <  1:2934601/148‑1 (MQ=255)
                                                                                            gTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAgg                                               <  1:2883221/148‑1 (MQ=255)
                                                                                            gTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAgg                                               <  2:2056209/148‑1 (MQ=255)
                                                                                             ttGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAggg                                              <  2:3420613/148‑1 (MQ=255)
                                                                                             ttGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAggg                                              <  1:3018275/148‑1 (MQ=255)
                                                                                              tgatgaCGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAATCATGTTGATCAAAGAGGgc                                             >  2:2012240/1‑148 (MQ=255)
                                                                                                atgaCGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGgcgc                                           <  2:3229778/148‑1 (MQ=255)
                                                                                                   aCGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGa                                        >  1:675778/1‑148 (MQ=255)
                                                                                                      aaaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaa                                     >  1:1901063/1‑148 (MQ=255)
                                                                                                       aaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaa                                    <  1:3386888/148‑1 (MQ=255)
                                                                                                       aaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaa                                    <  1:1578007/148‑1 (MQ=255)
                                                                                                       aaCGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGaaaaa                                    <  2:1791121/148‑1 (MQ=255)
                                                                                                         cGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAAt                                  >  2:2931862/1‑148 (MQ=255)
                                                                                                            tgttTGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAg                               >  1:2564365/1‑148 (MQ=255)
                                                                                                              tttGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCa                             >  2:1747801/1‑148 (MQ=255)
                                                                                                               ttGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCaa                            >  2:64793/1‑148 (MQ=255)
                                                                                                                tGAGTTCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGAAAg                           <  1:665342/148‑1 (MQ=255)
                                                                                                                    ttCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTgg                       <  2:427263/148‑1 (MQ=255)
                                                                                                                    ttCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTgg                       >  1:2809871/1‑148 (MQ=255)
                                                                                                                     tCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGa                      <  2:2533301/148‑1 (MQ=255)
                                                                                                                     tCTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGa                      >  1:677870/1‑148 (MQ=255)
                                                                                                                      cTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGAt                     <  1:1761764/148‑1 (MQ=255)
                                                                                                                      cTACGACCAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGAt                     <  1:3654364/148‑1 (MQ=255)
                                                                                                                            ccAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCg               >  1:216698/1‑148 (MQ=255)
                                                                                                                             cAGCGCATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGa              >  1:3670365/1‑148 (MQ=255)
                                                                                                                                 gcATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACtt          <  2:1073118/148‑1 (MQ=255)
                                                                                                                                 gcATCAGACACTATGGAATCTCCGCTCGCCACTTCGACAGCTGTTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACGTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACtt          <  1:2600936/148‑1 (MQ=255)
                                                                                                                                  cATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTc         >  1:828696/1‑148 (MQ=255)
                                                                                                                                  cATCAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGGGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTc         <  2:2394546/148‑1 (MQ=255)
                                                                                                                                     cAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGa                                        >  2:1288364/1‑114 (MQ=255)
                                                                                                                                     cAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGa                                        >  2:1007075/1‑114 (MQ=255)
                                                                                                                                     cAGCCACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGa                                        <  1:1007075/114‑1 (MQ=255)
                                                                                                                                        ccACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTGGCAt   >  2:3176923/1‑148 (MQ=255)
                                                                                                                                        ccACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTGGCAt   <  1:745260/148‑1 (MQ=255)
                                                                                                                                         cACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTGGCATc  <  1:3343012/148‑1 (MQ=255)
                                                                                                                                         cACTATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTGGCATc  <  1:3487110/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
GGCAGACGCGTCACGCATTCTTCCGTGAAAACCTGAAACTACGGGCGGAAGTAGAAGAGCTGGAACACAAATCACGTCGCCGCGATATTCTGGTTGATGACGAAACGTTGTTTGAGTTCTACGACCAGCGCATCAGCCACGATGTAATCTCCGCTCGCCACTTCGACAGCTGGTGGAAAAAAGTCAGCCGCGAAACGCCTGATTTGCTCAACTTTGAAAAAAGCATGTTGATCAAAGAGGGCGCAGAAAAAATCAGCAAGCTGGATTACCCGAACTTCTGGCATC  >  CP009273/1479593‑1479877

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