Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 361,806 G→A 93.7% 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,8060GA93.7% 15035.9 / inf 4571intergenic (‑63/+14)lacI/mhpRDNA‑binding transcriptional repressor/mhp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (124/161);  new base A (2147/2130);  total (2274/2297)
Fisher's exact test for biased strand distribution p-value = 3.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

Only 100 of 4577 total aligned reads displayed.

GGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCT  >  CP009273/361659‑361954
                                                                                                                                                   |                                                                                                                                                    
 gctggcCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCAcc                                                                                                                                                    <  2:3382803/149‑1 (MQ=255)
  ctggcCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCa                                                                                                                                                   <  1:197665/149‑1 (MQ=255)
   tggcCTGGTTCACCACGAGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCAt                                                                                                                                                  <  2:142947/149‑1 (MQ=255)
   aggcCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCAt                                                                                                                                                  <  2:4648801/148‑1 (MQ=255)
    ggcCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTTAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCAtt                                                                                                                                                 <  2:2797441/149‑1 (MQ=255)
     gcCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTc                                                                                                                                                <  2:2191128/149‑1 (MQ=255)
     gcCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTc                                                                                                                                                <  2:434224/149‑1 (MQ=255)
        tGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGAt                                                                                                                                             <  1:2847736/149‑1 (MQ=255)
        tGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGAt                                                                                                                                             <  1:4521378/149‑1 (MQ=255)
        tGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGAt                                                                                                                                             >  1:3315622/1‑149 (MQ=255)
        tGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCGTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGAt                                                                                                                                             <  1:4058239/149‑1 (MQ=255)
          gTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATgg                                                                                                                                           <  2:2639142/149‑1 (MQ=255)
              cccacGCGGGAAACGGTCAGAGAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTAACATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTc                                                                                                                                       <  2:2939333/148‑1 (MQ=255)
              accacGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTc                                                                                                                                       <  2:898342/149‑1 (MQ=255)
               ccacGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCa                                                                                                                                      >  2:1162751/1‑149 (MQ=255)
                  cgcgGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAtt                                                                                                                                   >  1:2753482‑M2/1‑147 (MQ=255)
                    cgGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaa                                                                                                                                 <  1:46739‑M2/149‑5 (MQ=255)
                     gggAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTTAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaat                                                                                                                                >  2:4737584‑M2/1‑144 (MQ=255)
                     gggAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaat                                                                                                                                >  2:1595013‑M2/1‑144 (MQ=255)
                       gAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTTAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaatta                                                                                                                              <  2:1507734‑M2/149‑8 (MQ=255)
                        aaaCGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaa                                                                                                                             >  1:2150373‑M2/1‑141 (MQ=255)
                          aCGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAatgca                                                                                                                           <  1:4742768/149‑1 (MQ=255)
                            ggTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggc                                                                                                                         <  1:972686‑M2/149‑13 (MQ=255)
                             gTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggca                                                                                                                        <  1:706209‑M2/149‑14 (MQ=255)
                             gTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTTAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggca                                                                                                                        >  2:3053815‑M2/1‑136 (MQ=255)
                             gTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggca                                                                                                                        <  1:2356759‑M2/149‑14 (MQ=255)
                             gTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggca                                                                                                                        >  1:1154612‑M2/1‑136 (MQ=255)
                              tCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcat                                                                                                                       >  2:36604‑M2/1‑135 (MQ=255)
                              tCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGcc                                                                                                                       >  1:613666/1‑149 (MQ=255)
                                  aTAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaa                                                                                                                   <  2:1245866‑M2/149‑19 (MQ=255)
                                   tAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaat                                                                                                                  >  1:3388158‑M2/1‑130 (MQ=255)
                                    aaGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaata                                                                                                                 >  2:2279174‑M2/1‑129 (MQ=255)
                                      gagaCACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaa                                                                                                               >  1:3782756‑M2/1‑127 (MQ=255)
                                       agaCACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTTAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaa                                                                                                              <  1:395071‑M2/149‑24 (MQ=255)
                                       agaCACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaa                                                                                                              >  1:3635923‑M2/1‑126 (MQ=255)
                                        gaCACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaac                                                                                                             >  1:3133342‑M2/1‑125 (MQ=255)
                                           acCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaa                                                                                                          >  2:2712879‑M2/1‑122 (MQ=255)
                                             cGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaag                                                                                                        >  2:785434‑M2/1‑120 (MQ=255)
                                             cGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaag                                                                                                        >  1:1000044‑M2/1‑120 (MQ=255)
                                               gCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTTAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggc                                                                                                      <  2:3563847‑M2/149‑32 (MQ=255)
                                                 aTACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctc                                                                                                    >  1:4299524‑M2/1‑116 (MQ=255)
                                                  tACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATGCCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaggggtca                                                                                                   <  1:2448542‑M2/149‑35 (MQ=255)
                                                  tACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctca                                                                                                   <  1:4457800‑M2/149‑35 (MQ=255)
                                                  tACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctca                                                                                                   <  2:529280‑M2/149‑35 (MQ=255)
                                                    ctctGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCGTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagt                                                                                                 >  1:140508‑M2/1‑113 (MQ=255)
                                                      ctGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcg                                                                                               >  1:2919147‑M2/1‑111 (MQ=255)
                                                            cATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaatcgaaaggctcagtcgaaagac                                                                                         <  2:3892643‑M2/149‑45 (MQ=255)
                                                             aTCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagact                                                                                        <  1:3704800‑M2/149‑46 (MQ=255)
                                                               cGTATAACGTTACTGGTTTCACATTCACCACCCTTAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgg                                                                                      >  2:629999‑M2/1‑102 (MQ=255)
                                                                    aaCGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctt                                                                                 >  1:2969650‑M2/1‑97 (MQ=255)
                                                                     aCGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggccttt                                                                                <  1:2250453‑M2/149‑54 (MQ=255)
                                                                          aCTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttt                                                                           >  2:3754910‑M2/1‑91 (MQ=255)
                                                                           cTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattacaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgtttt                                                                          >  2:3475908‑M2/1‑90 (MQ=255)
                                                                            tGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgtttta                                                                         <  2:4160717‑M2/149‑61 (MQ=255)
                                                                             ggTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCcttctt                                                                        <  1:187657/149‑1 (MQ=255)
                                                                             ggTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttat                                                                        >  2:4403671‑M2/1‑88 (MQ=255)
                                                                              gTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatc                                                                       >  2:2103141‑M2/1‑87 (MQ=255)
                                                                               tttCACATTCACCACCCTTAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctgtcgttgtatct                                                                      <  2:3287810‑M2/149‑64 (MQ=255)
                                                                                   acaTTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctctgtcgaaagactgggcctttcgttttatctgttg                                                                  <  1:1240150‑M2/149‑68 (MQ=255)
                                                                                   acaTTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttg                                                                  >  2:1840309‑M2/1‑82 (MQ=255)
                                                                                     aTTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGt                                                                <  1:1074252/149‑1 (MQ=255)
                                                                                     aTTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtt                                                                >  1:3564774‑M2/1‑80 (MQ=255)
                                                                                      ttCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgttt                                                               >  2:2014421‑M2/1‑79 (MQ=255)
                                                                                      gtcACCACCCAGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgttt                                                               >  2:2346057‑M2/2‑79 (MQ=255)
                                                                                       tCACCCCCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCAGACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttg                                                              >  2:3218010‑M2/1‑78 (MQ=255)
                                                                                       tCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttg                                                              >  1:3250340‑M2/1‑78 (MQ=255)
                                                                                       tCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttg                                                              >  2:3139092‑M2/1‑78 (MQ=255)
                                                                                        caccacCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATGCGATGGTGTCAAttaattaaaggcatcaaataaaacgacaggctcagtcgaaagactgggcctttcgttttatctgttgtttgt                                                             >  2:1945081‑M2/1‑77 (MQ=255)
                                                                                          ccaccCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGCTGGTGTCAAttaattaaaggcatcaaataactcgaaaggctcagtcgaaagactgggcctttcgttttatctgttggttggcg                                                           <  2:2887674‑M2/149‑75 (MQ=255)
                                                                                             cccTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtg                                                        >  1:1306211‑M2/1‑72 (MQ=255)
                                                                                              ccTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtga                                                       <  2:1708547‑M2/149‑79 (MQ=255)
                                                                                               cTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTAccc                                                      <  1:1712490/149‑1 (MQ=255)
                                                                                               cTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaa                                                      >  2:3822950‑M2/1‑70 (MQ=255)
                                                                                                 gAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATGCGAGGGTGGCAAttacttaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggccgtgcggtttatctgttgtttgtcggtgaacg                                                    <  2:2472579‑M2/149‑82 (MQ=255)
                                                                                                  aaTTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattcaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgc                                                   <  2:2565004‑M2/149‑83 (MQ=255)
                                                                                                  aaTTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgc                                                   <  2:4565688‑M2/149‑83 (MQ=255)
                                                                                                   aTTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgct                                                  <  1:1571868‑M2/149‑84 (MQ=255)
                                                                                                       aCTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCTCCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcc                                              <  2:3963440‑M2/149‑88 (MQ=255)
                                                                                                       aCTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcc                                              >  1:3473859‑M2/1‑62 (MQ=255)
                                                                                                         tctctTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctg                                            >  1:789494‑M2/1‑60 (MQ=255)
                                                                                                          ctcgTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaacggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctga                                           >  1:4219‑M2/1‑59 (MQ=255)
                                                                                                             ttCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagta                                        <  2:1331218‑M2/149‑94 (MQ=255)
                                                                                                             ttCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagta                                        <  2:1064075‑M2/149‑94 (MQ=255)
                                                                                                                cGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtagga                                     >  1:884848‑M2/1‑53 (MQ=255)
                                                                                                                cGGGCGCGATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtagga                                     >  1:2547124‑M2/1‑53 (MQ=255)
                                                                                                                 gggCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggac                                    >  1:2675122‑M2/1‑52 (MQ=255)
                                                                                                                 gggCGCTAACATGCCATACCGCGAAAGGTTATGCACCCAACGATGGTGTCCAttaaaaaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttgtatctgtagtttgtcggttaacgctctcctcagtcggac                                    <  1:2063599‑M2/149‑98 (MQ=255)
                                                                                                                      cTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatc                               >  1:3061437‑M2/1‑47 (MQ=255)
                                                                                                                      cTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatc                               <  1:1478512‑M2/149‑103 (MQ=255)
                                                                                                                       tATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatcc                              <  1:1640797‑M2/149‑104 (MQ=255)
                                                                                                                          cATGCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctggtgtttgtcggtgaacgctctcctgagtaggacaaatccgcc                           <  2:4242657‑M2/149‑107 (MQ=255)
                                                                                                                             gCCATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatccgccgcc                        <  1:1805908‑M2/149‑110 (MQ=255)
                                                                                                                              ccATACCTCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggaatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatccgccgccc                       >  2:699941‑M2/1‑39 (MQ=255)
                                                                                                                              ccATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatccgccgccc                       >  2:405769‑M2/1‑39 (MQ=255)
                                                                                                                               cATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatccgccgccct                      <  1:52421‑M2/149‑112 (MQ=255)
                                                                                                                               cATACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagaatgggcctttcgttttatctgttgtttttctgtgaacgctctactgagtaggcaaaatacgccgccct                      <  1:43866‑M2/149‑112 (MQ=255)
                                                                                                                                 tACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatccgccgccctag                    <  1:4593913‑M2/149‑114 (MQ=255)
                                                                                                                                 tACCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctttcgttttatctgttgtttgtcggtgaacgctctcctgagtaggacaaatccgccgccctag                    <  1:4144453‑M2/149‑114 (MQ=255)
                                                                                                                                  aCCGCGAAAGGTTTTGCACCATTCGATGGTGTCAAttaattaaaggcatcaaataaaacgaaaggctcagtcgaaagactgggcctgtcgttttatctggtgtttgtcggtgaacgctctcctgagtaggacaaatccgccgccctaga                   <  2:4476064‑M2/149‑115 (MQ=255)
                                                                                                                                             ttttGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCAt        >  1:3870176/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCAACGTAAATGCATGCCGCTTCGCCTTCCGGCCACCAGAATAGCCTGCGATTCAACCCCTTCTTCGATCTGTTTTGCTACCCGTTGTAGCGCCGGAAGATGCTTTTCCGCTGCCTGTTCAATGGTCATTGCGCT  >  CP009273/361659‑361954

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