Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,128,128 T→A 90.3% 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,1280TA90.3% 607.1 / 44.1 202K336* (AAA→TAA) wzccolanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (7/13);  new base A (72/110);  total (79/123)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

TGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTT  >  CP009273/2127983‑2128272
                                                                                                                                                 |                                                                                                                                                
 gctgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGgt                                                                                                                                               >  1:622206/1‑148 (MQ=255)
  ctgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGgtg                                                                                                                                              <  1:3031548/148‑1 (MQ=255)
   tgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGgtgt                                                                                                                                             <  1:2551523/148‑1 (MQ=255)
     aCCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTAc                                                                                                                                           >  2:164601/1‑148 (MQ=255)
     aCCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTAc                                                                                                                                           <  2:2017023/148‑1 (MQ=255)
     aCCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTAc                                                                                                                                           <  2:819740/148‑1 (MQ=255)
     aCCAGACTCGACAACGCGGGTCAGACCCACAATCACCAGCTGGGTTTTCGGCATCGCCGTTACGCCACCGTTAAGTATGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTATGTGTAc                                                                                                                                           >  2:649926/1‑148 (MQ=255)
       cAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAg                                                                                                                                         <  2:2649951/148‑1 (MQ=255)
       cAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAg                                                                                                                                         <  2:1776237/148‑1 (MQ=255)
        aGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGc                                                                                                                                        >  1:1266400/1‑148 (MQ=255)
         gACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCt                                                                                                                                       <  1:291423/148‑1 (MQ=255)
          aCTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCtt                                                                                                                                      >  1:3197113/1‑148 (MQ=255)
           cTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTg                                                                                                                                     >  2:1653473/1‑148 (MQ=255)
            tCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTgg                                                                                                                                    >  1:1638470/1‑148 (MQ=255)
            tCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTgg                                                                                                                                    <  1:706944/148‑1 (MQ=255)
            tCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTgg                                                                                                                                    >  2:706944/1‑148 (MQ=255)
              gACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGgag                                                                                                                                  <  2:3104442/148‑1 (MQ=255)
               aCATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGgaga                                                                                                                                 <  1:1293637/148‑1 (MQ=255)
                  tCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAttt                                                                                                                              <  2:1213393/148‑1 (MQ=255)
                      gggTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGc                                                                                                                          <  1:1407011/148‑1 (MQ=255)
                          cAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCctct                                                                                                                      <  1:19846/148‑1 (MQ=255)
                            gACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCttt                                                                                                                    >  1:2988197/1‑148 (MQ=255)
                               ggACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTaaa                                                                                                                 <  2:444717/148‑1 (MQ=255)
                               ggACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTaaa                                                                                                                 <  2:3255357/148‑1 (MQ=255)
                               ggACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTaaa                                                                                                                 <  2:2109889/148‑1 (MQ=255)
                                gACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAg                                                                                                                >  2:1028268/1‑148 (MQ=255)
                                gACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAg                                                                                                                <  1:2750762/148‑1 (MQ=255)
                                 aCAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAgg                                                                                                               >  2:2867930/1‑148 (MQ=255)
                                 aCAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAgg                                                                                                               >  2:933837/1‑148 (MQ=255)
                                 aCAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAgg                                                                                                               >  1:2450862/1‑148 (MQ=255)
                                 aCAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTAGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAgg                                                                                                               >  1:1121477/1‑148 (MQ=255)
                                  cAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGt                                                                                                              <  1:3210302/148‑1 (MQ=255)
                                   aaTCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTc                                                                                                             <  2:1858354/148‑1 (MQ=255)
                                     tctcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAg                                                                                                           >  2:128689/1‑148 (MQ=255)
                                     tctcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAg                                                                                                           <  1:1342497/148‑1 (MQ=255)
                                     tctcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAg                                                                                                           >  1:692352/1‑148 (MQ=255)
                                      ctcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgt                                                                                                          <  2:1507223/148‑1 (MQ=255)
                                      ctcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgt                                                                                                          >  1:2678628/1‑148 (MQ=255)
                                      ctcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgt                                                                                                          <  1:1226624/148‑1 (MQ=255)
                                       tcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgtt                                                                                                         <  2:2250727/148‑1 (MQ=255)
                                         ctgctgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttcg                                                                                                       <  1:1472556/148‑1 (MQ=255)
                                         ctgctgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttcg                                                                                                       <  1:3137440/148‑1 (MQ=255)
                                           gctgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttcgtt                                                                                                     >  1:331854/1‑148 (MQ=255)
                                             tgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCa                                                                                                   >  1:2649518/1‑148 (MQ=255)
                                               ggTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAAc                                                                                                 <  2:2664651/148‑1 (MQ=255)
                                                        cATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCg                                                                                        >  2:480795/1‑148 (MQ=255)
                                                              cGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTc                                                                                  <  1:2316984/148‑1 (MQ=255)
                                                                ttACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCAc                                                                                >  1:151987/1‑148 (MQ=255)
                                                                 tACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCAcc                                                                               >  1:3192415/1‑148 (MQ=255)
                                                                 tACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCAcc                                                                               <  1:2024432/148‑1 (MQ=255)
                                                                 tACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCAcc                                                                               <  1:164601/148‑1 (MQ=255)
                                                                  aCGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCa                                                                              <  2:2701736/148‑1 (MQ=255)
                                                                            tAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAg                                                                    >  2:1601089/1‑148 (MQ=255)
                                                                                tttGGCTTTTTCGTCTTCCAGCGCCTGAAGTTTATACAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCGTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCAcc                                                                >  2:611303/1‑148 (MQ=255)
                                                                                 ttGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCg                                                               >  1:472943/1‑148 (MQ=255)
                                                                                 ttGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCg                                                               >  1:2136804/1‑148 (MQ=255)
                                                                                   ggCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCt                                                             <  1:2771721/148‑1 (MQ=255)
                                                                                    gCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCtt                                                            >  2:3298529/1‑148 (MQ=255)
                                                                                    gCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCtt                                                            <  2:2578375/148‑1 (MQ=255)
                                                                                    gCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCtt                                                            >  1:237976/1‑148 (MQ=255)
                                                                                     cTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCttt                                                           <  1:1659848/148‑1 (MQ=255)
                                                                                        tttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGc                                                        >  1:1019567/1‑148 (MQ=255)
                                                                                        tttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGc                                                        <  1:460340/148‑1 (MQ=255)
                                                                                         ttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCt                                                       >  2:1951069/1‑148 (MQ=255)
                                                                                              cTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAg                                                  >  2:1268746/1‑148 (MQ=255)
                                                                                              cTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAg                                                  <  2:1121472/148‑1 (MQ=255)
                                                                                                tCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCg                                                >  2:1568459/1‑148 (MQ=255)
                                                                                                tCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCg                                                >  1:3080592/1‑148 (MQ=255)
                                                                                                tCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCg                                                >  2:2568493/1‑148 (MQ=255)
                                                                                                     cgcCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGa                                           >  2:787780/1‑148 (MQ=255)
                                                                                                     cgcCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGa                                           >  2:2442502/1‑148 (MQ=255)
                                                                                                       ccTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATc                                         >  2:1958201/1‑148 (MQ=255)
                                                                                                          gACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAAc                                      <  2:3025680/148‑1 (MQ=255)
                                                                                                               ttCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTTGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAAt                                 >  1:2831855/1‑148 (MQ=255)
                                                                                                               ttCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAAt                                 <  2:1515002/148‑1 (MQ=255)
                                                                                                                tctcCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATc                                <  1:2081224/148‑1 (MQ=255)
                                                                                                                  tcCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCtt                              <  1:1571972/148‑1 (MQ=255)
                                                                                                                   ccAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCttt                             >  1:3206503/1‑148 (MQ=255)
                                                                                                                   ccAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCttt                             <  2:502722/148‑1 (MQ=255)
                                                                                                                     agcagcGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTAt                           <  1:933837/148‑1 (MQ=255)
                                                                                                                      gcagcGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATc                          <  1:1861374/148‑1 (MQ=255)
                                                                                                                       cagcGTGCCGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTACACCATCGAATCGAGCACCGCTTTAGCTTCCAGCGGCAGATCAACAGAATCTTTATcc                         >  2:2481600/1‑148 (MQ=255)
                                                                                                                          cGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGa                      >  1:603732/1‑148 (MQ=255)
                                                                                                                          cGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGa                      <  2:10568/148‑1 (MQ=255)
                                                                                                                            tGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACg                    >  1:387985/1‑148 (MQ=255)
                                                                                                                              cGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGa                  <  1:2899882/148‑1 (MQ=255)
                                                                                                                              cGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGa                  <  2:2295478/148‑1 (MQ=255)
                                                                                                                               ggTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGaa                 <  1:1386560/148‑1 (MQ=255)
                                                                                                                                   cGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa             <  2:3216200/148‑1 (MQ=255)
                                                                                                                                   cGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa             <  2:1576872/148‑1 (MQ=255)
                                                                                                                                   cGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa             >  2:1145490/1‑148 (MQ=255)
                                                                                                                                   cGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa             >  2:841297/1‑148 (MQ=255)
                                                                                                                                     ccGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCAtt           <  1:480790/148‑1 (MQ=255)
                                                                                                                                     ccGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCAtt           <  1:481897/148‑1 (MQ=255)
                                                                                                                                      cGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTc          <  1:1203103/148‑1 (MQ=255)
                                                                                                                                           gAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgttt     >  2:746659/1‑148 (MQ=255)
                                                                                                                                            aaCTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCCTTCAgtttg    >  2:2952313/1‑148 (MQ=255)
                                                                                                                                             aCTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtttgt   <  2:3263985/148‑1 (MQ=255)
                                                                                                                                             aCTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtttgt   <  2:1356037/148‑1 (MQ=255)
                                                                                                                                              cTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGGTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtttgtt  <  1:1012275/148‑1 (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: