Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,157 T→G 8.0% 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 8.0% 543.7 / 6.2 253intergenic (‑91/+504)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (68/157);  new base G (18/10);  total (86/167)
Fisher's exact test for biased strand distribution p-value = 5.80e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

GCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGA  >  CP009273/360009‑360305
                                                                                                                                                    |                                                                                                                                                    
gCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                                                                                                                                      <  2:2674458/149‑1 (MQ=255)
 cACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGTTAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                     <  1:2774961/149‑1 (MQ=255)
  aCTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTc                                                                                                                                                    >  2:483636/1‑149 (MQ=255)
  aCTTCACCGCTTGCCAGAGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTATATGGCTGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTc                                                                                                                                                    >  1:3506522/1‑149 (MQ=255)
    ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGc                                                                                                                                                  >  2:4449096/1‑149 (MQ=255)
     tCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTTGCa                                                                                                                                                 >  1:2036065/1‑149 (MQ=255)
      ccccGCTTGCCAGCGGCTTACCATCCAGCGCCACCAGCCAGTGCAGCAGCTCCGATGATACGCTGACAAACATGATAATTCCGATCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                                <  1:3485616/147‑1 (MQ=37)
      caggagagcgttcaccgacaaacaacagataaaacgaaacgcccagtctttcgactgagcctttcgttttattttttgcctggggaattAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                                >  2:2384694‑M1/90‑149 (MQ=255)
      caggagagcgtacaccgacaaacaacagataaaacgacaggcccagtctttcgacttagcatttctttttatttgttgcctggggaattATTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                                >  2:1761631‑M1/90‑149 (MQ=255)
      cACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                                >  1:2488035/1‑149 (MQ=255)
       aggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCACCGGGCCGTGGCTGAGCAAc                                                                                                                                               <  2:2004893‑M1/61‑1 (MQ=255)
         cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                                                                                             <  2:2971412/149‑1 (MQ=255)
         cGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                                                                                             >  2:3447057/1‑149 (MQ=255)
         cGCTTCCCCGCGGCTTACCACCCCGCGCCACCACCCAGGGCCGGAGCTCCGATGATAAGCTGACAAACAGGAGAATTCCGATCCCCAATCCCGGGCTGTTTTTGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                                                                                             <  1:3843822/149‑1 (MQ=25)
          agagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCATGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                            <  2:906276‑M1/64‑1 (MQ=255)
           tagcgtaccccgacaaactacagatacaacgaccggcccagtcgttcgtctgagcctttcgtttaatgtgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                           >  2:4327913‑M1/85‑149 (MQ=255)
           cTTGCCAGCGGCATCCACTCCAGCCCCACCATCCCTTGACGGTGCTCATATGATAAGTTGTCAAAAATCAGATTTCGGATCCCATATTCCTTGATGTTTATTGCGGGGGTCATGCCAGCCTCCCTCCGGGCCGTTGCTTCGCAACGTTc                                                                                                                                           <  1:1130314/149‑1 (MQ=12)
               ccAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                       <  2:1950526/149‑1 (MQ=255)
                cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                                      <  2:1562215/149‑1 (MQ=255)
                 tcaccgacaaccaacagcaacaacgaaaggcccagtctttcgaatgagccgttgtttttatttgatgcctggggaattAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                     >  2:3284562‑M1/79‑149 (MQ=255)
                   cGGCTTACCAGCCCGCGCCACCATCCAGTGCTGGAGCTCCGATTATAAGCTGTAAAACATGAGAATTCCGACCCCCAATTCCTGGATGTTTTGGGCGGGCGTCGTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                                   <  1:2918511/149‑1 (MQ=25)
                   acccacaaacaaccgctaaaacgacaggcccaggctttcgactgagcctttcgatttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                                   >  2:3097163‑M1/77‑149 (MQ=255)
                    ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                  >  1:580707/1‑149 (MQ=255)
                    ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCCACCTTCAAAACCGcg                                                                                                                                  >  1:4678573/1‑148 (MQ=255)
                    ccgacaaacaacagataaaacgaaaggccctgtctttcgactgagcctttctttttatttgatgccgggggaattAATTCCTGGCATTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                  >  1:3534444‑M1/76‑149 (MQ=255)
                      gacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTcc                                                                                                                                <  2:262341‑M1/76‑1 (MQ=255)
                      cTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTcc                                                                                                                                <  2:1938658/149‑1 (MQ=255)
                       acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                               <  2:4611686‑M1/77‑1 (MQ=255)
                       acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                               <  1:2722865‑M1/77‑1 (MQ=255)
                            caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                          <  2:274249‑M1/82‑1 (MQ=255)
                            caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttagttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                          <  1:4351098‑M1/82‑1 (MQ=255)
                            caacagataaaacgaaaggcccagtctttcgactgagcctttcgttttagttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                          <  1:3192033‑M1/82‑1 (MQ=255)
                             tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCACCcggcggc                                                                                                                         >  1:4381440/1‑148 (MQ=255)
                                aGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                      >  1:2538986/1‑149 (MQ=255)
                                   gcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACCTGAGAATTCCGATCCCCCATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                   <  2:2517679/149‑1 (MQ=255)
                                   gcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                   >  2:1982526/1‑149 (MQ=255)
                                      aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                                                                                                <  2:3846510‑M1/92‑1 (MQ=255)
                                       acgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCGTAGCAACCTTCAAAACCGCGACCCGCGGATTtgtccgac                                                                                                               <  1:1709586‑M1/93‑4 (MQ=255)
                                          acctgggCAGGAGCTCCGATTATAAGCTGACAAAAATTACAATTCCGACCCCCACTTCCTGGCAGTTTTTGGGGGGCGTCGTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                            <  1:2816468/143‑1 (MQ=17)
                                           ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                           <  2:1495115/149‑1 (MQ=255)
                                                ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                      <  2:1016148‑M1/102‑1 (MQ=255)
                                                   ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                   <  2:4543672/149‑1 (MQ=255)
                                                    tctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGCGAGCGtt                                                                                                  <  2:2826236‑M1/106‑1 (MQ=255)
                                                    gagttcgactgagactttctttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                  >  2:4678568‑M1/44‑149 (MQ=255)
                                                    gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                  >  2:2293967/1‑149 (MQ=255)
                                                     ctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGTCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTc                                                                                                 <  2:484755‑M1/107‑1 (MQ=255)
                                                      gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCACCGGGCAAATCCGCGCCCCGCGGCTTTGTCCTACTCAGGAGAGCGtgag                                                                                                >  1:881977/1‑146 (MQ=18)
                                                          gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTCCTCAGGAGAGCGTTCACCGa                                                                                            <  2:1174997‑M1/112‑1 (MQ=255)
                                                          cGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCGTAAAATCCGCAACCCGCCGCTTTGTCCAACTCAGGTGAGCCTTCACCGa                                                                                            >  1:4094301/1‑149 (MQ=18)
                                                           gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                           <  1:2618312/149‑1 (MQ=255)
                                                            atgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGGTTGGCAAAGTTAAAAACACCACCCCCCGGGTTTGTGCTACACAGGGGAGAGGTAAGCGaca                                                                                          >  1:2183251/1‑149 (MQ=18)
                                                             tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                         <  2:779842‑M1/115‑1 (MQ=255)
                                                             tgagcctttcgttttatttgatgcctggggaattAATACATGGCAGTTTATGTCGGGCGTCCTGCCCGACACCCACCGGTCCGTTGCTTCTCAACCTTCAAATCCGCTCACCGCGGATTTGTCCTACTCAGGAGAGCGGTCAccgaaaa                                                                                         <  2:3578532‑M1/115‑4 (MQ=255)
                                                              gagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACCTTAAAATCCGCTACCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                        <  2:2227489‑M1/116‑1 (MQ=255)
                                                              gagcctttcgatctacatgatgcccggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCATCCGGTCCGTTGCTTAGCAACTATCAAATCCGCTCACGTCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaa                                                                                        >  2:3530260‑M1/34‑149 (MQ=255)
                                                               atAAGCTGTCAAACATGAGAATTCCGATCCCCAATGCCTGGCAGTTTCTGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                       <  2:2756262/149‑1 (MQ=255)
                                                                  aGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAaacaac                                                                                    <  2:4105121/149‑1 (MQ=255)
                                                                   gCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                   <  2:4427850/149‑1 (MQ=255)
                                                                     tGTCAAACATGAGAAGTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                                 <  1:1862134/149‑1 (MQ=255)
                                                                          aaCATGTGATTTCCTATCACCAATTCGTGGGAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTGCTACTCAGGAGAGCGTTCACCTAGATACAAAAGATataa                                                                            <  2:2912721/149‑1 (MQ=14)
                                                                           tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                           <  2:519941‑M1/129‑1 (MQ=255)
                                                                           aCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAc                                                                           <  2:3083804/149‑1 (MQ=255)
                                                                            atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAAAACCCAACCGGCGGATTTGTCCTACACAGGAGATCGTTCACCCCCAAACAAAAGataaaaac                                                                          <  1:4641480‑M1/130‑3 (MQ=255)
                                                                              tGAGAATTCCGATCCCCAATTCGTGGGTGTTTTTGGCGGGCGTCCTGCCCGCCTCCCTCCGGGCCGTTGCTTCGCATCGTTCAAATCCGCTCCCGGCGGTTTTGTCCTACTCAGGAGAGCGTTCACCGACATACAACAGATAAAACGaa                                                                        <  2:2880319/149‑1 (MQ=14)
                                                                              tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                        <  2:850524/149‑1 (MQ=255)
                                                                              tGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAACTTCAAATACGCTCCCCGCGGATTTTTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGaa                                                                        >  1:1020639/1‑149 (MQ=14)
                                                                               gagaATTCCGATCCCCATTTCGTGGCTGTTTTTGGCGGGGGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                       <  1:1786042/149‑1 (MQ=14)
                                                                               gagaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaaa                                                                       <  2:4245510/149‑1 (MQ=255)
                                                                                   aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                   >  2:3504501/1‑149 (MQ=255)
                                                                                    ttCCGATCCCGAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                  <  2:2538986/149‑1 (MQ=255)
                                                                                    ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                  <  2:2709960/149‑1 (MQ=255)
                                                                                      ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                                >  2:86147/1‑149 (MQ=40)
                                                                                         aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                             <  2:3590867/149‑1 (MQ=38)
                                                                                         aTCCCAAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                             <  1:2985446/149‑1 (MQ=25)
                                                                                          gaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATACCCTCCCGGCGGATTTTTCCTACTCAGGCGAGCGTTCACCCACAAACAACACATAAAAACAAATGGCCACTCtt                                                            <  1:3835019‑M1/144‑1 (MQ=255)
                                                                                           caccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCttt                                                           <  1:2293967/147‑1 (MQ=21)
                                                                                             ttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                         <  2:3951133‑M1/147‑1 (MQ=255)
                                                                                             ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTACGAATCAGGAGAGCGTTGACCGACAAACAACAGATAAAACGAAAGGCCCCGTCgtttg                                                         >  1:1179310/1‑147 (MQ=11)
                                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAATGCCCCGTCTTTCGAc                                                       >  1:1550023/1‑149 (MQ=11)
                                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGGGGATTTGTCCGACTCAGGAGAGCGTTCACCGGCAAACAACAGATAAAACGAAAGGCCCAGTGTTTGGAc                                                       >  2:2107195/1‑149 (MQ=11)
                                                                                                aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTGGGCAAAGGTGCAATCCGCTCCCCGGGGCTTTGTACTGATCAGGTGTGCGGTGCCCGACAAACAACAGGTAAAACGAGTGGCCCAGTGGTTCGGCt                                                      >  1:1378236/1‑149 (MQ=2)
                                                                                                  tCCTGGCATTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                    <  1:124110/149‑1 (MQ=11)
                                                                                                  tCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                    <  1:4774920/149‑1 (MQ=11)
                                                                                                  gcgtGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                    <  1:838618/146‑1 (MQ=0)
                                                                                                                        gccggCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTACAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGcc                              <  1:483636/145‑1 (MQ=0)
                                                                                                                            gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                          <  1:4738670/149‑1 (MQ=11)
                                                                                                                            gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTGGCGGTTTATTTGATGCCTggg                          >  2:1855572/1‑149 (MQ=2)
                                                                                                                               cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                       <  2:1206818/149‑1 (MQ=17)
                                                                                                                                 ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                     <  1:765286/149‑1 (MQ=34)
                                                                                                                                 ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                     <  1:3156228/149‑1 (MQ=34)
                                                                                                                                 ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                     <  2:2463399/149‑1 (MQ=34)
                                                                                                                                   aCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATGTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGCTAAAACGACAGGCCCAGTCTTTCGTCTGAGCCGTTCGTTTTTTTTGCTGCGTGggggattgg                   >  2:2616023/1‑149 (MQ=11)
                                                                                                                                   aCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCAACCGGCGGATTTGTCCTACACAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCACTCTGTCGACTGAGCCGGTCGTTGTATTTGATGCCTGgggaattgg                   >  2:1305550/1‑149 (MQ=11)
                                                                                                                                      cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                    >  2:2781088/1‑145 (MQ=35)
                                                                                                                                       tCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaat               <  2:721723/149‑1 (MQ=39)
                                                                                                                                        ccGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaatt              <  2:4787058/149‑1 (MQ=255)
                                                                                                                                             ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTc         >  1:1973409/1‑149 (MQ=255)
                                                                                                                                               gTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt       <  2:639613/149‑1 (MQ=255)
                                                                                                                                                 tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGa     >  2:1660794/1‑149 (MQ=255)
                                                                                                                                                 tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTGTGATTTGATGCCTGGGGACTTGGGACTGAATTCTTGa     >  1:918846/1‑149 (MQ=14)
                                                                                                                                                    |                                                                                                                                                    
GCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGA  >  CP009273/360009‑360305

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: