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.7% 987.3 / 1.1 299K336* (AAA→TAA) wzccolanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base A (130/164);  total (132/167)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
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 299 total aligned reads displayed.

CTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTT  >  CP009273/2127982‑2128274
                                                                                                                                                  |                                                                                                                                                  
ctgctgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGgt                                                                                                                                                 >  2:3284866/1‑150 (MQ=255)
ctgctgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGgt                                                                                                                                                 <  2:4054702/150‑1 (MQ=255)
 tgctgACCAGACTCGACCTCGCGGGTCAGCCGGACAATCTCCAGCTGGGTGTGCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTGTGCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGGACGCCGGGTGAACTTAGgtg                                                                                                                                                <  2:1521017/150‑1 (MQ=255)
 tgctgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGgtg                                                                                                                                                >  1:3474676/1‑150 (MQ=255)
     gACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACa                                                                                                                                            >  1:1653242/1‑150 (MQ=255)
     gACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACa                                                                                                                                            <  1:1576690/150‑1 (MQ=255)
          gACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTg                                                                                                                                       <  1:3159511/150‑1 (MQ=255)
           aCTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTgg                                                                                                                                      <  1:3517666/150‑1 (MQ=255)
            cTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGga                                                                                                                                     <  2:1344917/150‑1 (MQ=255)
            cTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGga                                                                                                                                     >  2:2807813/1‑150 (MQ=255)
               gACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAt                                                                                                                                  >  2:2436492/1‑150 (MQ=255)
               gACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAt                                                                                                                                  >  2:2276837/1‑150 (MQ=255)
                aCATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAtt                                                                                                                                 <  1:1103529/150‑1 (MQ=255)
                 cATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAttt                                                                                                                                >  1:3841311/1‑150 (MQ=255)
                    cgcgGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCg                                                                                                                             <  1:262787/150‑1 (MQ=255)
                       gggTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCTACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCct                                                                                                                          >  1:7345/1‑150 (MQ=255)
                       gggTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCct                                                                                                                          <  1:4452968/150‑1 (MQ=255)
                       gggTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCct                                                                                                                          >  2:3714431/1‑150 (MQ=255)
                       gggTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCct                                                                                                                          >  1:3623820/1‑150 (MQ=255)
                       gggTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCct                                                                                                                          >  2:346317/1‑150 (MQ=255)
                          tCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCtt                                                                                                                       >  1:2254802/1‑150 (MQ=255)
                          tCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCtt                                                                                                                       <  1:2046874/150‑1 (MQ=255)
                           cAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCttt                                                                                                                      <  1:4469040/150‑1 (MQ=255)
                            aGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTa                                                                                                                     >  2:2703602/1‑150 (MQ=255)
                            aGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTa                                                                                                                     <  2:174400/150‑1 (MQ=255)
                                ggACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAgg                                                                                                                 >  1:43698/1‑150 (MQ=255)
                                ggACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAgg                                                                                                                 >  1:943392/1‑150 (MQ=255)
                                 gACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGt                                                                                                                >  2:475959/1‑150 (MQ=255)
                                  aCAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTc                                                                                                               >  2:2868001/1‑150 (MQ=255)
                                   cAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCa                                                                                                              <  1:1581761/150‑1 (MQ=255)
                                   cAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCa                                                                                                              <  1:1228456/150‑1 (MQ=255)
                                      tctcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgtt                                                                                                           >  1:4192586/1‑150 (MQ=255)
                                      tctcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgtt                                                                                                           <  2:1419966/150‑1 (MQ=255)
                                       ctcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttc                                                                                                          <  1:3003914/150‑1 (MQ=255)
                                          ctgctgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttcgtc                                                                                                       <  1:3740979/150‑2 (MQ=255)
                                            gctgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCa                                                                                                     <  2:4218728/150‑1 (MQ=255)
                                             ctgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCaa                                                                                                    >  1:2359840/1‑150 (MQ=255)
                                               gggTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACt                                                                                                  >  2:3214171/1‑150 (MQ=255)
                                                ggTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTg                                                                                                 <  2:1114592/150‑1 (MQ=255)
                                                 gTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTgc                                                                                                <  2:1806457/150‑1 (MQ=255)
                                                     tCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCAt                                                                                            >  1:3662075/1‑150 (MQ=255)
                                                      cGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATc                                                                                           <  1:3406480/150‑1 (MQ=255)
                                                       ggCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCg                                                                                          >  2:4276108/1‑150 (MQ=255)
                                                          aTCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATg                                                                                       >  2:802182/1‑150 (MQ=255)
                                                                 ttACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCa                                                                                <  1:665185/150‑1 (MQ=255)
                                                                  tACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTTTACAGCTTGGAGATTTCAGCCTCTTTAAAGGTCAGTTAGTTCAACTGCGCATCGATGTTCACCat                                                                               >  1:921240/1‑150 (MQ=255)
                                                                   aCGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatc                                                                              <  2:662917/150‑1 (MQ=255)
                                                                    cgcgACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcg                                                                             <  1:4030457/150‑1 (MQ=255)
                                                                    cgcgACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTAGTTCAACTGCGCATCGATGTTCACCatcg                                                                             >  1:1696068/1‑150 (MQ=255)
                                                                     gcgACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcga                                                                            <  1:5159165/150‑1 (MQ=255)
                                                                       gACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcgaat                                                                          >  2:1204450/1‑150 (MQ=255)
                                                                         ccGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcgaatcg                                                                        <  1:4245007/150‑1 (MQ=255)
                                                                            ttAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGc                                                                     >  1:823830/1‑150 (MQ=255)
                                                                            ttAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGc                                                                     <  1:3080236/150‑1 (MQ=255)
                                                                              aaGGTTGGCTTGTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCAc                                                                   <  2:3897049/150‑1 (MQ=255)
                                                                               aGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCAcc                                                                  >  2:1640621/1‑150 (MQ=255)
                                                                               aGTTTGGCTTTTTAGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCAcc                                                                  >  2:3120887/1‑150 (MQ=255)
                                                                                gTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCg                                                                 >  1:1420891/1‑150 (MQ=255)
                                                                                 tttGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGc                                                                >  1:4733587/1‑150 (MQ=255)
                                                                                  ttGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCt                                                               >  2:1929053/1‑150 (MQ=255)
                                                                                    ggCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCttt                                                             <  2:3645019/150‑1 (MQ=255)
                                                                                     gCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCtttt                                                            <  2:4073059/150‑1 (MQ=255)
                                                                                     gCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCtttt                                                            <  2:4487275/150‑1 (MQ=255)
                                                                                      cTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTg                                                           >  2:2920995/1‑150 (MQ=255)
                                                                                        ttttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTCGGTGTACAGCTTGGAGATTTCCGCCGCTTTAAAGGTCAGTTCGTTCAACTGCGCCTCGATGTTCACCATCGAATCGAGCACCGCTTTTGCt                                                         >  2:541372/1‑150 (MQ=255)
                                                                                        ttttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGGTCACCATCGAATCGAGCACCGCTTTTGCt                                                         <  2:4962234/150‑1 (MQ=255)
                                                                                         tttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCtt                                                        >  2:897450/1‑150 (MQ=255)
                                                                                           tCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTcc                                                      <  2:1270470/150‑1 (MQ=255)
                                                                                              tCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGc                                                   <  2:2525511/150‑1 (MQ=255)
                                                                                              tCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGc                                                   <  1:1845515/150‑1 (MQ=255)
                                                                                               cTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCg                                                  >  2:4874260/1‑150 (MQ=255)
                                                                                                 tCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGc                                                >  1:2510718/1‑150 (MQ=255)
                                                                                                    aGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTACAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGa                                             <  2:2567904/150‑1 (MQ=255)
                                                                                                      cgcCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATc                                           <  2:4346853/150‑1 (MQ=255)
                                                                                                        ccTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCaa                                         <  2:3582408/150‑1 (MQ=255)
                                                                                                         cTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAAc                                        <  2:5069813/150‑1 (MQ=255)
                                                                                                          tGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACa                                       <  1:2798079/150‑1 (MQ=255)
                                                                                                                   tcCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTa                              <  1:1000786/150‑1 (MQ=255)
                                                                                                                   tcCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTa                              <  1:1514024/150‑1 (MQ=255)
                                                                                                                         agcGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGa                        >  1:3266545/1‑150 (MQ=255)
                                                                                                                              gCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGaa                   <  2:83591/150‑1 (MQ=255)
                                                                                                                              gCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGaa                   <  1:4614047/150‑1 (MQ=255)
                                                                                                                              gCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGaa                   <  1:2658236/150‑1 (MQ=255)
                                                                                                                               cGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAg                  <  2:3759301/150‑1 (MQ=255)
                                                                                                                               cGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAg                  <  1:3232505/150‑1 (MQ=255)
                                                                                                                                  tACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa               >  1:4940007/1‑150 (MQ=255)
                                                                                                                                  tACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa               <  1:5040350/150‑1 (MQ=255)
                                                                                                                                   aCGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCAt              <  2:3382620/150‑1 (MQ=255)
                                                                                                                                    cGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCAtt             <  1:3923961/150‑1 (MQ=255)
                                                                                                                                     gCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTc            <  1:883300/150‑1 (MQ=255)
                                                                                                                                        gggTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgt         <  1:2781078/150‑1 (MQ=255)
                                                                                                                                         ggTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtt        <  2:44323/150‑1 (MQ=255)
                                                                                                                                         ggTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtt        <  1:5114139/150‑1 (MQ=255)
                                                                                                                                          gTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgttt       <  2:3876675/150‑1 (MQ=255)
                                                                                                                                          gTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgttt       <  1:2115/150‑1 (MQ=255)
                                                                                                                                          gTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgttt       <  2:392823/150‑1 (MQ=255)
                                                                                                                                          gTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATAGAGCACCGCTTTTGATTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgttt       >  1:4118184/1‑150 (MQ=255)
                                                                                                                                          gTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGGTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgttt       <  2:1836795/150‑1 (MQ=255)
                                                                                                                                               cTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTATAGGTCAGTTCGTTCAACTGCGCATCGATGGTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGtttt  <  2:5138982/150‑1 (MQ=255)
                                                                                                                                               cTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGtttt  <  1:101451/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTT  >  CP009273/2127982‑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: