Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,932,296 G→T 10.8% D131Y (GAT→TAT)  pykA → pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,932,2960GT10.8% 659.3 / 36.8 193D131Y (GAT→TAT) pykApyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base G (82/90);  new base T (9/12);  total (91/102)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 193 total aligned reads displayed.

AAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACAAA  >  CP009273/1932150‑1932442
                                                                                                                                                  |                                                                                                                                                  
aagaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGa                                                                                                                                                   <  2:2933545/148‑1 (MQ=255)
aagaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGa                                                                                                                                                   <  1:1923699/148‑1 (MQ=255)
 agaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCTGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGAt                                                                                                                                                  >  2:1918889/1‑148 (MQ=255)
 agaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGAt                                                                                                                                                  <  1:2580889/148‑1 (MQ=255)
  gaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATg                                                                                                                                                 >  2:890244/1‑148 (MQ=255)
   aagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATgg                                                                                                                                                >  1:406357/1‑148 (MQ=255)
   aagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATgg                                                                                                                                                <  1:1522233/148‑1 (MQ=255)
   aagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATgg                                                                                                                                                >  1:2196213/1‑148 (MQ=255)
    agGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGt                                                                                                                                               <  2:43611/148‑1 (MQ=255)
      gCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTcg                                                                                                                                             <  1:975280/148‑1 (MQ=255)
        aaaGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTcgcg                                                                                                                                           <  1:1471978/148‑1 (MQ=255)
                 cTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTaa                                                                                                                                  >  2:2320628/1‑148 (MQ=255)
                    aaTATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAg                                                                                                                               <  1:1879857/148‑1 (MQ=255)
                     atatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGt                                                                                                                              >  1:2289772/1‑148 (MQ=255)
                     atatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGt                                                                                                                              >  1:2289622/1‑148 (MQ=255)
                      tatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTa                                                                                                                             >  1:2554373/1‑148 (MQ=255)
                      tatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTa                                                                                                                             <  1:2817092/148‑1 (MQ=255)
                       atTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTAc                                                                                                                            <  1:793385/148‑1 (MQ=255)
                         tGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTg                                                                                                                          >  2:3280376/1‑148 (MQ=255)
                         tGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTg                                                                                                                          >  2:691030/1‑148 (MQ=255)
                            ggATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGaa                                                                                                                       >  2:2674281/1‑148 (MQ=255)
                                aaaTTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTc                                                                                                                   >  2:2207499/1‑148 (MQ=255)
                                  aTTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAg                                                                                                                 >  1:3476679/1‑148 (MQ=255)
                                    tCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAggg                                                                                                               >  1:593342/1‑148 (MQ=255)
                                      cTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCa                                                                                                             <  2:2105927/148‑1 (MQ=255)
                                       tGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCAt                                                                                                            >  1:1146422/1‑148 (MQ=255)
                                        gCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATg                                                                                                           >  2:685761/1‑148 (MQ=255)
                                          tCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGaa                                                                                                         <  2:129875/148‑1 (MQ=255)
                                           cGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCCTGaaa                                                                                                        >  2:839726/1‑148 (MQ=255)
                                              cGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAgtg                                                                                                     <  2:2020294/148‑1 (MQ=255)
                                               gCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAgtgt                                                                                                    <  2:439256/148‑1 (MQ=255)
                                                  aaCCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCa                                                                                                 >  1:899314/1‑148 (MQ=255)
                                                    ccTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCAcc                                                                                               <  2:3260688/148‑1 (MQ=255)
                                                     cTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGTTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAATATTTCACCg                                                                                              >  1:498800/1‑148 (MQ=255)
                                                      tGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGa                                                                                             >  1:1845321/1‑148 (MQ=255)
                                                        ggTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAg                                                                                           <  2:1713370/148‑1 (MQ=255)
                                                        ggTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAg                                                                                           >  2:1309575/1‑148 (MQ=255)
                                                        ggTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAg                                                                                           <  1:1222388/148‑1 (MQ=255)
                                                          tAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTc                                                                                         <  1:3221889/148‑1 (MQ=255)
                                                          tAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTc                                                                                         >  2:504865/1‑148 (MQ=255)
                                                           aaaGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCa                                                                                        <  2:519888/148‑1 (MQ=255)
                                                           aaaGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCa                                                                                        >  1:246832/1‑148 (MQ=255)
                                                              ggTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCg                                                                                     >  2:1835309/1‑148 (MQ=255)
                                                              ggTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCg                                                                                     >  2:962869/1‑148 (MQ=255)
                                                               gTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCTCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGt                                                                                    <  2:464738/148‑1 (MQ=255)
                                                                   aGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggt                                                                                <  1:2189994/148‑1 (MQ=255)
                                                                        aCAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTcc                                                                           <  1:3194795/148‑1 (MQ=255)
                                                                        aCAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTcc                                                                           >  2:3483686/1‑148 (MQ=255)
                                                                         cAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTccc                                                                          <  2:3064889/148‑1 (MQ=255)
                                                                          aaaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTcccc                                                                         >  2:2591205/1‑148 (MQ=255)
                                                                           aaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCct                                                                        <  1:715305/148‑1 (MQ=255)
                                                                           aaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACAGTCGGTGGTCCCct                                                                        >  2:2133892/1‑148 (MQ=255)
                                                                             gAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCctct                                                                      >  1:2382188/1‑148 (MQ=255)
                                                                             gAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCctct                                                                      <  1:3280376/148‑1 (MQ=255)
                                                                              aaaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCtctc                                                                     >  2:1770471/1‑148 (MQ=255)
                                                                               aaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTcc                                                                    >  2:1561088/1‑148 (MQ=255)
                                                                                 aaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcaa                                                                  <  2:3453181/148‑1 (MQ=255)
                                                                                 aaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcaa                                                                  >  1:567906/1‑148 (MQ=255)
                                                                                  aGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcaac                                                                 <  2:1771131/148‑1 (MQ=255)
                                                                                   gTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcaaca                                                                <  2:1183068/148‑1 (MQ=255)
                                                                                       gTATCGACTACCAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATaa                                                            <  1:3360264/148‑1 (MQ=255)
                                                                                       gTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATaa                                                            >  1:175251/1‑148 (MQ=255)
                                                                                        tATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATaaa                                                           >  1:774783/1‑148 (MQ=255)
                                                                                         aTCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAg                                                          <  2:2600627/148‑1 (MQ=255)
                                                                                           cGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGt                                                        <  2:972954/148‑1 (MQ=255)
                                                                                               tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATca                                                    <  1:2737436/148‑1 (MQ=255)
                                                                                               tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATca                                                    <  1:547404/148‑1 (MQ=255)
                                                                                                aCAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaa                                                   >  1:2607893/1‑148 (MQ=255)
                                                                                                aCAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaa                                                   >  2:3028251/1‑148 (MQ=255)
                                                                                                 cAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaac                                                  >  2:2712471/1‑148 (MQ=255)
                                                                                                 cAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaac                                                  <  1:691014/148‑1 (MQ=255)
                                                                                                   aaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaacaa                                                <  2:1775399/148‑1 (MQ=255)
                                                                                                       cctgcctgCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACtt                                            >  2:2472530/1‑148 (MQ=255)
                                                                                                       cctgcctgCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACtt                                            >  1:2913554/1‑148 (MQ=255)
                                                                                                              gctgACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGgcggcg                                     <  1:2852213/148‑1 (MQ=255)
                                                                                                               ctgACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGcggcgg                                    <  2:1845563/148‑1 (MQ=255)
                                                                                                                tgACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGt                                   <  2:593342/148‑1 (MQ=255)
                                                                                                                  aCGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGttt                                 <  2:1039945/148‑1 (MQ=255)
                                                                                                                   cgtcgtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTg                                >  1:1009166/1‑148 (MQ=255)
                                                                                                                   cgtcgtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTg                                >  2:3316810/1‑148 (MQ=255)
                                                                                                                            tGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGaa                       <  2:350032/148‑1 (MQ=255)
                                                                                                                             ggTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAg                      <  2:856781/148‑1 (MQ=255)
                                                                                                                             ggTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTGAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAg                      <  2:2070692/148‑1 (MQ=255)
                                                                                                                               tGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcg                    <  1:3102535/148‑1 (MQ=255)
                                                                                                                               tGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcg                    <  1:3021213/148‑1 (MQ=255)
                                                                                                                               tGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcg                    <  1:1250391/148‑1 (MQ=255)
                                                                                                                               tGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcg                    <  1:2434534/148‑1 (MQ=255)
                                                                                                                                gACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcgc                   <  2:175251/148‑1 (MQ=255)
                                                                                                                                gACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcgc                   <  1:2181467/148‑1 (MQ=255)
                                                                                                                                     ccTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAAAATGGCGGCGGTTTGTCGGCGGAAGCGCTGAcc              >  1:1023228/1‑148 (MQ=255)
                                                                                                                                     ccTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCTACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGAcc              >  2:714926/1‑148 (MQ=255)
                                                                                                                                         ctgctgGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGaaa          >  1:805341/1‑148 (MQ=255)
                                                                                                                                         ctgctgGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGaaa          <  1:1770471/148‑1 (MQ=255)
                                                                                                                                           gctgGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGaaaaa        >  2:1367374/1‑148 (MQ=255)
                                                                                                                                            ctgGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAg       <  1:2033849/148‑1 (MQ=255)
                                                                                                                                            ctgGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAg       <  2:762430/148‑1 (MQ=255)
                                                                                                                                               gACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACa    >  2:2369995/1‑148 (MQ=255)
                                                                                                                                               gACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACa    <  2:2554373/148‑1 (MQ=255)
                                                                                                                                                 cGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACaaa  <  1:3462521/148‑1 (MQ=255)
                                                                                                                                                 cGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACaaa  >  2:1167761/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACAAA  >  CP009273/1932150‑1932442

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: ‑