Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,157 T→G 6.3% intergenic (‑91/+504) lacZ ← / ← lacI pseudogene, truncated/DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273360,1570TG 6.3% 630.2 / 4.0 260intergenic (‑91/+504)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (72/162);  new base G (14/11);  total (86/174)
Fisher's exact test for biased strand distribution p-value = 1.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 272 total aligned reads displayed.

GGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGAC  >  CP009273/360008‑360306
                                                                                                                                                     |                                                                                                                                                     
ggCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                                                                                                                                       >  2:2550886/1‑150 (MQ=255)
ggCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                                                                                                                                       <  2:1836299/150‑1 (MQ=255)
ggCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                                                                                                                                       >  2:574547/1‑150 (MQ=255)
  cACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTc                                                                                                                                                     >  2:4080694/1‑150 (MQ=255)
    cTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGc                                                                                                                                                   >  2:2926343/1‑150 (MQ=255)
     ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                                  <  2:2653650/150‑1 (MQ=255)
     ctcaggagagcgtacaccgacacacaccagctaaaccgcaaggcccagtcgttcgacggagcctttcgttttatttgatgccgggggaattAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                                  >  2:3334428‑M1/92‑150 (MQ=255)
      tcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                                 <  2:69134‑M1/60‑1 (MQ=255)
      tcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                                 <  2:1384988‑M1/60‑1 (MQ=255)
       cACCGATTGCCAGCGGCTTACCATCCAGCGCAACCATCCGGTGCAGGAGCTCCGATTATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                                                                                <  2:1352172/150‑1 (MQ=38)
        aggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgggttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACg                                                                                                                                               <  2:4305892‑M1/62‑1 (MQ=255)
        aggagagcgtacaccgacaaccaacagatcaaacgaaagccccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                               >  2:1098377‑M1/89‑150 (MQ=255)
         ggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgccgggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCTACGt                                                                                                                                              <  2:4624045‑M1/63‑1 (MQ=255)
            gagcttacaccgacaaacaacagaaaaacccaacggcccagtctttcgcctgagcctttcgttttttttgatgcgtgggtaattAATTCCTGGCAGTTTATGGCGGGCGACCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                                           >  1:1563242‑M1/85‑150 (MQ=255)
            cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                                           <  2:490589/150‑1 (MQ=255)
             ttGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGAGCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                          <  2:710372/150‑1 (MQ=255)
              tGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAACTTCaaa                                                                                                                                         >  1:3863275/1‑150 (MQ=255)
               cgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTGGCAACGGTCAAAt                                                                                                                                        <  2:2757038‑M1/69‑1 (MQ=255)
                gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                                       >  2:2140364‑M1/81‑150 (MQ=255)
                gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                                       <  2:1421551‑M1/70‑1 (MQ=255)
                 cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAATcc                                                                                                                                      >  2:73479/1‑150 (MQ=255)
                 cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCAACCTTCAAAtac                                                                                                                                      >  1:2261695/1‑148 (MQ=255)
                 cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCTACGATGAAAtac                                                                                                                                      >  1:1522143/1‑148 (MQ=255)
                  aGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGTTGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTGGCAACGTGGAAATTCg                                                                                                                                     >  1:2154450/1‑150 (MQ=37)
                     ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTc                                                                                                                                  >  1:2297450/1‑150 (MQ=255)
                      cgacaaacaccagataaaacgaaaggcccagtctttcgactgagcctttcgttttattttatgcctggggaattAATTCCTGGAATTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTcc                                                                                                                                 >  2:234814‑M1/75‑150 (MQ=255)
                        ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                               >  2:3508726/1‑150 (MQ=255)
                              aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                         >  2:3775031‑M1/67‑150 (MQ=255)
                               cccgcgcCACCATCCAGTGCAGGAGCTCCGATTATAAGCTGTCAAACATGACAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                        <  1:3882723/147‑1 (MQ=255)
                               ccacccccacaaTCGTGTGCTGGAGCTCCGATGATACGCTGTAAACCATGTGAATTCGTACCCCCACGTCCTGGCTGTTTTTGGGGGGGGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                        <  1:3732792/144‑1 (MQ=12)
                               ccAGCGCCACCATCCAGGGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATACCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                                                                                                        <  1:3118534/150‑1 (MQ=255)
                                ccgctaaaacgaccggcccagtctttcgactgagcctttcgttttattggatgcctggggaattAATTCCGGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                       >  2:3997008‑M1/65‑150 (MQ=255)
                                     aaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                                  <  2:1520344‑M1/91‑1 (MQ=255)
                                      caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTAGCAACCGTAAAAAACGCTCCCCGCGGCTTTGTACTa                                                                                                                 >  1:2624810/1‑150 (MQ=25)
                                      caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGGGTCGCAACGGTCAAATACGCGCCCCGCGGCTTTGtacga                                                                                                                 >  1:249283/1‑146 (MQ=25)
                                       accaGCCAGTGCTGGAGCTCCTTTTATACGCTGTAAAACATGATAATTCCGATCCCCATTGCGGGGCTGTTTTTGGGGGGCGTCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                <  1:1309241/150‑1 (MQ=17)
                                       aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                <  1:3123765‑M1/93‑1 (MQ=255)
                                         caTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                              <  2:186887/150‑1 (MQ=255)
                                           tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                            >  2:3623854/1‑150 (MQ=255)
                                              aGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                         <  2:4029978/150‑1 (MQ=255)
                                               gTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTAGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                        <  2:1265776/150‑1 (MQ=255)
                                                 ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                      <  2:2532333‑M1/103‑1 (MQ=255)
                                                  ctggACCTCCGATTATAAGCTGACAAAAATGATATTTCCTACCCCCAAGTCCTGGATTTTTATGGCGGGCGACCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                     <  1:4445303/148‑1 (MQ=14)
                                                  cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                     <  2:4401135/150‑1 (MQ=255)
                                                   aGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAAAGTTAAAATACGCTCCCGGCGGATTTGTCCAACTCAGGAGAGCGt                                                                                                    >  1:4343415/1‑150 (MQ=21)
                                                    gtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                   <  2:3050368‑M1/106‑1 (MQ=255)
                                                    gtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                   <  2:2491485‑M1/106‑1 (MQ=255)
                                                     tttttcgagtgagcatttgtttttatttgttgcctggggaattATTTCCTGGATGTTTTTGGGGGGGGTCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                  >  1:683177‑M1/44‑150 (MQ=255)
                                                     tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGGTCCCGGCGGGTTTGTCCTACTCAGGTGagcggtc                                                                                                  <  2:2432146‑M1/107‑4 (MQ=255)
                                                        cgccgcagatACGCTGTCAAACATGAGACCTCCGATCCCCAATACCCGGCAGTTTATGCCGGGCGTCCTGCCCGCCACCCCCCGGGCCGTTGCTTCGCAACGTTCCCATCCGCTCCCGGCGGATTTGTCCTCCTCAGGAGAGCGGTcccc                                                                                               <  2:3670097/143‑1 (MQ=17)
                                                        cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                               <  2:3671683/150‑1 (MQ=255)
                                                         tCCGATGATACGCTGTCAAACATGAGAAGTCCGATCCCCAAGTCCGGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCg                                                                                              <  2:139406/150‑1 (MQ=17)
                                                          tgcctgagcctttcgttttatttgatgtcagggtaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGa                                                                                             >  2:1121727‑M1/39‑150 (MQ=255)
                                                          cgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCGCCCCGCGGATTTGTCCAACACAGGAGAGCGTTAACCGa                                                                                             <  1:3518241‑M1/112‑1 (MQ=255)
                                                           gactgagcctttcgttttatttgatgcctggggaattCATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                            >  2:2674886‑M1/38‑150 (MQ=255)
                                                               gttAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                        <  1:1826919/148‑1 (MQ=255)
                                                                 tAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                      >  2:2718382/1‑150 (MQ=255)
                                                                 tAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                      >  2:3900612/1‑150 (MQ=255)
                                                                 tAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                      <  2:3649698/150‑1 (MQ=255)
                                                                 tAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAAAACCCACCCCGCGGCTTTGTCCTACACAGGAGAGCGGTCACCGccaaacca                                                                                      >  1:296312/1‑148 (MQ=18)
                                                                 gcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACcgacaaacca                                                                                      <  2:4735104‑M1/119‑3 (MQ=255)
                                                                  aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTACTAATCATGAGAGCGTTCACCGACAaacaac                                                                                     >  2:2267420/1‑150 (MQ=255)
                                                                    gCTGTCAACCATGCGAATTCCGATCCCCAATTCCTGGCTGTTTTTGGCGGGCGGCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                   <  1:1026517/150‑1 (MQ=17)
                                                                    gCTGTCAAACCTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                   <  1:3169804/150‑1 (MQ=255)
                                                                    gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                   <  2:2586486/150‑1 (MQ=255)
                                                                      tttCAAACATGAGAATTCCGACCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                 <  1:3003757/148‑1 (MQ=255)
                                                                        tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                               <  2:1941285/150‑1 (MQ=255)
                                                                        tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                               <  2:886846/150‑1 (MQ=255)
                                                                         ttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                      >  2:1780776‑M1/24‑78 (MQ=255)
                                                                           ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                            <  2:4708555‑M1/129‑1 (MQ=255)
                                                                           aaCATGAGAATTCCTCACCCCATTTCGTGTCTGTTTTTGGGGGGGGTCCGCCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTCGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                            <  1:506885/150‑1 (MQ=11)
                                                                            tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAATACCCTCCCCGCGGATTTGTCCTACACAGGGGAGCGTTGACCGACAAACCACACATAAAacc                                                                           <  1:4704589‑M1/130‑2 (MQ=255)
                                                                            tatttgatgcctggggaattAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                           >  1:2757038‑M1/21‑150 (MQ=255)
                                                                              ttttttgcgtgggtaattAATGCTGGGCTGTTTTGGGGGGGCGCCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                         >  1:94623‑M1/19‑150 (MQ=255)
                                                                                tgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                       <  2:3438389‑M1/134‑1 (MQ=255)
                                                                                tgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATACGCTCCCGGCGGGTTTGTGCTACGCAGGAGTGCGTTCACCGTCAAAAAACAGATTAATCGAAAg                                                                       <  2:1563242‑M1/134‑1 (MQ=255)
                                                                                   aaTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                    <  2:3863275/150‑1 (MQ=255)
                                                                                        ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCGTCAAAACCGCACCCCGCGGCTTTGTCCAACACAGGAGAGCGTTCACCGACAAACAACACATAAAACGAAAGGcccactc                                                               <  2:1689630‑M1/142‑4 (MQ=255)
                                                                                         gATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                              <  2:2998649/150‑1 (MQ=38)
                                                                                          aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                             <  2:1883290/150‑1 (MQ=38)
                                                                                            aattATTTCGTGGCTGTTTTTGGGGGGCGTCGTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                           >  1:2550886‑M1/5‑150 (MQ=255)
                                                                                            aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                           >  1:1102366‑M1/5‑150 (MQ=255)
                                                                                             cccAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                          <  2:3635631/150‑1 (MQ=17)
                                                                                               taaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGTTAAAAAACGCTCCCGGCGGAATTTGCCAACTCCGGCGAGGGGTGACCGGCAACCACCAGGTAAAACGAATGGCCCAGTCTTTCGAc                                                        >  1:2106000/2‑150 (MQ=2)
                                                                                               tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATGCGCTCCCGGCGGCGTTGTCCTACTCCGGAGAGCGGTCACCGACACACAACCGATACATCGACAGGCCCCGTCTGTCgag                                                        <  2:3124403‑M1/149‑2 (MQ=255)
                                                                                                aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                       >  2:785902/1‑150 (MQ=32)
                                                                                                  ttCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                     <  1:1928530/150‑1 (MQ=14)
                                                                                                  ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                     <  2:659810/150‑1 (MQ=14)
                                                                                                             ttttGGCGGGCGTCCGCCCCCCCACCCTCCGGGCCGTTGATTCGCAACGTTCAAATCCGCTCCCGGCGGATTGGTCCTACTCAGTAGAGCGTTCACCGACAAAAAACATATAAAACTAAAGGCCCAGTCTTTCGACTGAGCCTTTCGttt                                          <  1:3641233/147‑1 (MQ=1)
                                                                                                                             gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCAGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                          <  1:3920104/150‑1 (MQ=14)
                                                                                                                              cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                         <  1:1051768/150‑1 (MQ=32)
                                                                                                                              cccGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                         <  2:2932873/150‑1 (MQ=17)
                                                                                                                               ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCACGAGAGCGTACACCGACAAACAACCGATAAAACGAAAGGCCCAGTCTGTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggac                        <  1:3454151/150‑2 (MQ=11)
                                                                                                                               ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                        <  2:2297450/150‑1 (MQ=32)
                                                                                                                                 cccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                      <  1:1234070/149‑1 (MQ=21)
                                                                                                                                            ggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaatt           <  1:4755440/150‑1 (MQ=255)
                                                                                                                                            gccggTTGCTTAGAATCGTTCAATGCCGCTCGCGGCGGATTTGTCCTACTCATGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaatt           <  1:3412028/146‑1 (MQ=14)
                                                                                                                                                gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTg       >  2:1198543/1‑150 (MQ=255)
                                                                                                                                                   gCTTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTACTCAGTAGAGCGTTCACCGAAAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAg    <  1:1692351/150‑1 (MQ=255)
                                                                                                                                                     ttCGCAACGTTCACATCCGATCCCGGCGGTTGTGTCCTACTCAGGAGAGCGTTCCCCGACCAACAACAGATACAACGACAGGCCCAGTCTGTCGACTGAGCCTTTAGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGAc  <  1:4180931/150‑1 (MQ=17)
                                                                                                                                                     |                                                                                                                                                     
GGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGAC  >  CP009273/360008‑360306

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