Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,157 T→G 7.2% 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 7.2% 538.2 / 2.8 217intergenic (‑91/+504)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (60/135);  new base G (13/9);  total (73/144)
Fisher's exact test for biased strand distribution p-value = 1.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 226 total aligned reads displayed.

GCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAA  >  CP009273/360009‑360303
                                                                                                                                                    |                                                                                                                                                  
 ctactcaggagtgcgttcaccgacaaacaacagatacaccgaaaggccctgtgtttctactgaccctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                   >  2:4148552‑M1/95‑149 (MQ=255)
 cACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATACGCTGTCAACCATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                   <  2:2052502/149‑1 (MQ=255)
 cACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTAAAACATGAGAATTCCGATCCCCACTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                   <  1:1417936/149‑1 (MQ=255)
    ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGc                                                                                                                                                >  1:1695901/1‑149 (MQ=255)
    ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGc                                                                                                                                                >  2:2758972/1‑149 (MQ=255)
    ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGc                                                                                                                                                >  1:3844638/1‑149 (MQ=255)
     tcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                               <  2:2835273‑M1/59‑1 (MQ=255)
     tcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                               <  2:635492‑M1/59‑1 (MQ=255)
      cACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                              <  2:3233341/149‑1 (MQ=255)
        ggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTAGCAACg                                                                                                                                            <  1:3508728‑M1/62‑1 (MQ=255)
          agagcgttccccgacacacaacagctaaaccgacaggcccagtctttcgcctgcgcctgtcgttttatttgatgcctggggaattATTTCCGGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                          >  2:863417‑M1/86‑149 (MQ=255)
           cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTTCGCCACGTTc                                                                                                                                         >  2:13248/1‑149 (MQ=255)
            ttGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGTTAGCAAccataa                                                                                                                                        >  1:1480688/1‑144 (MQ=255)
             gcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                       <  2:273198‑M1/67‑1 (MQ=255)
               gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTATATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCCACGGTGAAAt                                                                                                                                     <  1:1215769‑M1/69‑1 (MQ=255)
                cAGCGGCTTACCAACCCGCGCCACCATCCCGTGCAGGAGCTCCTTTTATATGGTGAAAAACATGCTAATTCCGACCCCAAATCCCTGGCATTTTTGGGGGGGCGTCCGGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                                                                                                    <  1:1296124/149‑1 (MQ=16)
                  caccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                  <  2:4375161‑M1/72‑1 (MQ=255)
                    ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                                >  2:1785335/1‑149 (MQ=255)
                      gacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATACGCTcc                                                                                                                              <  1:1734280‑M1/76‑1 (MQ=255)
                       ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                             <  2:3199621/149‑1 (MQ=255)
                             tCCAGCGCCACCATCCAGGGCAGGAGCTCCGATGATACGCTGTCAAACATGAGAAGTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                       <  1:3886140/149‑1 (MQ=255)
                             tCCAGCGCCACCAGCCAGGGCCGGAGCTCCTATGATAAGCTGTCACCCATGAGATTTCCGCGCCCCAATGCCTGGCTGTTTATGGGGGGCGTCGGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                       <  2:2969132/149‑1 (MQ=18)
                              acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGGTCAAATACGCTcccggcggct                                                                                                                      <  2:2669947‑M1/84‑3 (MQ=255)
                                 gcgcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAACGTTCAAAAACGCTCCCGGCGGCTTTg                                                                                                                   >  1:4271106/1‑149 (MQ=25)
                                   taaaacgaatggcccagtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGCTGCGCTACGTGCCATTCCGCTCCcggcggagtggtc                                                                                                                 <  2:3227992‑M1/89‑7 (MQ=255)
                                   gcCACCATCCCGTGCAGGAGCTCCGATGATCAGCTGTCAAACTTGAGAATTCCGATCCACAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                 <  2:1028321/149‑1 (MQ=21)
                                   gcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATACGCTCCCGGCGGCTTTGTc                                                                                                                 >  2:3692810/1‑149 (MQ=255)
                                    ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                                                                                                <  2:1770437/149‑1 (MQ=255)
                                    ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGGTGAAATCCGCTCCCGGCGGATTTGtgc                                                                                                                >  1:2100904/1‑147 (MQ=255)
                                    aaaacgaaaggcccagtctttcgactgagcctgtcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCCCACCCGGCGGAGTTGTcc                                                                                                                <  1:3261061‑M1/90‑1 (MQ=255)
                                       ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAACCCTGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                             <  2:3026828/149‑1 (MQ=255)
                                       ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                             <  2:612541/149‑1 (MQ=255)
                                       ccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAACCTTCAAATACGCTCCCCGCGGAGTTGTCCTAc                                                                                                             >  1:1382203/1‑149 (MQ=25)
                                       acgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCCACGGTCAAATACGCGCCCGGCGGctttgtgctga                                                                                                             <  1:3183388‑M1/93‑6 (MQ=255)
                                         gaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAACGTTCAAAAACGCTCCCCGCGGCTTGGACCTACTc                                                                                                           <  2:2210540‑M1/95‑1 (MQ=255)
                                         aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCCATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                           <  2:436872/149‑1 (MQ=255)
                                          aaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTAGTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCAACCGGGCCGTGGCTTAGCAACGTTAAATTGCGCTCCCGGGGGAGTTGTCctacgca                                                                                                          <  2:253133‑M1/96‑4 (MQ=255)
                                               tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCAACGATCAAATAACCTCCCGGCGGATATGTAATACTCAGgaga                                                                                                     >  1:257990/1‑149 (MQ=21)
                                               tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                     <  2:2384372/149‑1 (MQ=255)
                                               cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCGCCCGGCGGATTTGTCCTACTCAGgaga                                                                                                     <  2:1017812‑M1/101‑1 (MQ=255)
                                                gccggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTGCGCAACCTTCAAAAACGCACCCCGCGGATTTGTCCGACGCAGgaggg                                                                                                    >  1:3866592/4‑147 (MQ=21)
                                                gCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGAGCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGagag                                                                                                    <  2:2097912/149‑1 (MQ=255)
                                                  agtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                  <  1:634229‑M1/104‑1 (MQ=255)
                                                   ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCATTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                 <  1:306066/149‑1 (MQ=255)
                                                   ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACCGTCAAATCCCCGCCCGGCGGATTTGTCCTACTCAGGAGAGccg                                                                                                 >  2:1268364/1‑147 (MQ=21)
                                                    gAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                >  2:128899/1‑149 (MQ=255)
                                                     aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                               <  1:1415476/149‑1 (MQ=255)
                                                      tttcgactgagcctttcgttttatgtgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTTTCCTACTCAGGAGAGCGTTCa                                                                                              <  1:2483490‑M1/108‑1 (MQ=255)
                                                       ttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                             >  1:1821947‑M1/41‑149 (MQ=255)
                                                       cTCCGATGTTAAGCTGTCAAACATGAGAATTCCGAGCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                             <  2:2612404/149‑1 (MQ=255)
                                                          gactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGa                                                                                          <  1:4315509‑M1/112‑1 (MQ=255)
                                                           actgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGac                                                                                         <  2:3883477‑M1/113‑1 (MQ=255)
                                                             tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                       <  1:681318/149‑1 (MQ=255)
                                                             tgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                       >  2:1723844‑M1/35‑149 (MQ=255)
                                                               agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGacaaac                                                                                     <  2:295492‑M1/117‑1 (MQ=255)
                                                                    ttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                                <  1:4007224‑M1/122‑1 (MQ=255)
                                                                     tcgttttatttgatgcctggggaattAATTCCGGGCAGTTTATGGCGGGCGTCCTGACCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGATCGTTCACCGACAAACAACAGa                                                                               >  2:818787‑M1/27‑149 (MQ=255)
                                                                      gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                              <  2:3241591/149‑1 (MQ=255)
                                                                        cAAACTTGAGATTTCCGATCCCCATTTCCTGGCTGTTTTGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                            <  1:2138529/149‑1 (MQ=14)
                                                                         aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTGGCAACCTTCAAATACGCGCCCGCCGGAGTTGTACAAATCACGAGAGCGTTCACCCACAaaaaacccaaaaa                                                                           >  1:3896391/1‑142 (MQ=17)
                                                                          aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTAGCAAACTTAAAATACCCACCCGGCGGTGTTTTGCGAAAACAGTGAGCGTTAACCCACAAAAAACAGATaaaa                                                                          >  2:752226/1‑149 (MQ=12)
                                                                            cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAAAGTTCAAAAAAGCGCCCGGCCGGTTTTTGCTACACAGGAGAGGGTTCACCGACAAACAACAGATaaaaag                                                                        >  1:1004222/1‑147 (MQ=17)
                                                                            atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAATCCGCTCCCGGCCGAGTTGGCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                        <  2:324589‑M1/130‑1 (MQ=255)
                                                                            atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCGTCAAATCCCCTCCCCGCCGATTTGTCCGACACACGAGAGCGTTCACCCACAAACAACACAAAAAacc                                                                        <  1:2912872‑M1/130‑2 (MQ=255)
                                                                             tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGa                                                                       <  2:1950090‑M1/131‑1 (MQ=255)
                                                                                gatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                             <  2:1373184‑M1/109‑1 (MQ=255)
                                                                                 gaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGGGTGGCCACGGTGAAAAAACCTACCCGCGGGTTTGTGCGAATGAGGCGAGCGTTGACCGAGCAAAAACACAAAAAACGgaatg                                                                   >  1:3739242/1‑147 (MQ=9)
                                                                                 atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAgg                                                                   <  2:325217‑M1/135‑1 (MQ=255)
                                                                                 atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCGGCGGGTTTGTCCTACGCAGGGGCGCGGTCACCGACAACCAACAGATAAAACGAAAgg                                                                   <  2:326174‑M1/135‑1 (MQ=255)
                                                                                  tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCGCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGc                                                                  >  1:3069253‑M1/14‑149 (MQ=255)
                                                                                  aaTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGc                                                                  <  2:3843094/149‑1 (MQ=255)
                                                                                   aTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                 <  2:2786968/149‑1 (MQ=255)
                                                                                    gatgtgtaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                >  2:1310773‑M1/12‑149 (MQ=255)
                                                                                     tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                               <  2:4117811/149‑1 (MQ=40)
                                                                                     tCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                               <  1:1660775/149‑1 (MQ=40)
                                                                                        gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACGTTCAAAACCGCGACCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAAACAACAGATAAAAAGAaaggcccccta                                                            <  1:3252538‑M1/142‑5 (MQ=255)
                                                                                          tCCCCAATTCCTGGCTTTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                          <  1:2777244/149‑1 (MQ=255)
                                                                                          gaactAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                          >  2:3160060‑M1/6‑149 (MQ=255)
                                                                                           aattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTGCACCGCCACACAACAGATAAAACGAAAGGCCCCGTCttt                                                         <  2:3953815‑M1/145‑1 (MQ=255)
                                                                                            cccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                        <  2:343007/149‑1 (MQ=35)
                                                                                            attaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAAACGTCAAATCCGCACCCCGCGGATTTGTACTACAACGGAGAGCGGTCACCGCCCAACCACACAAAAAAAGAAAGGCCCCCTCttta                                                        >  1:3230458/4‑148 (MQ=2)
                                                                                            attAATCCCTGGCAGTTTTGGGCGGGGGGCCGGCCCGCCACCCCCCGGGCCGTTGCTTCGCAACGTTCACATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                        >  1:326174‑M1/4‑149 (MQ=255)
                                                                                             ccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGGTCAAAAACGCTCCCGGCGGGTTTGTCCAAAAAAGGAGAGCGTTGAACGACAAACAACAGATAAAAAGAAAGGCCCAGTCTTTCg                                                       >  1:305100/1‑149 (MQ=11)
                                                                                              tgaggcgtGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGa                                                      <  1:2706809/142‑1 (MQ=11)
                                                                                               aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTACTACTCAGGAGAGCGTTAACCCACAAACAACAAAAAAAACGAAAGGCCCAGTCTTTAGAc                                                     >  1:2849131/1‑149 (MQ=11)
                                                                                                aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                    <  2:1182195/149‑1 (MQ=32)
                                                                                                aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                    >  2:2878698/1‑149 (MQ=32)
                                                                                                  tCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                  <  1:4350363/149‑1 (MQ=11)
                                                                                                  tCCTGGCAGTGTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGa                                                  <  2:3074470/149‑1 (MQ=11)
                                                                                                       gctgTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCtt                                             <  1:201700/146‑1 (MQ=1)
                                                                                                               gggggggcgTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAACCAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTa                                     <  1:1268364/145‑1 (MQ=0)
                                                                                                                            gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTCCTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTggg                        <  2:1480688/149‑1 (MQ=11)
                                                                                                                              ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGggga                      <  2:2452853/149‑1 (MQ=32)
                                                                                                                                gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                    <  2:908274/149‑1 (MQ=33)
                                                                                                                                 ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                   <  2:1895108/149‑1 (MQ=34)
                                                                                                                                  cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                  <  2:540351/149‑1 (MQ=35)
                                                                                                                                          gggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCGGGGGACGTGGGaattaa          >  2:1038082/1‑149 (MQ=255)
                                                                                                                                           ggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaat         <  1:1407629/149‑1 (MQ=255)
                                                                                                                                             ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACACACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTGGATTTGATGCCGGGGGCAGTGGGCATTAATTc       >  1:836115/1‑149 (MQ=11)
                                                                                                                                               ggtGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTACACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCtt     <  1:2291214/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAA  >  CP009273/360009‑360303

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: