Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,325,492 A→T 44.4% E63V (GAA→GTA)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,325,4920AT44.4% 86.2 / inf 230E63V (GAA→GTA) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base A (62/66);  new base T (53/49);  total (115/115)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 233 total aligned reads displayed.

CCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGAC  >  CP009273/1325345‑1325637
                                                                                                                                                   |                                                                                                                                                 
ccAAAACGATCAACAAGTATCTGGTTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGa                                                                                                                                                   <  2:1929915/148‑1 (MQ=255)
ccAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGa                                                                                                                                                   <  1:255015/148‑1 (MQ=255)
    aaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGt                                                                                                                                               <  1:3398745/148‑1 (MQ=255)
    aaCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGt                                                                                                                                               <  2:703492/148‑1 (MQ=255)
     aCGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTg                                                                                                                                              >  1:1245548/1‑148 (MQ=255)
          caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGct                                                                                                                                         >  1:2201232/1‑148 (MQ=255)
          caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGct                                                                                                                                         >  1:1427851/1‑148 (MQ=255)
          caacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGct                                                                                                                                         >  1:2889/1‑148 (MQ=255)
           aacaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGctc                                                                                                                                        <  1:691485/148‑1 (MQ=255)
            acaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGctct                                                                                                                                       >  1:1430337/1‑148 (MQ=255)
            acaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGctct                                                                                                                                       <  2:1123327/148‑1 (MQ=255)
            acaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGctct                                                                                                                                       >  1:2640350/1‑148 (MQ=255)
            acaaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGctct                                                                                                                                       >  1:2638609/1‑148 (MQ=255)
             caaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCtctc                                                                                                                                      >  1:1861974/1‑148 (MQ=255)
             caaGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCtctc                                                                                                                                      >  2:2544720/1‑148 (MQ=255)
               aGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGt                                                                                                                                    >  2:2718350/1‑148 (MQ=255)
                     tGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCg                                                                                                                              <  2:820125/148‑1 (MQ=255)
                        gTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTAt                                                                                                                           <  2:3057519/148‑1 (MQ=255)
                            tGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATgggg                                                                                                                       >  1:266785/1‑148 (MQ=255)
                              aCTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGt                                                                                                                     <  2:1413330/148‑1 (MQ=255)
                                tACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTCCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTg                                                                                                                   <  2:2675000/148‑1 (MQ=255)
                                     ggtgAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCg                                                                                                              >  2:1000423/1‑148 (MQ=255)
                                     ggtgAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCg                                                                                                              <  2:1340920/148‑1 (MQ=255)
                                      gtgAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGt                                                                                                             >  1:2108189/1‑148 (MQ=255)
                                        gAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTgg                                                                                                           <  1:470306/148‑1 (MQ=255)
                                            tCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGcaca                                                                                                       >  1:486746/1‑148 (MQ=255)
                                             ccAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACaa                                                                                                      >  2:616429/1‑148 (MQ=255)
                                              cAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAAt                                                                                                     <  1:199808/148‑1 (MQ=255)
                                                   tCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGa                                                                                                >  1:3469222/1‑148 (MQ=255)
                                                     ggTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACACTTGGGAgg                                                                                              >  2:3450241/1‑148 (MQ=255)
                                                      gTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGgc                                                                                             >  1:2823524/1‑148 (MQ=255)
                                                       tCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGgcg                                                                                            >  1:913001/1‑148 (MQ=255)
                                                              cgcgATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATg                                                                                     <  1:2399688/148‑1 (MQ=255)
                                                              cgcgATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATg                                                                                     <  1:1973841/148‑1 (MQ=255)
                                                                  aTTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAgt                                                                                 >  2:2324944/1‑148 (MQ=255)
                                                                  aTTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAgt                                                                                 >  1:1818150/1‑148 (MQ=255)
                                                                       ccGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGc                                                                            <  2:1587647/148‑1 (MQ=255)
                                                                          aCCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTg                                                                         >  1:200448/1‑148 (MQ=255)
                                                                              gTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTaa                                                                     <  2:1430337/148‑1 (MQ=255)
                                                                              gTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTaa                                                                     <  2:2638609/148‑1 (MQ=255)
                                                                                ggCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAag                                                                   <  2:2956885/148‑1 (MQ=255)
                                                                                  cTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAagag                                                                 <  1:509076/148‑1 (MQ=255)
                                                                                       cTGCAAAAAAGATTGCCGACTCTACCTCCACCAAGACGGATAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGgt                                                            <  2:106145/148‑1 (MQ=255)
                                                                                        tGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGgtc                                                           <  1:2744069/148‑1 (MQ=255)
                                                                                          ccAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGgtcgt                                                         >  2:910975/1‑148 (MQ=255)
                                                                                          ccAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGgtcgt                                                         <  2:2229877/148‑1 (MQ=255)
                                                                                            aaaaaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGtct                                                       <  2:2169932/148‑1 (MQ=255)
                                                                                            aaaaaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGtct                                                       <  2:2784740/148‑1 (MQ=255)
                                                                                              aaaGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTctct                                                     <  2:212648/148‑1 (MQ=255)
                                                                                               aaGAGTGCCGACTCTACCTCGACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctg                                                    >  2:2834014/1‑148 (MQ=255)
                                                                                                 gagTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaa                                                  <  2:1982403/148‑1 (MQ=255)
                                                                                                 gagTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaa                                                  <  2:3013863/148‑1 (MQ=255)
                                                                                                    tGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaactg                                               >  1:1870414/1‑148 (MQ=255)
                                                                                                     gCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaactga                                              >  2:2533126/1‑148 (MQ=255)
                                                                                                     gCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaactga                                              >  2:2244674/1‑148 (MQ=255)
                                                                                                      ccGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaactgaa                                             <  2:376113/148‑1 (MQ=255)
                                                                                                      ccGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaactgaa                                             <  2:2640739/148‑1 (MQ=255)
                                                                                                      ccGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaactgaa                                             <  2:829364/148‑1 (MQ=255)
                                                                                                      ccGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTctgaactgaa                                             <  2:860242/148‑1 (MQ=255)
                                                                                                       cGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGaaa                                            <  2:1861974/148‑1 (MQ=255)
                                                                                                       cGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGaaa                                            <  1:3038332/148‑1 (MQ=255)
                                                                                                       cGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGaaa                                            <  1:1267712/148‑1 (MQ=255)
                                                                                                          ctctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaacaa                                         >  1:1036229/1‑148 (MQ=255)
                                                                                                          ctctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaacaa                                         >  1:425481/1‑148 (MQ=255)
                                                                                                           tctACCTCCACCAAGACGGCTAAAAAGCCTAAACAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaacaac                                        <  2:3409492/148‑1 (MQ=255)
                                                                                                           tctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAaacaac                                        <  2:1242782/148‑1 (MQ=255)
                                                                                                            ctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACt                                       <  1:1482825/148‑1 (MQ=255)
                                                                                                            ctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACt                                       <  1:3249857/148‑1 (MQ=255)
                                                                                                            ctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGACGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACt                                       >  2:1528658/1‑148 (MQ=255)
                                                                                                            ctACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACt                                       >  1:3004092/1‑148 (MQ=255)
                                                                                                              aCCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGACGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTgg                                     >  2:2418682/1‑148 (MQ=255)
                                                                                                               ccTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGc                                    >  2:719592/1‑148 (MQ=255)
                                                                                                                 tCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTg                                  >  2:2856909/1‑148 (MQ=255)
                                                                                                                 tCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTg                                  >  2:1066647/1‑148 (MQ=255)
                                                                                                                   caccaAGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGaa                                >  1:1112759/1‑148 (MQ=255)
                                                                                                                    accaAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGaaa                               <  1:2117384/148‑1 (MQ=255)
                                                                                                                     ccaAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGaaaa                              <  2:348437/148‑1 (MQ=255)
                                                                                                                       aaGACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAg                            <  2:3158862/148‑1 (MQ=255)
                                                                                                                       aaGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAg                            >  1:3221107/1‑148 (MQ=255)
                                                                                                                        aGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGc                           >  2:3143400/1‑148 (MQ=255)
                                                                                                                         taCGGCTAATAAGACTAAATAGGATGTACGTGGCGCTCTCGTCACCCGTTTGGGGGTTGACCGTTGGCACATTTGTGCTGCGCACTATGAAGTGTTGCCTGGTAATGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGcc                          <  2:3276428/147‑1 (MQ=255)
                                                                                                                         gACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGcc                          >  2:1289544/1‑148 (MQ=255)
                                                                                                                         gACGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGcc                          >  1:1033634/1‑148 (MQ=255)
                                                                                                                          aCGGCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCg                         >  1:2591669/1‑148 (MQ=255)
                                                                                                                             gCTAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTTACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGAcc                      >  1:2363299/1‑148 (MQ=255)
                                                                                                                               tAAAAAGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACcac                    <  1:3505414/148‑1 (MQ=255)
                                                                                                                               tAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACcac                    >  1:852428/1‑148 (MQ=255)
                                                                                                                                aaaaaGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACcaca                   >  2:1250628/1‑148 (MQ=255)
                                                                                                                                aaaaaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACcaca                   >  1:1204845/1‑148 (MQ=255)
                                                                                                                                 aaaaGCCTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACat                  >  2:921166/1‑148 (MQ=255)
                                                                                                                                  aaaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatc                 <  1:2835330/148‑1 (MQ=255)
                                                                                                                                  aaaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatc                 >  2:747772/1‑148 (MQ=255)
                                                                                                                                  aaaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTATCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatc                 <  2:502168/148‑1 (MQ=255)
                                                                                                                                   aaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatct                >  2:1571177/1‑148 (MQ=255)
                                                                                                                                   aaGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatct                >  1:23009/1‑148 (MQ=255)
                                                                                                                                    aGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatcta               >  1:570286/1‑148 (MQ=255)
                                                                                                                                    aGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatcta               <  1:1325864/148‑1 (MQ=255)
                                                                                                                                      ccTAAAAAGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatctatc             >  1:2113734/1‑148 (MQ=255)
                                                                                                                                      ccTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACatctatc             <  2:3173764/148‑1 (MQ=255)
                                                                                                                                           aaaGGATGTACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCa        <  2:1427851/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTATGAAGTGTTGCCTGGTAAAGAGAAGGTCGTCTCTGAACTGAAACAACTGGCTGAAAAAGCCGACCACATCTATCTCGCAACCGAC  >  CP009273/1325345‑1325637

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: