Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,680 C→A 6.6% 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,6800CA 6.6% 49272.9 / inf 13037R355L (CGA→CTA) lacIDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (6098/6063);  new base A (435/439);  total (6534/6503)
Fisher's exact test for biased strand distribution p-value = 8.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 13199 total aligned reads displayed.

GGGTCGGGTGAATCTTGTTGATGGCGATATGGATGTGCAGGTTGTCGGTGTCGTGATGCACGGCACTGACGCGCTGATGCTCGGCGAAGCCAAGCCCAGCGCAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGC  >  CP009273/360531‑360829
                                                                                                                                                     |                                                                                                                                                     
ttacgacctgcacagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTc                                                                                                                                                       <  2:2499765‑M1/29‑1 (MQ=255)
ttacgacctgcacagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTc                                                                                                                                                       >  1:1585120‑M1/122‑150 (MQ=255)
 tacgacctgcacagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCg                                                                                                                                                      <  1:133020‑M1/30‑1 (MQ=255)
   cgacctgcacagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCggg                                                                                                                                                    >  2:1712551‑M1/119‑150 (MQ=255)
    gacctgcacagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGa                                                                                                                                                   <  2:2886781‑M1/33‑1 (MQ=255)
      cctgcacagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGaaa                                                                                                                                                 >  2:1627650‑M1/116‑150 (MQ=255)
         gcacagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCt                                                                                                                                              >  1:2384067‑M1/113‑150 (MQ=255)
           acagccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGt                                                                                                                                            >  2:625343‑M1/111‑150 (MQ=255)
             agccataccacagcttccgatggctgcctgacgccaggagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCg                                                                                                                                          <  1:1236978‑M1/42‑1 (MQ=255)
             agccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCg                                                                                                                                          >  1:3580859‑M1/109‑150 (MQ=255)
             agccataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCg                                                                                                                                          <  2:1131075‑M1/42‑1 (MQ=255)
                 ataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgctgtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGcc                                                                                                                                      <  1:3276062‑M1/46‑1 (MQ=255)
                 ataccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGcc                                                                                                                                      <  1:751128‑M1/46‑1 (MQ=255)
                  taccacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCa                                                                                                                                     >  2:764580‑M1/104‑150 (MQ=255)
                   accacagcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAg                                                                                                                                    <  2:3405212‑M1/48‑1 (MQ=255)
                        agcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCa                                                                                                                               <  1:496298‑M1/53‑1 (MQ=255)
                         gcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCAt                                                                                                                              <  1:1901323‑M1/54‑1 (MQ=255)
                         gcttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCAt                                                                                                                              <  2:1555482‑M1/54‑1 (MQ=255)
                          cttccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCAtt                                                                                                                             <  1:4147987‑M1/55‑1 (MQ=255)
                            tccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTaa                                                                                                                           <  1:2298151‑M1/57‑1 (MQ=255)
                            tccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTaa                                                                                                                           >  2:2432276‑M1/94‑150 (MQ=255)
                             ccgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAAt                                                                                                                          >  2:4146371‑M1/93‑150 (MQ=255)
                              cgatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATg                                                                                                                         >  2:3742198‑M1/92‑150 (MQ=255)
                               gatggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGa                                                                                                                        <  1:1100333‑M1/60‑1 (MQ=255)
                                 tggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAAt                                                                                                                      >  2:2277160‑M1/89‑150 (MQ=255)
                                 tggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAAt                                                                                                                      <  2:2598881‑M1/62‑1 (MQ=255)
                                 tggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAAt                                                                                                                      >  1:449866‑M1/89‑150 (MQ=255)
                                  ggctgcctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATc                                                                                                                     <  2:2123146‑M1/63‑1 (MQ=255)
                                        ctgacgccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCaa                                                                                                               >  1:2746093‑M1/82‑150 (MQ=255)
                                             gccagaagcattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAAcgcgc                                                                                                          <  2:2975503‑M1/74‑1 (MQ=255)
                                                   agcattggagcaccgtgcagtcgatgataagctgacaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGgag                                                                                                    <  2:3666381‑M1/80‑1 (MQ=255)
                                                     cattggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGagag                                                                                                  <  1:3741258‑M1/82‑1 (MQ=255)
                                                      attggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAgg                                                                                                 >  2:839295‑M1/68‑150 (MQ=255)
                                                      attggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAgg                                                                                                 >  1:4253625‑M1/68‑150 (MQ=255)
                                                       ttggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGc                                                                                                >  2:3826266‑M1/67‑150 (MQ=255)
                                                       ttggtgcaccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGc                                                                                                >  1:3968718‑M1/67‑150 (MQ=255)
                                                             caccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTg                                                                                          >  2:1872395‑M1/61‑150 (MQ=255)
                                                              accgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGc                                                                                         <  2:3902832‑M1/91‑1 (MQ=255)
                                                               ccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCg                                                                                        >  2:929984‑M1/59‑150 (MQ=255)
                                                               ccgtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCg                                                                                        <  1:2900313‑M1/92‑1 (MQ=255)
                                                                 gtgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTa                                                                                      <  2:195352‑M1/94‑1 (MQ=255)
                                                                  tgcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTAt                                                                                     >  2:552186‑M1/56‑150 (MQ=255)
                                                                   gcagtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTAtt                                                                                    >  2:2036593‑M1/55‑150 (MQ=255)
                                                                     agtcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTgg                                                                                  >  2:4324075‑M1/53‑150 (MQ=255)
                                                                       tcgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGc                                                                                >  2:2208900‑M1/51‑150 (MQ=255)
                                                                        cgatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGgcg                                                                               >  1:828950‑M1/50‑149 (MQ=255)
                                                                         gatgataagctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGggcgc                                                                              <  1:1175210‑M1/102‑3 (MQ=255)
                                                                                agctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtg                                                                       >  2:4062747‑M1/42‑150 (MQ=255)
                                                                                 gctgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGggtgg                                                                      >  2:2789002‑M1/41‑150 (MQ=255)
                                                                                   tgtcaaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGtt                                                                    <  1:3812434‑M1/112‑1 (MQ=255)
                                                                                      caaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGttttt                                                                 >  2:476254‑M1/36‑150 (MQ=255)
                                                                                       aaaccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttc                                                                <  2:3133100‑M1/116‑1 (MQ=255)
                                                                                          ccagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttt                                                             <  2:3052192‑M1/119‑1 (MQ=255)
                                                                                           cagatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttctttt                                                            >  1:2823088‑M1/31‑150 (MQ=255)
                                                                                            agatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttttc                                                           <  2:1296962‑M1/121‑1 (MQ=255)
                                                                                            agatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttttc                                                           <  2:350210‑M1/121‑1 (MQ=255)
                                                                                            agatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTttttcttttc                                                           >  1:3065586‑M1/30‑150 (MQ=255)
                                                                                             gatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCa                                                          >  1:3511395‑M1/29‑150 (MQ=255)
                                                                                             gatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCa                                                          <  2:698984‑M1/122‑1 (MQ=255)
                                                                                             gatcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCa                                                          <  1:2132781‑M1/122‑1 (MQ=255)
                                                                                              atcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGCGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAc                                                         <  1:1670390‑M1/123‑1 (MQ=255)
                                                                                              atcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAc                                                         >  1:3135655‑M1/28‑150 (MQ=255)
                                                                                              atcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAc                                                         <  2:3583653‑M1/123‑1 (MQ=255)
                                                                                              atcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAc                                                         >  1:1412161‑M1/28‑150 (MQ=255)
                                                                                               tcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAcc                                                        <  2:3982858‑M1/124‑1 (MQ=255)
                                                                                               tcaattcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTAGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCAcc                                                        <  1:4060345‑M1/124‑1 (MQ=255)
                                                                                                  attcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                     <  1:3655248‑M1/127‑1 (MQ=255)
                                                                                                  attcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                     >  2:1047298‑M1/24‑150 (MQ=255)
                                                                                                    tcgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                   <  1:388732‑M1/129‑1 (MQ=255)
                                                                                                     cgcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgag                                                  >  2:2694834‑M1/21‑150 (MQ=255)
                                                                                                      gcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                 >  1:3433340‑M1/20‑150 (MQ=255)
                                                                                                      gcgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                 >  2:1796610‑M1/20‑150 (MQ=255)
                                                                                                       gATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                >  1:4330758/1‑150 (MQ=255)
                                                                                                       cgctaacttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATAGGGCGCCAGGGTGGATTTTCTTTTCACCAGAGAGAc                                                <  2:1460446‑M1/132‑1 (MQ=255)
                                                                                                              ttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAc                                         <  1:2486552‑M1/139‑1 (MQ=255)
                                                                                                              ttacattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAc                                         <  1:658761‑M1/139‑1 (MQ=255)
                                                                                                                 cattaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCACCGCGCGGGGAGAGGCGGTTTGCGTCTTGGGCGCCAGGGTGGTTTTTCTTTTCGCCAGTGAGACGGGCAACAGc                                      <  2:1377836‑M1/142‑1 (MQ=255)
                                                                                                                  attaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCt                                     <  1:2061357‑M1/143‑1 (MQ=255)
                                                                                                                  attaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCt                                     <  2:4244195‑M1/143‑1 (MQ=255)
                                                                                                                   ttaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTg                                    >  1:1000565‑M1/7‑150 (MQ=255)
                                                                                                                   ttaattGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTg                                    <  1:2654718‑M1/144‑1 (MQ=255)
                                                                                                                      attGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt                                 <  2:2363658‑M1/147‑1 (MQ=255)
                                                                                                                      attGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt                                 <  1:1498369‑M1/147‑1 (MQ=255)
                                                                                                                      attGCGTTGCGCTCACTGCCCGCTTTCCAGTAGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTGTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt                                 <  1:4414875‑M1/147‑1 (MQ=255)
                                                                                                                       ttGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTg                                >  2:4493269‑M1/3‑150 (MQ=255)
                                                                                                                           gTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCt                            >  1:965623/1‑150 (MQ=255)
                                                                                                                            ttGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCtt                           <  1:3357547/150‑1 (MQ=255)
                                                                                                                               ctcTCACTGCGCGCTTTCCTGTCTGGTAACCTGTCGTGCCAGCTGCATTAATGAATCTGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCTCCAGTGAGACGGGCAACAGCTGATTGCCCTTCAc                        <  2:4412658/148‑1 (MQ=255)
                                                                                                                               cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCAc                        <  1:892525/150‑1 (MQ=255)
                                                                                                                                 cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCg                      <  1:2566327/150‑1 (MQ=255)
                                                                                                                                  tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGc                     >  1:574582/1‑150 (MQ=255)
                                                                                                                                    aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTAGTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                   >  1:1832045/1‑150 (MQ=255)
                                                                                                                                      tGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGAGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTgg                 >  2:3235419/1‑150 (MQ=255)
                                                                                                                                       gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGc                >  1:2984986/1‑150 (MQ=255)
                                                                                                                                       gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGc                >  1:2205184/1‑150 (MQ=255)
                                                                                                                                        cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGcc               >  2:1215257/1‑150 (MQ=255)
                                                                                                                                           gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTg            <  1:3203884/150‑1 (MQ=255)
                                                                                                                                            cTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTga           >  1:1742593/1‑150 (MQ=255)
                                                                                                                                             tttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgag          <  2:332035/150‑1 (MQ=255)
                                                                                                                                                     cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgc  >  1:1979865/1‑150 (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: