Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,128,128 T→A 100% K336* (AAA→TAA)  wzc ← colanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,128,1280TA95.2% 725.8 / 19.8 229K336* (AAA→TAA) wzccolanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (3/8);  new base A (95/123);  total (98/131)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 229 total aligned reads displayed.

TGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTT  >  CP009273/2127983‑2128272
                                                                                                                                                 |                                                                                                                                                
tgctgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAgg                                                                                                                                                <  2:2679555/148‑1 (MQ=255)
 ggtgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGgt                                                                                                                                               >  2:1954413/3‑148 (MQ=255)
  ctgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGgtg                                                                                                                                              <  2:2341388/148‑1 (MQ=255)
    gACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTa                                                                                                                                            <  2:2900047/148‑1 (MQ=255)
    gACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTa                                                                                                                                            <  2:1901272/148‑1 (MQ=255)
    gACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTa                                                                                                                                            <  1:1398136/148‑1 (MQ=255)
     aCCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTAc                                                                                                                                           <  2:394336/148‑1 (MQ=255)
     aCCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTAc                                                                                                                                           <  1:267419/148‑1 (MQ=255)
     aCCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTAc                                                                                                                                           <  1:1728355/148‑1 (MQ=255)
     aCCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTAc                                                                                                                                           >  1:820374/1‑148 (MQ=255)
       cAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAg                                                                                                                                         >  1:1359951/1‑148 (MQ=255)
       cAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAg                                                                                                                                         <  1:1973682/148‑1 (MQ=255)
         gACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCt                                                                                                                                       >  1:1280226/1‑148 (MQ=255)
         gACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCt                                                                                                                                       <  2:3354134/148‑1 (MQ=255)
         gACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCt                                                                                                                                       >  1:793396/1‑148 (MQ=255)
          aCTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCtt                                                                                                                                      >  1:762186/1‑148 (MQ=255)
           cTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTg                                                                                                                                     <  1:1378213/148‑1 (MQ=255)
           cTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTg                                                                                                                                     <  1:880533/148‑1 (MQ=255)
           cTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTg                                                                                                                                     <  1:3147191/148‑1 (MQ=255)
           cTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTg                                                                                                                                     <  1:3147117/148‑1 (MQ=255)
              gACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGgag                                                                                                                                  <  1:870869/148‑1 (MQ=255)
              gACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGgag                                                                                                                                  <  1:3470696/148‑1 (MQ=255)
              gACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGgag                                                                                                                                  <  1:107074/148‑1 (MQ=255)
                cATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGAt                                                                                                                                <  1:3191528/148‑1 (MQ=255)
                  tCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAttt                                                                                                                              <  2:1499547/148‑1 (MQ=255)
                  tCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAttt                                                                                                                              <  2:2224243/148‑1 (MQ=255)
                  tCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAttt                                                                                                                              >  1:412543/1‑148 (MQ=255)
                  tCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAttt                                                                                                                              >  1:412729/1‑148 (MQ=255)
                     cgGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCg                                                                                                                           >  2:941508/1‑148 (MQ=255)
                      gggTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGc                                                                                                                          <  1:2706926/148‑1 (MQ=255)
                       ggTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGcc                                                                                                                         >  2:2595357/1‑148 (MQ=255)
                         tCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCctc                                                                                                                       <  1:3001591/148‑1 (MQ=255)
                             aCGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTa                                                                                                                   >  2:1836874/1‑148 (MQ=255)
                             aCGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTa                                                                                                                   >  2:2678918/1‑148 (MQ=255)
                                gACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAg                                                                                                                <  2:1978808/148‑1 (MQ=255)
                                gACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAg                                                                                                                >  2:3096663/1‑148 (MQ=255)
                                 aCAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAgg                                                                                                               >  2:3228231/1‑148 (MQ=255)
                                  cAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGt                                                                                                              <  2:1446497/148‑1 (MQ=255)
                                  cAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGt                                                                                                              <  2:966199/148‑1 (MQ=255)
                                   aaTCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTc                                                                                                             >  1:465647/1‑148 (MQ=255)
                                      ctcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgt                                                                                                          >  2:1568765/1‑148 (MQ=255)
                                        ccTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttc                                                                                                        <  1:541741/148‑1 (MQ=255)
                                         ctgctgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGTTTGGAGATTTCCGCCTCTTTAAAGGTCAgttcg                                                                                                       <  2:2173697/148‑1 (MQ=255)
                                            ctgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttcgttc                                                                                                    <  2:784363/148‑1 (MQ=255)
                                              gggTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCaa                                                                                                  <  2:871940/148‑1 (MQ=255)
                                              gggTTTTCGGCATAGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCaa                                                                                                  >  2:600240/1‑148 (MQ=255)
                                                 ttttCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTg                                                                                               >  2:175270/1‑148 (MQ=255)
                                                  tttCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTgc                                                                                              >  2:2827907/1‑148 (MQ=255)
                                                    tCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTgcgc                                                                                            <  2:824479/148‑1 (MQ=255)
                                                        cATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCCGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCg                                                                                        >  1:1029402/1‑148 (MQ=255)
                                                          tCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGAt                                                                                      <  1:1135543/148‑1 (MQ=255)
                                                          tCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCACCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTACGCCTCGTTAAAGGTCAGTTCGTTCAACAGCGCATCGAt                                                                                      >  1:3283858/1‑148 (MQ=255)
                                                             ccGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGtt                                                                                   <  1:1106014/148‑1 (MQ=255)
                                                             ccGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGtt                                                                                   <  2:1137194/148‑1 (MQ=255)
                                                             ccGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGtt                                                                                   <  2:1137508/148‑1 (MQ=255)
                                                              cGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTc                                                                                  >  1:982584/1‑148 (MQ=255)
                                                                 tACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCAcc                                                                               <  1:1682702/148‑1 (MQ=255)
                                                                 tACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCAcc                                                                               <  1:1839340/148‑1 (MQ=255)
                                                                 tACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCAcc                                                                               <  1:38773/148‑1 (MQ=255)
                                                                      gACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcga                                                                          <  1:1464438/148‑1 (MQ=255)
                                                                      gACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcga                                                                          <  2:3268156/148‑1 (MQ=255)
                                                                       aCCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcgaa                                                                         >  2:2215269/1‑148 (MQ=255)
                                                                        ccGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCATTAAAGGTCAGTTCGTTCAACTGCGCATCGAAGTTCACCatcgaat                                                                        >  1:2699423/1‑148 (MQ=255)
                                                                            tAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAg                                                                    >  2:1354466/1‑148 (MQ=255)
                                                                            tAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAg                                                                    >  1:2025336/1‑148 (MQ=255)
                                                                                   ggCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCt                                                             <  1:1906665/148‑1 (MQ=255)
                                                                                    gCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCtt                                                            >  2:2475605/1‑148 (MQ=255)
                                                                                      tttttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCtttt                                                          >  2:241414/1‑148 (MQ=255)
                                                                                       ttttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTg                                                         <  2:912111/148‑1 (MQ=255)
                                                                                       ttttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTg                                                         >  2:43916/1‑148 (MQ=255)
                                                                                         ttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCt                                                       >  1:1237985/1‑148 (MQ=255)
                                                                                         ttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCt                                                       <  2:3053175/148‑1 (MQ=255)
                                                                                         ttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCt                                                       >  2:1666050/1‑148 (MQ=255)
                                                                                          tCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCtt                                                      >  1:3289969/1‑148 (MQ=255)
                                                                                           cGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTc                                                     >  2:1112500/1‑148 (MQ=255)
                                                                                            gTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTcc                                                    <  2:491150/148‑1 (MQ=255)
                                                                                              cTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAg                                                  >  1:898900/1‑148 (MQ=255)
                                                                                              cTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAg                                                  >  1:1740099/1‑148 (MQ=255)
                                                                                                tCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCg                                                >  2:2764448/1‑148 (MQ=255)
                                                                                                    gcgcCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAg                                            >  1:511407/1‑148 (MQ=255)
                                                                                                     cgcCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTACCAGCGGCAGa                                           >  1:1360223/1‑148 (MQ=255)
                                                                                                        cTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCa                                        <  1:584419/148‑1 (MQ=255)
                                                                                                             gTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGa                                   <  2:1029402/148‑1 (MQ=255)
                                                                                                              tttCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGaa                                  <  1:1188275/148‑1 (MQ=255)
                                                                                                               ttCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAAt                                 >  2:2294835/1‑148 (MQ=255)
                                                                                                                 ctcCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCt                               <  2:302884/148‑1 (MQ=255)
                                                                                                                   ccAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCttt                             <  2:469267/148‑1 (MQ=255)
                                                                                                                       cagcGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATcc                         <  2:1830343/148‑1 (MQ=255)
                                                                                                                            tGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACg                    <  2:1016518/148‑1 (MQ=255)
                                                                                                                              cGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGa                  <  1:3063374/148‑1 (MQ=255)
                                                                                                                                 tACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAgg               >  1:3402130/1‑148 (MQ=255)
                                                                                                                                  aCGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGGTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGc              <  1:755746/148‑1 (MQ=255)
                                                                                                                                   cGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa             <  1:1533309/148‑1 (MQ=255)
                                                                                                                                   cGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa             <  1:2104150/148‑1 (MQ=255)
                                                                                                                                   cGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa             <  1:311289/148‑1 (MQ=255)
                                                                                                                                     ccGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCAtt           >  2:223047/1‑148 (MQ=255)
                                                                                                                                          tGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtt      <  1:2595357/148‑1 (MQ=255)
                                                                                                                                           gAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgttt     >  1:957188/1‑148 (MQ=255)
                                                                                                                                            aaCTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtttg    >  2:621822/1‑148 (MQ=255)
                                                                                                                                             aCTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtttgt   >  2:2505823/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTT  >  CP009273/2127983‑2128272

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: