Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,157 T→G 9.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 9.3% 499.0 / 11.0 244intergenic (‑91/+504)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (45/166);  new base G (21/9);  total (67/177)
Fisher's exact test for biased strand distribution p-value = 2.08e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

GGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGAC  >  CP009273/360008‑360306
                                                                                                                                                     |                                                                                                                                                      
 gCACTTCACCGCTGGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                       <  2:4293535/150‑1 (MQ=255)
   tactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCg                                                                                                                                                     <  2:2221220‑M1/57‑1 (MQ=255)
      tCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCaa                                                                                                                                                 >  1:868165/1‑150 (MQ=255)
       caggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                                                                                <  1:2497959‑M1/61‑1 (MQ=255)
       aaacGCTTGACAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAATCTGTAAAACATGAGAATTCCGATCCCCAATTCCTCGAAGTTTATGGCGGGCGTCCTGCCAGCCACCCTCCGGGCCGTTTCATAGCaaa                                                                                                                                                >  2:4057193/4‑149 (MQ=38)
        aCCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                               >  2:4254167/1‑150 (MQ=255)
          cGCTTGCCCGCGGCTTACCAGCCAGCGCCACCAACCCGTGCAGGAGCTCCGATGATAAGCTTTAAAACATGATAATTCCGATCCCAATTTCCTGGCATTTTTGGGCGGGCGCCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                             <  1:3862936/150‑1 (MQ=25)
           gCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCCAACGtc                                                                                                                                             >  1:3388454/1‑149 (MQ=255)
           gCTTGCCAGCGGCTTACCACCCCGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGAAAAACATGAGAATTCCGATCCCCAATCCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                            <  2:2545200/150‑1 (MQ=37)
            cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCa                                                                                                                                           >  1:3237356/1‑150 (MQ=255)
             ttGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                          <  2:2771984/150‑1 (MQ=255)
               gCCAGGGGCTTACCATCCCGCGCCACCATCCCGGGCTGGAGCTCCGTGGATACGCTGTCAAACATGAGAATGCCGATCCCCAAGTCCTGGATTTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                        <  1:2370998/150‑1 (MQ=21)
               gCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                        <  2:3983422/150‑1 (MQ=255)
                 cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTGGCGTCGCAACCATCAAATcc                                                                                                                                      >  1:4216966/1‑150 (MQ=38)
                 cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                      >  1:3586821/1‑150 (MQ=255)
                 cAGCGGATTACCATCCCGCGCCACCCTCCTGTGCTGGAGCTCCGATTATAAGGTGACAAACATCAGAATCCCGACCCCCACTGCCGGGCTGTTTATGGGGGGCGACCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                      <  1:945466/150‑1 (MQ=18)
                  aGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                     <  2:1695711/150‑1 (MQ=255)
                    cGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                   <  2:3071305/150‑1 (MQ=255)
                     ccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGGTTCGCcacctgcaaaacaggaa                                                                                                                                  <  1:3672537‑M1/75‑8 (MQ=255)
                      cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTAAAATCCGCTcc                                                                                                                                 <  2:1453350‑M1/76‑1 (MQ=255)
                       cTTACCATCCAGCGCCACCATCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                                <  2:3889060/150‑1 (MQ=255)
                         caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTGAGCAACGTTCAaatacgcagcccg                                                                                                                              <  1:1181588‑M1/79‑7 (MQ=255)
                         caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCCACGTTCAAATCCGCTCCcgg                                                                                                                              <  2:3533864‑M1/79‑1 (MQ=255)
                         caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGCTTGGCCACCTTCAAAACCGCTCCcgg                                                                                                                              <  2:67226‑M1/79‑1 (MQ=255)
                             aTCCCGCGCCCCCAGCCTGTGCTGTAGCTCCGATGATAAGCTTTCAAACATGACATTTCCGACCCCCATTGCGTGGATGTTTTGGGGGGGGGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                          <  1:1868303/150‑1 (MQ=12)
                             aTCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                          <  2:4134265/150‑1 (MQ=255)
                              tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGATAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                         <  1:2816749/150‑1 (MQ=255)
                              tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                         <  2:4211919/150‑1 (MQ=255)
                              aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                         <  2:1766945‑M1/84‑1 (MQ=255)
                                      caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTTCGCAACCGTCAAATAAGCACCCCGCGGCtttttgcaa                                                                                                                 >  1:2679555/1‑144 (MQ=25)
                                      aaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                 <  1:3891289‑M1/92‑1 (MQ=255)
                                       accaTCCAGTGCTGGAGCTCCGATGATAAGCTTTAAAACATGTTAATTCCGATCCCCAATTCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                                <  1:4297577/150‑1 (MQ=21)
                                          gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattATTTCCTGTCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTCAAATCCGCTCCCGGCGGCTTTGTCCTCCTCa                                                                                                             <  2:798844‑M1/96‑1 (MQ=255)
                                          gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGGTCAAATGCGCACCCGGCGGATTTGTGCGACtcc                                                                                                             <  2:2972219‑M1/96‑2 (MQ=255)
                                          aTCCTGTGCTGGAGCTCCGATGATATGCTGACAAACATGAGAATTCCGACCCCCAATTCCTGGCAGTTTTGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                             <  1:2318901/150‑1 (MQ=18)
                                              ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                     >  2:2499876‑M1/51‑122 (MQ=255)
                                                  cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                     >  2:3012410/1‑150 (MQ=255)
                                                   agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                    <  2:4260219‑M1/105‑1 (MQ=255)
                                                   aGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                    <  2:3544064/150‑1 (MQ=255)
                                                    ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                   <  2:1972622/150‑1 (MQ=255)
                                                    ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                   <  2:4162673/150‑1 (MQ=255)
                                                      ctttcgactgagcctttggttttatttgatgcgtggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                                 >  1:2134028‑M1/43‑150 (MQ=255)
                                                      ctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGGTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                                 <  2:2306637‑M1/108‑1 (MQ=255)
                                                      aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                                 <  2:3813323/150‑1 (MQ=255)
                                                       gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCATTTTCTGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                                <  1:1322095/150‑1 (MQ=255)
                                                            gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTGGGAACCGTTAAATACCCCACCGGCGGATTTGTGCTAAGCCGGAGTGCGGTAACCGaca                                                                                           >  1:433083/1‑150 (MQ=18)
                                                             ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                                                                                                                                        <  1:4268958‑M1/54‑1 (MQ=255)
                                                             ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                                                                                                                                        >  2:4268958‑M1/36‑89 (MQ=255)
                                                              tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCAACCTTAAAATACGGACCCGGCGGCTTTTTGGTAAACAGGTGGGCGTTAACACacaaa                                                                                         >  2:370003/1‑150 (MQ=18)
                                                               gatACGCTTTAAAACATGATAAGTCCGATCCCAAATTCCTGGCTGTTTAGGGGGGGCGCCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                        <  1:457494/150‑1 (MQ=14)
                                                                agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                       <  2:2504106‑M1/118‑1 (MQ=255)
                                                                 tAAGCTTTAAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                      <  1:4138830/150‑1 (MQ=255)
                                                                  cctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCAAaacaac                                                                                     <  2:1570209‑M1/120‑1 (MQ=255)
                                                                    gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                   <  2:92621/150‑1 (MQ=255)
                                                                    gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                   <  1:1674890/150‑1 (MQ=255)
                                                                        tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                               <  2:3915038/150‑1 (MQ=255)
                                                                        aaaaaaATCAGAATTCCTACCCCAAATACCTGGAATTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                               <  1:3059258/148‑1 (MQ=11)
                                                                         cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTAAAATCCGCACCCCGCGGATTTGTCCTACTCAGGAGATCGTTAAACGACACACAACAGATaaa                                                                              >  2:3641361/1‑150 (MQ=17)
                                                                           aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                            <  2:242535/150‑1 (MQ=255)
                                                                           aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                            <  2:2671572/150‑1 (MQ=255)
                                                                           aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAAACTTAAAAAACGCGACCCGCGGAGTTTTGCAACTAAGGAGAGCGTTAAACGACCAACAACACAAAAAAc                                                                            >  1:2253553/1‑150 (MQ=12)
                                                                             cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGCGCGTTCACCGACAAACAACACATAAAACGa                                                                          >  2:1523453/1‑150 (MQ=255)
                                                                             atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGTATTGTCCTACTCAGGGGAGCGTTGACCGACACACAACAGTTAAAACGa                                                                          <  2:3736905‑M1/131‑1 (MQ=255)
                                                                                tgatgcctgggtaattAATTCGTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                       >  2:3657642‑M1/17‑150 (MQ=255)
                                                                                  gaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCCCAACGTTCAAATCCCCTCCCGGCGGATTTGTCCTACGCAGGCGCGCGGTCTCCGCCAAACCCCAGATAAGAAGAAAGGc                                                                     >  1:2860867/1‑150 (MQ=14)
                                                                                  atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCCGCGGCTTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGAAAGGc                                                                     <  2:3372534‑M1/136‑1 (MQ=255)
                                                                                     ttCCGATCCCCAATTCCTGGCAGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                  <  2:564014/150‑1 (MQ=255)
                                                                                      ctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                                 <  2:311023‑M1/140‑1 (MQ=255)
                                                                                       tgggtaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                                           >  1:3506708‑M1/10‑75 (MQ=255)
                                                                                       tggggaattAATCCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGCTTTGTCCTAATCAGGAGAGGGTTAACCGACAAACAACAGATAAAAAGAAAGGCCCAGt                                                                <  1:2463217‑M1/141‑1 (MQ=255)
                                                                                        ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                >  1:1051097‑M1/9‑133 (MQ=255)
                                                                                        ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                                <  2:1051097‑M1/125‑1 (MQ=255)
                                                                                         gggaattAATTCCTGGCAGTTTATGGTGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                              >  1:3978897‑M1/8‑150 (MQ=255)
                                                                                         gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                              >  1:1766945‑M1/8‑150 (MQ=255)
                                                                                         gATCCCCAATTCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                              <  2:1105129/150‑1 (MQ=255)
                                                                                          aTCCCCATTTCGTGGCTGTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                             <  1:2666175/150‑1 (MQ=11)
                                                                                          aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                             <  2:2155885/150‑1 (MQ=38)
                                                                                          aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCCACGTTCAAATCCGCTCCCGGCGGCTTTGTGCTAATCAGGGGGGGGGTCCCCGACCAACAACAGATTAAACGGAAGGCCCCGTCtt                                                             >  2:3679632/1‑150 (MQ=11)
                                                                                            ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                           <  1:3563201/150‑1 (MQ=35)
                                                                                            aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGGTAACCGACAAACAACAGATAAAAACAAAAGCCCACTcttta                                                           <  2:2563305‑M1/146‑2 (MQ=255)
                                                                                               tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGGTCAAATCCGCGCCCCGCGGCTTTGGCCTACACAGGAGAGCGTTCACCCACCAACAACAGATAAAACGAAAGGCCCAGTCATTCGAc                                                        <  2:418177‑M1/149‑1 (MQ=255)
                                                                                               cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                        <  1:2882341/150‑1 (MQ=33)
                                                                                                atttCCTGGCTTTTTCTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                       <  1:843478/148‑1 (MQ=9)
                                                                                                aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                       <  2:1819434/150‑1 (MQ=32)
                                                                                                aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAAAGGTCAAAACCGCGCCCGGCGGGTTTGTCCTACACAGGAGGGCGTTCACCGCCCAACCACAGATAAAACGCAAGGGCCAGTCTTTCGTCt                                                       >  1:1999055/1‑150 (MQ=11)
                                                                                                 aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                      >  2:2412019/1‑150 (MQ=32)
                                                                                                 aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                      >  2:532174/1‑150 (MQ=32)
                                                                                                          tgttttgggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCg                                             <  1:4174501/143‑1 (MQ=1)
                                                                                                                            tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                           <  1:2418664/150‑1 (MQ=11)
                                                                                                                             gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAGCCGCTCCCGGCGGATTGGTCCTCCTCAGGAGAGCGTTCACCGACAAACAACAGATAAACCGAAAGGCCCAGTCTTTCGACTGAGCCGTTCGTTTTATTTGATGCCTgggg                          <  2:2931970/150‑1 (MQ=11)
                                                                                                                               cccccACCCTCCGGGCCGGTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                        <  1:2718110/147‑1 (MQ=12)
                                                                                                                                cGCCACCCGCCGGGCCGTTGCATCGCAACGTACAACTCCGCGCCCGGCGGATTGGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGATAACACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                       <  2:1903363/150‑1 (MQ=11)
                                                                                                                                 gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGACTGGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATACAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                      <  2:2094991/150‑1 (MQ=11)
                                                                                                                                 gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                      <  1:3453626/150‑1 (MQ=34)
                                                                                                                                 gCCACCCGCCGGGCCGTTGCTTCGCAACGTTCACAGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGATAACACGAAAGGCCCAGGCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                      <  2:4216966/150‑1 (MQ=11)
                                                                                                                                 cccACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                      <  1:3293328/149‑1 (MQ=21)
                                                                                                                                     cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggga                  <  1:2979274/150‑1 (MQ=38)
                                                                                                                                                  tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                     >  1:3339450/1‑134 (MQ=35)
                                                                                                                                                  tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                     <  2:3339285/134‑1 (MQ=35)
                                                                                                                                                     ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGAATTAATTCTTGTAGAc  >  1:292809/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
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: