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,1280TA98.9% 1017.9 / ‑0.8 309K336* (AAA→TAA) wzccolanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (4/0);  new base A (122/183);  total (126/183)
Fisher's exact test for biased strand distribution p-value = 2.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 309 total aligned reads displayed.

CCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTT  >  CP009273/2127981‑2128274
                                                                                                                                                   |                                                                                                                                                  
   gctgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGgtgt                                                                                                                                               <  1:892809/150‑1 (MQ=255)
    ctgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTa                                                                                                                                              <  1:1301463/150‑1 (MQ=255)
     tgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTAc                                                                                                                                             <  1:2814365/150‑1 (MQ=255)
      gACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACa                                                                                                                                            >  2:5216506/1‑150 (MQ=255)
      gACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACa                                                                                                                                            >  2:2069686/1‑150 (MQ=255)
         cAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCt                                                                                                                                         >  2:2681910/1‑150 (MQ=255)
          aGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCtt                                                                                                                                        >  2:862585/1‑150 (MQ=255)
          aGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCtt                                                                                                                                        <  2:3008384/150‑1 (MQ=255)
           gACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTg                                                                                                                                       >  1:3050457/1‑150 (MQ=255)
             cTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCTACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGga                                                                                                                                     <  2:4113006/150‑1 (MQ=255)
              tCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGgag                                                                                                                                    <  1:5095744/150‑1 (MQ=255)
                 aCATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAtt                                                                                                                                 <  2:4546045/150‑1 (MQ=255)
                  cATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAttt                                                                                                                                >  1:559656/1‑150 (MQ=255)
                   aTCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTc                                                                                                                               >  2:4469847/1‑150 (MQ=255)
                   aTCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTc                                                                                                                               <  2:262231/150‑1 (MQ=255)
                   aTCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTc                                                                                                                               >  2:1259540/1‑150 (MQ=255)
                    tCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTcc                                                                                                                              <  2:2316179/150‑1 (MQ=255)
                    tCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTcc                                                                                                                              <  1:5196523/150‑1 (MQ=255)
                        gggTCAGACGGACAATCTCCTGCTGGGTTTTCTGCATCGACGTTACGCGACCGTTAAGATTGGCTTTTTCGTCTTCCAGCGCCTGAAGTTTCTCCAGCATCGTGCGGTACTCCGGGTGAACTTAGGTGTTCAGCTTGTATATTTCCGCct                                                                                                                          <  2:508375/150‑1 (MQ=255)
                           tCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTCAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCtt                                                                                                                       <  1:183383/150‑1 (MQ=255)
                            cAGACGGACAATCTCCTGCTGGGTGTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCttt                                                                                                                      <  2:3204154/150‑1 (MQ=255)
                             aGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTa                                                                                                                     <  2:3856948/150‑1 (MQ=255)
                             aGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTa                                                                                                                     >  2:4672517/1‑150 (MQ=255)
                             aGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTGCCAGCGCCGGACGTTTCTCCCGCAGCGTGCGGTACGCCGGGTGAACGGAGGTGTACAGCTTGGAGATTTCCGCctcggta                                                                                                                     >  2:824141/1‑146 (MQ=255)
                                  gACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGt                                                                                                                >  1:2078759/1‑150 (MQ=255)
                                   aCAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTc                                                                                                               <  1:1091008/150‑1 (MQ=255)
                                    cAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCa                                                                                                              <  1:3257988/150‑1 (MQ=255)
                                     aaTCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAg                                                                                                             >  2:2530268/1‑150 (MQ=255)
                                     aaTCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAg                                                                                                             >  1:4095518/1‑150 (MQ=255)
                                     aaTCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAg                                                                                                             >  2:592835/1‑150 (MQ=255)
                                     aaTCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAg                                                                                                             >  1:4583633/1‑150 (MQ=255)
                                        ctcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttc                                                                                                          >  2:4356258/1‑150 (MQ=255)
                                         tcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttcg                                                                                                         >  2:2781785/1‑150 (MQ=255)
                                           ctgctgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttcgtt                                                                                                       >  2:4433242/1‑150 (MQ=255)
                                           ctgctgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttcgtt                                                                                                       >  2:3254785/1‑150 (MQ=255)
                                                 ggTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTg                                                                                                 >  1:1384820/1‑150 (MQ=255)
                                                 ggTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTg                                                                                                 <  1:2366349/150‑1 (MQ=255)
                                                 ggTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTg                                                                                                 <  1:3563748/150‑1 (MQ=255)
                                                 ggTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTg                                                                                                 <  1:2628008/150‑1 (MQ=255)
                                                 ggTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTg                                                                                                 <  1:2939939/150‑1 (MQ=255)
                                                          cATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTACAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTACGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGAt                                                                                        >  1:2289827/1‑150 (MQ=255)
                                                               ccGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCa                                                                                   >  1:4654461/1‑150 (MQ=255)
                                                               ccGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCa                                                                                   <  2:4925573/150‑1 (MQ=255)
                                                                cGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCAc                                                                                  >  1:1478103/1‑150 (MQ=255)
                                                                     cgcgACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcg                                                                             <  1:5281487/150‑1 (MQ=255)
                                                                       cgACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcgaa                                                                           >  1:675365/1‑150 (MQ=255)
                                                                         aCCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcgaatc                                                                         <  2:4028613/150‑1 (MQ=255)
                                                                         aCCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcgaatc                                                                         <  1:4890292/150‑1 (MQ=255)
                                                                          ccGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcgaatcg                                                                        <  1:1843835/150‑1 (MQ=255)
                                                                            gTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAg                                                                      >  2:1070972/1‑150 (MQ=255)
                                                                            gTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAg                                                                      >  1:2601716/1‑150 (MQ=255)
                                                                             ttAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGc                                                                     >  2:749119/1‑150 (MQ=255)
                                                                             ttAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGc                                                                     <  1:5012810/150‑1 (MQ=255)
                                                                                aGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCAcc                                                                  >  2:630857/1‑150 (MQ=255)
                                                                                   ttGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCt                                                               <  2:2270926/150‑1 (MQ=255)
                                                                                    tGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCtt                                                              >  2:3682507/1‑150 (MQ=255)
                                                                                     ggCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCttt                                                             >  2:4495856/1‑150 (MQ=255)
                                                                                       cTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTg                                                           <  2:5115971/150‑1 (MQ=255)
                                                                                         ttttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCt                                                         >  1:219625/1‑150 (MQ=255)
                                                                                          tttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCtt                                                        <  1:5237747/150‑1 (MQ=255)
                                                                                          tttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCtt                                                        >  1:4791565/1‑150 (MQ=255)
                                                                                          tttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCtt                                                        >  1:1972776/1‑150 (MQ=255)
                                                                                               tCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGc                                                   <  2:1610374/150‑1 (MQ=255)
                                                                                                cTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCg                                                  >  1:743210/1‑150 (MQ=255)
                                                                                                    cAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAg                                              >  1:140666/1‑150 (MQ=255)
                                                                                                     aGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGa                                             <  2:1555719/150‑1 (MQ=255)
                                                                                                     aGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGa                                             >  1:1140185/1‑150 (MQ=255)
                                                                                                     aGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGa                                             >  1:2774044/1‑150 (MQ=255)
                                                                                                     aGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGa                                             >  2:386195/1‑150 (MQ=255)
                                                                                                     aGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGATCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGa                                             >  1:2100230/1‑150 (MQ=255)
                                                                                                          cTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAAc                                        <  1:4418131/150‑1 (MQ=255)
                                                                                                           tGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACa                                       <  2:4263421/150‑1 (MQ=255)
                                                                                                                 ttCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCt                                 <  1:1641924/150‑1 (MQ=255)
                                                                                                                  tctcCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCtt                                <  1:1609967/150‑1 (MQ=255)
                                                                                                                  tctcCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCtt                                <  2:1897980/150‑1 (MQ=255)
                                                                                                                    tcCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTa                              <  2:4098879/150‑1 (MQ=255)
                                                                                                                     ccAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTAt                             <  2:1454054/150‑1 (MQ=255)
                                                                                                                     ccAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTAt                             <  2:4729237/150‑1 (MQ=255)
                                                                                                                       agcagcGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATcc                           <  1:3472484/150‑1 (MQ=255)
                                                                                                                        gcagcGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCt                          <  1:3855460/150‑1 (MQ=255)
                                                                                                                        gcagcGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCt                          <  2:2108932/150‑1 (MQ=255)
                                                                                                                        gcagcGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCt                          <  2:1111387/150‑1 (MQ=255)
                                                                                                                            cGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACg                      >  1:4491228/1‑150 (MQ=255)
                                                                                                                              tGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGa                    >  1:4342115/1‑150 (MQ=255)
                                                                                                                                cGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAg                  <  1:3403461/150‑1 (MQ=255)
                                                                                                                                 ggTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAgg                 >  2:66343/1‑150 (MQ=255)
                                                                                                                                 ggTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAgg                 <  1:3293656/150‑1 (MQ=255)
                                                                                                                                   tACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa               <  2:3194976/150‑1 (MQ=255)
                                                                                                                                   tACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATGTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTGTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa               <  2:3602883/150‑1 (MQ=255)
                                                                                                                                     cGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACAGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCTGCAGATCAACAGAATCTTTATCCTGACGGAAGGCAtt             >  1:4246929/1‑150 (MQ=255)
                                                                                                                                        cGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAg          <  2:3131186/150‑1 (MQ=255)
                                                                                                                                         gggTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgt         >  2:1799150/1‑150 (MQ=255)
                                                                                                                                          ggTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtt        <  1:3682507/150‑1 (MQ=255)
                                                                                                                                          ggTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtt        <  2:4999443/150‑1 (MQ=255)
                                                                                                                                           gTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgttt       <  2:342729/150‑1 (MQ=255)
                                                                                                                                             gAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtttgt     >  2:2442919/1‑150 (MQ=255)
                                                                                                                                             gAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtttgt     >  2:3090767/1‑150 (MQ=255)
                                                                                                                                              aaCTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtttgtt    >  1:1034439/1‑150 (MQ=255)
                                                                                                                                               tcTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtttgttt   <  1:642680/149‑1 (MQ=255)
                                                                                                                                                cTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGtttt  <  1:2387953/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTT  >  CP009273/2127981‑2128274

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: