Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,680 C→A 16.2% R355L (CGA→CTA)  lacI ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273360,6800CA16.2% 36439.4 / inf 11219R355L (CGA→CTA) lacIDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (4724/4665);  new base A (908/920);  total (5633/5586)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 11363 total aligned reads displayed.

GGGTCGGGTGAATCTTGTTGATGGCGATATGGATGTGCAGGTTGTCGGTGTCGTGATGCACGGCACTGACGCGCTGATGCTCGGCGAAGCCAAGCCCAGCGCAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGC  >  CP009273/360531‑360829
                                                                                                                                                     |                                                                                                                                                     
 tacgacctgcacagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCg                                                                                                                                                      >  2:1028526‑M1/121‑150 (MQ=255)
   cgacctgcacagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCggg                                                                                                                                                    <  2:3243151‑M1/32‑1 (MQ=255)
    gacctgcacagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGa                                                                                                                                                   <  1:3012817‑M1/33‑1 (MQ=255)
       ctgcacagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAAc                                                                                                                                                <  1:2337497‑M1/36‑1 (MQ=255)
       ctgcacagccataccacagcttccgatggctgcctgacgccagaagcataggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAAc                                                                                                                                                <  1:443795‑M1/36‑1 (MQ=255)
         gcacagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCt                                                                                                                                              >  2:2501970‑M1/113‑150 (MQ=255)
           acagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGt                                                                                                                                            >  2:1273686‑M1/111‑150 (MQ=255)
           acagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGt                                                                                                                                            >  1:2163945‑M1/111‑150 (MQ=255)
            cagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTc                                                                                                                                           >  2:840534‑M1/110‑150 (MQ=255)
            cagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTc                                                                                                                                           <  2:2416834‑M1/41‑1 (MQ=255)
              gccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGt                                                                                                                                         >  2:4201308‑M1/108‑150 (MQ=255)
              gccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGt                                                                                                                                         <  1:1849278‑M1/43‑1 (MQ=255)
                  taccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCa                                                                                                                                     >  1:2664259‑M1/104‑150 (MQ=255)
                     cacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCt                                                                                                                                  <  2:2327014‑M1/50‑1 (MQ=255)
                        agcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCa                                                                                                                               >  2:3010593‑M1/98‑150 (MQ=255)
                           ttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTa                                                                                                                            <  2:3851028‑M1/56‑1 (MQ=255)
                            tccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTaa                                                                                                                           <  1:1438275‑M1/57‑1 (MQ=255)
                               gatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGa                                                                                                                        <  2:3076493‑M1/60‑1 (MQ=255)
                               gatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGa                                                                                                                        >  1:3599226‑M1/91‑150 (MQ=255)
                                 tggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAAt                                                                                                                      <  2:2227382‑M1/62‑1 (MQ=255)
                                   gctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCg                                                                                                                    >  2:1096613‑M1/87‑150 (MQ=255)
                                     tgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGc                                                                                                                  >  1:1511976‑M1/85‑150 (MQ=255)
                                      gcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGcc                                                                                                                 >  2:1973971‑M1/84‑150 (MQ=255)
                                        ctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCaa                                                                                                               <  1:1162890‑M1/69‑1 (MQ=255)
                                         tgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAAc                                                                                                              <  2:4376990‑M1/70‑1 (MQ=255)
                                         tgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacataaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAAc                                                                                                              <  1:4508360‑M1/70‑1 (MQ=255)
                                         tgacgccaaaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAAc                                                                                                              <  1:2840736‑M1/70‑1 (MQ=255)
                                               cagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCgg                                                                                                        <  2:1336997‑M1/76‑1 (MQ=255)
                                                  aagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGga                                                                                                     >  2:1662420‑M1/72‑150 (MQ=255)
                                                     cattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGagag                                                                                                  >  2:1184920‑M1/69‑150 (MQ=255)
                                                      attggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAgg                                                                                                 <  2:597538‑M1/83‑1 (MQ=255)
                                                          gtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGt                                                                                             >  2:3776584‑M1/64‑150 (MQ=255)
                                                          cACGGCACTGACGCGCTGATGCTCGGCGAAGCCAAGCCCAGCGCAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGt                                                                                             >  2:4455493/1‑150 (MQ=255)
                                                             caccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTg                                                                                          <  1:2588983‑M1/90‑1 (MQ=255)
                                                                  tgcagtcgatgctaagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTAt                                                                                     >  2:1577234‑M1/56‑150 (MQ=255)
                                                                     agtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTgg                                                                                  >  2:760631‑M1/53‑150 (MQ=255)
                                                                     agtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTgg                                                                                  <  1:2005522‑M1/98‑1 (MQ=255)
                                                                      gtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTggg                                                                                 >  1:3522089‑M1/52‑150 (MQ=255)
                                                                       tcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGATAGTCGGTTTTCGTATTGGGc                                                                                >  2:3686099‑M1/51‑150 (MQ=255)
                                                                       tcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGc                                                                                >  1:3345005‑M1/51‑150 (MQ=255)
                                                                       tcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGc                                                                                >  1:2418121‑M1/51‑150 (MQ=255)
                                                                        cgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGgcg                                                                               <  2:1807347‑M1/101‑2 (MQ=255)
                                                                          atgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTgggcgcc                                                                             <  2:3510995‑M1/103‑4 (MQ=255)
                                                                             ataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATAGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAgg                                                                          >  1:1079119‑M1/45‑150 (MQ=255)
                                                                              taagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAggg                                                                         <  1:2750395‑M1/107‑1 (MQ=255)
                                                                               aagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggt                                                                        >  1:3091978‑M1/43‑150 (MQ=255)
                                                                               aagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggt                                                                        >  1:4350921‑M1/43‑150 (MQ=255)
                                                                                agctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtg                                                                       <  2:2766736‑M1/109‑1 (MQ=255)
                                                                                   tgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtt                                                                    >  1:575735‑M1/39‑150 (MQ=255)
                                                                                    gtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGttt                                                                   >  2:258632‑M1/38‑150 (MQ=255)
                                                                                     tcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtttt                                                                  >  1:177802‑M1/37‑150 (MQ=255)
                                                                                     tcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtttt                                                                  >  1:840359‑M1/37‑150 (MQ=255)
                                                                                      caaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGttttt                                                                 <  1:1346541‑M1/115‑1 (MQ=255)
                                                                                      caaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGttttt                                                                 <  2:2152365‑M1/115‑1 (MQ=255)
                                                                                       aaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttc                                                                <  1:4271830‑M1/116‑1 (MQ=255)
                                                                                       aaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttc                                                                >  1:1924404‑M1/35‑150 (MQ=255)
                                                                                        aaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttct                                                               <  2:461181‑M1/117‑1 (MQ=255)
                                                                                         accagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctt                                                              >  2:2679675‑M1/33‑150 (MQ=255)
                                                                                          ccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttt                                                             <  1:393930‑M1/119‑1 (MQ=255)
                                                                                           cagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctttt                                                            >  2:4279686‑M1/31‑150 (MQ=255)
                                                                                           cagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctttt                                                            >  2:4541479‑M1/31‑150 (MQ=255)
                                                                                            agatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttttc                                                           <  2:341261‑M1/121‑1 (MQ=255)
                                                                                            agatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttttc                                                           >  2:3145974‑M1/30‑150 (MQ=255)
                                                                                            agatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttttc                                                           >  1:208512‑M1/30‑150 (MQ=255)
                                                                                             gatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCa                                                          <  2:142126‑M1/122‑1 (MQ=255)
                                                                                               tcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAcc                                                        <  2:1904526‑M1/124‑1 (MQ=255)
                                                                                                caattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                       >  1:2506863‑M1/26‑150 (MQ=255)
                                                                                                  attcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                     >  2:910000‑M1/24‑150 (MQ=255)
                                                                                                  attcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                     <  1:3257230‑M1/127‑1 (MQ=255)
                                                                                                  attcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                     <  1:3774931‑M1/127‑1 (MQ=255)
                                                                                                    tcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                   >  1:4122224‑M1/22‑150 (MQ=255)
                                                                                                    tcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                   >  1:1708044‑M1/22‑150 (MQ=255)
                                                                                                       cgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                >  1:4196041‑M1/19‑150 (MQ=255)
                                                                                                        gctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACg                                               <  1:1608709‑M1/133‑1 (MQ=255)
                                                                                                         ctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTAGTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                              <  2:519491‑M1/134‑1 (MQ=255)
                                                                                                          taacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACggg                                             >  2:402684‑M1/16‑150 (MQ=255)
                                                                                                            acttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTTATGAATAGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCa                                           <  1:1000315‑M1/137‑1 (MQ=255)
                                                                                                              ttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAc                                         >  1:1250103‑M1/12‑150 (MQ=255)
                                                                                                                acattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAg                                       >  1:369074‑M1/10‑150 (MQ=255)
                                                                                                                acattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAg                                       >  1:31745‑M1/10‑150 (MQ=255)
                                                                                                                   ttaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTg                                    <  1:482995‑M1/144‑1 (MQ=255)
                                                                                                                     aattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAt                                  >  2:1498264‑M1/5‑150 (MQ=255)
                                                                                                                      attGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt                                 >  2:1734176‑M1/4‑150 (MQ=255)
                                                                                                                      attGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt                                 <  1:3845108‑M1/147‑1 (MQ=255)
                                                                                                                           gTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGcgga                            <  1:667471/150‑4 (MQ=255)
                                                                                                                             tGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTc                          <  1:2925436/150‑1 (MQ=255)
                                                                                                                              gcgcTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCa                         <  2:4119524/150‑1 (MQ=255)
                                                                                                                                  tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGc                     >  1:941135/1‑150 (MQ=255)
                                                                                                                                   cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGcc                    >  1:3442213/1‑150 (MQ=255)
                                                                                                                                    aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACTGGCAACAGCTGATTGCCCTTCACCGCCt                   <  1:2256953/150‑1 (MQ=255)
                                                                                                                                    aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                   <  1:459399/150‑1 (MQ=255)
                                                                                                                                    aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                   <  2:3601320/150‑1 (MQ=255)
                                                                                                                                     cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTg                  <  1:750975/150‑1 (MQ=255)
                                                                                                                                     cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTg                  <  2:2940781/150‑1 (MQ=255)
                                                                                                                                       gCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGc                >  2:3332579/1‑150 (MQ=255)
                                                                                                                                         ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCAATAATGAATAGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGccc              >  2:2854320/1‑150 (MQ=255)
                                                                                                                                          cGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCt             <  2:2072136/150‑1 (MQ=255)
                                                                                                                                              ttCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgaga         <  2:685663/150‑1 (MQ=255)
                                                                                                                                                  aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGt     >  2:3931925/1‑150 (MQ=255)
                                                                                                                                                    tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTg   <  1:4031340/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGTCGGGTGAATCTTGTTGATGGCGATATGGATGTGCAGGTTGTCGGTGTCGTGATGCACGGCACTGACGCGCTGATGCTCGGCGAAGCCAAGCCCAGCGCAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGC  >  CP009273/360531‑360829

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: