Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 361,806 G→A 100% intergenic (‑63/+14) lacI ← / ← mhpR DNA‑binding transcriptional repressor/mhp operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273361,8060GA96.6% 32908.4 / inf 9010intergenic (‑63/+14)lacI/mhpRDNA‑binding transcriptional repressor/mhp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (150/158);  new base A (4324/4372);  total (4476/4534)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 9010 total aligned reads displayed.

GCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGC  >  CP009273/361660‑361951
                                                                                                                                                  |                                                                                                                                                 
gctggcCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCAcc                                                                                                                                                 <  1:3327585/149‑1 (MQ=255)
  tggcCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCAt                                                                                                                                               <  2:1321389/149‑1 (MQ=255)
   ggcCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCAtt                                                                                                                                              >  1:1074025/1‑149 (MQ=255)
   ggcCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCAtt                                                                                                                                              <  1:4520510/149‑1 (MQ=255)
   ggcCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCAtt                                                                                                                                              <  1:1001969/149‑1 (MQ=255)
    gcCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTc                                                                                                                                             >  1:3424115/1‑149 (MQ=255)
    gcCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTc                                                                                                                                             >  2:124474/1‑149 (MQ=255)
       tGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGAt                                                                                                                                          <  1:1495417/149‑1 (MQ=255)
       tGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGAt                                                                                                                                          >  2:2084354/1‑149 (MQ=255)
         gTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTGTGCACCAGTCGAggg                                                                                                                                        >  1:475166/1‑149 (MQ=255)
            caccacGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGgtgt                                                                                                                                     <  2:3543078/149‑1 (MQ=255)
               cacGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCaa                                                                                                                                  <  1:3163262/149‑1 (MQ=255)
                aCGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAt                                                                                                                                 >  1:1879270‑M2/1‑148 (MQ=255)
                 cgcgGGAAACGGTCAGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAtt                                                                                                                                >  2:2606560‑M2/1‑147 (MQ=255)
                   cgGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaa                                                                                                                              >  1:134472‑M2/1‑145 (MQ=255)
                   cgGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaa                                                                                                                              >  2:2424981‑M2/1‑145 (MQ=255)
                   cgGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTATCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaa                                                                                                                              <  2:2011637‑M2/149‑5 (MQ=255)
                     ggAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaatt                                                                                                                            >  1:3506396‑M2/1‑143 (MQ=255)
                     ggAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaatt                                                                                                                            >  1:443544‑M2/1‑143 (MQ=255)
                        aaCGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAatgc                                                                                                                         >  1:3244920/1‑149 (MQ=255)
                        aaCGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTAGATGGTGTCAAttaattaaa                                                                                                                         <  1:4086638‑M2/149‑10 (MQ=255)
                         aCGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaag                                                                                                                        <  1:2914735‑M2/149‑11 (MQ=255)
                          cGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaagg                                                                                                                       <  1:4240898‑M2/149‑12 (MQ=255)
                             tCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcat                                                                                                                    >  1:4013659‑M2/1‑135 (MQ=255)
                             tCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcat                                                                                                                    <  2:1496470‑M2/149‑15 (MQ=255)
                               tGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatca                                                                                                                  <  2:531206‑M2/149‑17 (MQ=255)
                                gATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGCAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaa                                                                                                                 <  2:1679269‑M2/149‑18 (MQ=255)
                                 aTAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaa                                                                                                                >  1:2662980‑M2/1‑131 (MQ=255)
                                 aTAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaa                                                                                                                >  1:3851655‑M2/1‑131 (MQ=255)
                                     gagaCACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTCCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaa                                                                                                            >  1:2738391‑M2/1‑127 (MQ=255)
                                       gaCACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaac                                                                                                          >  2:2173377‑M2/1‑125 (MQ=255)
                                        acacCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacg                                                                                                         <  1:1266670‑M2/149‑26 (MQ=255)
                                            cGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaag                                                                                                     <  1:2842997‑M2/149‑30 (MQ=255)
                                            cGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaag                                                                                                     >  1:978356‑M2/1‑120 (MQ=255)
                                              gCATACTCTGCGACATAGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGAGGGTGTCAAttaattaaaggcatcaaataaaacgaaaggc                                                                                                   <  2:3707656‑M2/149‑32 (MQ=255)
                                               cATACTCTGCGACATCGTATCACGTTCCTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTGTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaccgaaaggct                                                                                                  >  2:3778764‑M2/1‑117 (MQ=255)
                                                aTACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctc                                                                                                 >  1:1720640‑M2/1‑116 (MQ=255)
                                                   ctctGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagt                                                                                              <  2:3953674‑M2/149‑37 (MQ=255)
                                                   ctctGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagt                                                                                              >  1:2322733‑M2/1‑113 (MQ=255)
                                                   ctctGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagt                                                                                              >  2:4381332‑M2/1‑113 (MQ=255)
                                                    tctGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtc                                                                                             >  1:1964729‑M2/1‑112 (MQ=255)
                                                          aCATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaaga                                                                                       <  1:1640676‑M2/149‑44 (MQ=255)
                                                            aTCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagact                                                                                     <  1:1802397‑M2/149‑46 (MQ=255)
                                                            aTCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagact                                                                                     >  1:575175‑M2/1‑104 (MQ=255)
                                                            aTCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagact                                                                                     <  1:751374‑M2/149‑46 (MQ=255)
                                                             tCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactg                                                                                    <  1:3759606‑M2/149‑47 (MQ=255)
                                                               gTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactggg                                                                                  <  1:4604530‑M2/149‑49 (MQ=255)
                                                               gTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactggg                                                                                  <  2:4213780‑M2/149‑49 (MQ=255)
                                                                  taACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcct                                                                               <  2:4658997‑M2/149‑52 (MQ=255)
                                                                  taACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcct                                                                               >  1:2566795‑M2/1‑98 (MQ=255)
                                                                   aaCGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctt                                                                              >  2:1592225‑M2/1‑97 (MQ=255)
                                                                   aaCGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctg                                                                              <  2:3857053‑M2/149‑53 (MQ=255)
                                                                     cGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttc                                                                            >  1:1163024‑M2/1‑95 (MQ=255)
                                                                      gTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcg                                                                           <  2:4113961‑M2/149‑56 (MQ=255)
                                                                       ttACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgt                                                                          <  2:2246836‑M2/149‑57 (MQ=255)
                                                                        tACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgtt                                                                         <  2:3066493‑M2/149‑58 (MQ=255)
                                                                         aCTGGTTACACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttt                                                                        >  2:836816‑M2/1‑91 (MQ=255)
                                                                          cTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgtttt                                                                       <  2:726996‑M2/149‑60 (MQ=255)
                                                                          cTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgtttt                                                                       <  2:139921‑M2/149‑60 (MQ=255)
                                                                          cTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGAGGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgtttt                                                                       <  2:1936058‑M2/149‑60 (MQ=255)
                                                                           tGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgtttta                                                                      <  1:3593933‑M2/149‑61 (MQ=255)
                                                                           tGGTTTCACATACACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCTATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgtttta                                                                      >  1:4321623‑M2/1‑89 (MQ=255)
                                                                            ggTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgacaggctcagtcgaaagactgggcctttcgttttat                                                                     >  2:298280‑M2/1‑88 (MQ=255)
                                                                             gTTTCACATTCACCACCCTGACTTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatc                                                                    >  2:1834980‑M2/1‑87 (MQ=255)
                                                                             gTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgacaggctcagtcgaaagactgggcctttcgttttatc                                                                    >  1:1432641‑M2/1‑87 (MQ=255)
                                                                             gTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatc                                                                    <  1:300107‑M2/149‑63 (MQ=255)
                                                                              tttCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatct                                                                   <  1:2034135‑M2/149‑64 (MQ=255)
                                                                                      tCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttg                                                           >  2:2518887‑M2/1‑78 (MQ=255)
                                                                                        accaccCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaacggctcagtcgaaagactgggcctttcgttttatctgttgtttgtc                                                         >  2:3465071‑M2/1‑76 (MQ=255)
                                                                                         ccaccCTGAATGGACTCGCGTCCGGGCGCTATCATGCCATACCGCGACAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcg                                                        >  2:3620786‑M2/1‑75 (MQ=255)
                                                                                           accCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggt                                                      >  1:2166318‑M2/1‑73 (MQ=255)
                                                                                             ccTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtga                                                    <  1:893644‑M2/149‑79 (MQ=255)
                                                                                              cTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaa                                                   <  2:4576328‑M2/149‑80 (MQ=255)
                                                                                                gAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGt                                                 <  1:4145736/149‑1 (MQ=255)
                                                                                                gAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacg                                                 <  2:278935‑M2/149‑82 (MQ=255)
                                                                                                gAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacg                                                 <  2:1752817‑M2/149‑82 (MQ=255)
                                                                                                   ttGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaacgactgggccttgcgttgtatctgtggtttgtcggtgaacgctc                                              <  2:3392026‑M2/149‑85 (MQ=255)
                                                                                                   ttGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctc                                              <  2:287427‑M2/149‑85 (MQ=255)
                                                                                                    tGACTCTCTTCTGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctct                                             >  2:4031211‑M2/1‑64 (MQ=255)
                                                                                                      aCTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcc                                           <  1:2400661‑M2/149‑88 (MQ=255)
                                                                                                       ctctctTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgtgttatctgttgtttgtcggtgaacgctctcct                                          <  2:1164418‑M2/149‑89 (MQ=255)
                                                                                                         ctctTCCGTGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggatcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctga                                        >  2:315481‑M2/1‑59 (MQ=255)
                                                                                                             tCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtag                                    <  1:2250225‑M2/149‑95 (MQ=255)
                                                                                                              ccGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatatgttgtttgtcggtgaaagctctcctgagtagg                                   <  1:2103818‑M2/149‑96 (MQ=255)
                                                                                                               cGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGa                                  <  2:3229459/149‑1 (MQ=255)
                                                                                                                gggCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggac                                 <  2:2694356‑M2/149‑98 (MQ=255)
                                                                                                                gggCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggac                                 <  1:2487338‑M2/149‑98 (MQ=255)
                                                                                                                  gcgcTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaa                               <  1:832423‑M2/149‑100 (MQ=255)
                                                                                                                   cgcTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaa                              >  2:916678‑M2/1‑49 (MQ=255)
                                                                                                                     cTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatc                            <  2:632736‑M2/149‑103 (MQ=255)
                                                                                                                     cTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggaatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatc                            >  1:3946306‑M2/1‑47 (MQ=255)
                                                                                                                       aTCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTc                          >  1:1566608/1‑149 (MQ=255)
                                                                                                                       aTCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatccg                          <  2:1066618‑M2/149‑105 (MQ=255)
                                                                                                                         cATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatccgcc                        <  1:3672462‑M2/149‑107 (MQ=255)
                                                                                                                          aTGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatccgccg                       <  2:3323736‑M2/149‑108 (MQ=255)
                                                                                                                            gCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatccgccgcc                     <  1:664113‑M2/149‑110 (MQ=255)
                                                                                                                               aTACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatccgccgcccta                  <  2:4305917‑M2/149‑113 (MQ=255)
                                                                                                                                tACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatccgccgccctag                 >  1:3086027‑M2/1‑36 (MQ=255)
                                                                                                                                 aCCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGt                >  2:4489188/1‑149 (MQ=255)
                                                                                                                                           gTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCa      <  2:2338851/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGC  >  CP009273/361660‑361951

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: