Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,157 T→G 8.1% 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.1% 655.0 / 9.8 279intergenic (‑91/+504)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (72/172);  new base G (19/13);  total (92/187)
Fisher's exact test for biased strand distribution p-value = 1.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 282 total aligned reads displayed.

GGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGA  >  CP009273/360008‑360302
                                                                                                                                                      |                                                                                                                                                 
ggCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCt                                                                                                                                                   >  1:810094/1‑150 (MQ=255)
  ctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttgtatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTc                                                                                                                                                 <  2:1399808‑M1/56‑1 (MQ=255)
   aCTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCg                                                                                                                                                >  2:278836/1‑150 (MQ=255)
    cTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGGTTCgg                                                                                                                                               >  1:4485752/1‑149 (MQ=255)
     ttCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCTGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                              <  2:1299711/150‑1 (MQ=255)
      tcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGAAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                             <  2:1127441‑M1/60‑1 (MQ=255)
       caggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                                                                            <  2:679144‑M1/61‑1 (MQ=255)
       caggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                                                                            <  2:4805177‑M1/61‑1 (MQ=255)
       cACCGCTTGCCAGCGGCTAACCAACCCGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTAAAACATTATAATTCCGAGCCCAACTTCCTGGCTTTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                                                                            <  1:4756966/150‑1 (MQ=25)
        aCCGCTTGCCCGGGGCTCCCCAGCCCGCCCCCCCAGCCTGTGCCGGAGCGCCTATGATACGCTGACAAACTTCAGACTCCCGAGCCCCATGTCCTGGCTGTTTTGGGGGGGGGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                           <  1:4733827/150‑1 (MQ=12)
          gagagcgttcaccgacaaacaacagataacacgaacggcccagtctttcgactgcgcctttcgttttattggcggccgggggaattAATGCCGGGCTGTTTAGGGCGGGCGGCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                         >  1:1888632‑M1/87‑150 (MQ=255)
             agcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaa                                                                                                                                      >  2:3140565‑M1/84‑150 (MQ=255)
               gCCAGCGGCTTACCATCCAGCGCCACCATCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATGCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                    <  2:3881344/150‑1 (MQ=255)
               gCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                    >  1:3424895/1‑150 (MQ=255)
                 cAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATcc                                                                                                                                  <  2:516205/150‑1 (MQ=255)
                    accgacaaccaacagagaaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                               >  2:1505435‑M1/77‑150 (MQ=255)
                    accgacaaacaccagataaaacgaacggccctgtctttcgaatgagcatttctttttatttgatgcatgggtaattAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCt                                                                                                                               >  1:3636928‑M1/77‑150 (MQ=255)
                      cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTcc                                                                                                                             <  2:3132888‑M1/76‑1 (MQ=255)
                       attACCATCCAGCGCAACAACCCAGGGCTGGAGCTCGGTTGATAAGCTGTCAAACATGTGAATGCCGAACCCAATTTCCTGGCATTTTTTGGCGGGCGTCGTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                            <  1:3132018/149‑1 (MQ=18)
                          aaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttattggatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                         <  2:587548‑M1/80‑1 (MQ=255)
                          aCCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggc                                                                                                                         >  2:2690011/1‑150 (MQ=255)
                           ccATCCAGCGCCACCTGCCAGTGCAGGAGCTCGGTTGATAAGCTGTCAAACATTTGAATTCCGATCCCAATTTCCTGGAAGTTTATGGCGGGCGTCCTGCCGGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcg                                                                                                                        <  2:4084527/150‑1 (MQ=21)
                           ccATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGCTTCGCAACGTTCAAATCCCCCCCcggcg                                                                                                                        >  2:4552647/1‑150 (MQ=255)
                           aacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggacttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCACCGGTCCAGTCCGCGCCcggcg                                                                                                                        <  2:4337781‑M1/81‑1 (MQ=255)
                            cATCCCGCGCCACCATCCTGTGCAGGAGCTCCTATGATATGCTGTCAAACATGATAATTCCGATCCCAAATTCCTGGCTTTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                                                                                                       <  2:1609256/150‑1 (MQ=21)
                            cATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTGGCCACCGTCAAATACCCGCCcggcgg                                                                                                                       >  1:2528018/1‑150 (MQ=37)
                                cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                   <  2:3195400/150‑1 (MQ=255)
                                  gcgcCACCATCCAGTCCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGt                                                                                                                 <  2:3365072/150‑1 (MQ=255)
                                      caccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAAACTTAAAAAACGCACCCCGCGGATTTGttcaa                                                                                                             >  1:1425771/1‑146 (MQ=25)
                                       accaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTAAAATACGCTACCGGCGGATTTGTCCTaa                                                                                                            >  2:3025782/1‑149 (MQ=255)
                                       aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                            <  2:4709746‑M1/93‑1 (MQ=255)
                                       aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                            >  1:2294506‑M1/58‑150 (MQ=255)
                                          aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCTATCCCCAATTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                         <  1:4413347/150‑1 (MQ=255)
                                          aTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                         <  1:1962779/150‑1 (MQ=255)
                                           tCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATGCGTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                        <  2:16174/150‑1 (MQ=255)
                                              aGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCATGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                     >  2:2656127/1‑150 (MQ=255)
                                                tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTTCGTAGCAACGGTCAAATACCCTACCCGCGGATTTTTGCTACACAGGagag                                                                                                   >  1:1062941/1‑150 (MQ=21)
                                                 ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                  <  2:1150188‑M1/103‑1 (MQ=255)
                                                   aGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                <  2:394922/150‑1 (MQ=255)
                                                    gtctttcgcctgagcttttttttttatttgttgcctgggtaattATTTCCGGGCTGTTTTTGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                               >  1:550460‑M1/45‑150 (MQ=255)
                                                    gtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAAGTCAAATACGCTACCCGCCGCTTTGTCCAAAACACGAGAGCGtt                                                                                               <  1:2743732‑M1/106‑1 (MQ=255)
                                                       tttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATACGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                            <  2:1884422‑M1/109‑1 (MQ=255)
                                                       gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                            <  1:1559056/150‑1 (MQ=255)
                                                        ctcgattgcgcttgtcgatttatttgaggcctgggtaattAATTCCGGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                           >  1:2977297‑M1/41‑150 (MQ=255)
                                                        cTCCGATGATATGGTTTCAAACATGAGAATTCCGACCCCCAAGTCCTGGCTTTTTTGGGCGGGCGTCCTCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGCTCAcc                                                                                           <  1:2116077/150‑1 (MQ=17)
                                                        cTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                           <  2:3047389/150‑1 (MQ=255)
                                                            gatgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaca                                                                                       >  2:2694202/1‑150 (MQ=255)
                                                             atgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                      <  1:3246520/150‑1 (MQ=255)
                                                                agcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                   <  2:289581‑M1/118‑1 (MQ=255)
                                                                 gcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGaccaacaa                                                                                  <  2:2593597‑M1/119‑1 (MQ=255)
                                                                    tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTGCAAATACGCTACCGGCGGCTTTTTCCAACTCAGGAGAGCGGTGACCGACAAACAACAg                                                                               <  1:2404236‑M1/122‑1 (MQ=255)
                                                                    gCTGTAAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAg                                                                               <  2:2074127/150‑1 (MQ=255)
                                                                     cTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGa                                                                              <  2:4277290/150‑1 (MQ=255)
                                                                      tGTCAAACATGAGATTTCCGATCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                             <  1:1466467/150‑1 (MQ=255)
                                                                      tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTTTGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                             <  1:1466350/150‑1 (MQ=255)
                                                                      tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCCACCTTAAAATCCGCTCCCGGCGGATTTGTCCGACGCAGGAGAGCGTTCACCGCCAAACAACAGaa                                                                             >  1:754184/1‑149 (MQ=17)
                                                                        gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCACTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                           >  2:2728159‑M1/25‑150 (MQ=255)
                                                                        gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATCCGCACCCGGCGGATTTGTCCAACTCAGGAGAGCGTTCACCGACAAACAAcagaaaa                                                                           <  2:2325804‑M1/126‑4 (MQ=255)
                                                                          aaaattGATAATGCCGATCCCCAATGCGTGGCTTTTTAGGGGGGGGGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                         <  1:902623/145‑1 (MQ=11)
                                                                          aaaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                         <  2:4241469/150‑1 (MQ=255)
                                                                           ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTAAAATACGCTCCCCGCGGAGTTGTGCTACAAAGGAGTGCGTTCACCCACAAACAACAcaaaaaaa                                                                        <  1:650219‑M1/129‑2 (MQ=255)
                                                                           ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCAAAGTTAAAATACCCACCCCGCGGATTTGTGCTACACAGGAGAGCGGTCACCGACAAACAAACGAtaaaaa                                                                        <  1:3984247‑M1/129‑2 (MQ=255)
                                                                            tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGCGCGTTCACCGACAAACAACAGATAAAACg                                                                       <  1:1117120‑M1/130‑1 (MQ=255)
                                                                            acttGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                       <  2:3300546/147‑1 (MQ=255)
                                                                            aCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                       <  2:3563818/150‑1 (MQ=255)
                                                                             cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGa                                                                      >  2:3473016/1‑150 (MQ=255)
                                                                             cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTGCCAAAGGTCAAATACGCACCCCGCGGATTTTTGCTAATGAGGAGATCGTTGACCGAAAAAAAACAGAAAAAAAGa                                                                      >  1:1024547/1‑150 (MQ=17)
                                                                                   tgcgtgggtaattAATTCGTGGCAGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGcc                                                                >  2:3251507‑M1/14‑150 (MQ=255)
                                                                                     ttCCGATCCCCAATTCCGGGCAGTTTTTGGGGGGCGTCAGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                              <  1:2779369/150‑1 (MQ=14)
                                                                                     cctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAAACTTAAAAAACGCTCCCCGCGGATTTGTCCTAAACAGGAGAGCGTTCACCGACCAAAAACAGAAAAAAAGAAAGGCCCa                                                              <  1:1141603‑M1/139‑1 (MQ=255)
                                                                                       ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                            <  2:902001/150‑1 (MQ=40)
                                                                                       ccGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGt                                                            <  2:1233568/150‑1 (MQ=40)
                                                                                         gcgccccAAGTCCGGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                          <  2:3657406/147‑1 (MQ=16)
                                                                                           tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTAAAATACCCTACCCGCGGATTTTTCCTACTCAGGAGTGCGTTAACAGACAAACAAAAGATAAAAAGAAACGACCAGTCttt                                                        >  1:1540992/1‑150 (MQ=11)
                                                                                            ccccAATGCCTGGATGTTTCTGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                       <  1:924144/150‑1 (MQ=11)
                                                                                             ctcAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                      <  1:3951504/148‑1 (MQ=18)
                                                                                             attaaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTAGCCACGTTCAAAACCGCAACCGGCGGAATTGGGCGACTCCGGAGAGCGGGGACCCACAAACAACAGATAAAACGAAAGGCCCCCTCTTTCg                                                      >  1:2687621/4‑150 (MQ=11)
                                                                                              ccAATTCCCGGGCAGTTTTTGGCGGGCGCCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                      <  1:1800981/150‑1 (MQ=11)
                                                                                               tacttCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCGACCCTCCGGGCCGTTGCGTCGCTAGGTGCAAATCCGCTGCCGGCGGTTTTGTCCTACTCTGGAGAGCGTTCACGGACAACCAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                    <  2:253501/147‑1 (MQ=11)
                                                                                               cAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                    <  2:1866234/150‑1 (MQ=33)
                                                                                                aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCAAATCCGCTCCCGGCGGATTTGTCCTACGCAGGAGAGCGTTGACCGACCAACAACAGATAAAACGCAAGGCCCAGTCTTTAGACt                                                   >  2:3736663/1‑150 (MQ=11)
                                                                                                 aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAc                                                                                                                                            >  2:402078/1‑60 (MQ=35)
                                                                                                  ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCAAAGGTCAAAAACGCTCCCGGGGGGTTTGTGCGGAACAGGCGCGCGGTCACCCCCCAACAACCGAAAAAAAGAAAGGCCCCCTCTTTGGACTgg                                                 >  2:2018372/1‑149 (MQ=2)
                                                                                                   tCCGGGCAGTTTCTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAg                                                <  2:4496041/150‑1 (MQ=11)
                                                                                                            gtttcgggcgggcgTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGt                                        <  1:2917341/144‑1 (MQ=0)
                                                                                                                           cccgCCCCCCACCCTCCGGGCCGTCGCCTCGCAACGGACAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTg                         <  2:1263996/147‑1 (MQ=1)
                                                                                                                             tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCGTTCGGTTTATTTGATGCCTggg                       >  2:3895945/1‑150 (MQ=11)
                                                                                                                                ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                    <  2:1540992/150‑1 (MQ=32)
                                                                                                                                 ccccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                   <  1:2866430/148‑1 (MQ=18)
                                                                                                                                 cGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCGGgggaat                   >  2:360981/1‑150 (MQ=18)
                                                                                                                                  gCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                  <  2:3523518/150‑1 (MQ=21)
                                                                                                                                      catcccGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCTGCGGATTTGTCCTACTCAGGAGTGCGTTCACCGACAAAAAACAGATAAAACGAAAGGCCCAGTCTTTCTACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggga              <  1:4175548/146‑1 (MQ=11)
                                                                                                                                            gggCCGTTGCTTCGCAACGTTCAAATCCCCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTACACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaat        <  1:4047950/150‑1 (MQ=255)
                                                                                                                                             ggCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaatt       <  1:2713642/150‑1 (MQ=255)
                                                                                                                                               ccGCTGCTTCGCAACGCTCAAATCCGCTCCCCGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGCTAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt     <  2:3984247/150‑1 (MQ=255)
                                                                                                                                                cGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGACTTGGGCCTTAATTCtt    >  2:1797515/1‑150 (MQ=255)
                                                                                                                                                  ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGACTTAATTCTTGa  >  2:4765679/1‑150 (MQ=255)
                                                                                                                                                  ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGa  >  2:4777559/1‑150 (MQ=255)
                                                                                                                                                     cTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                   >  1:4846406/1‑130 (MQ=34)
                                                                                                                                                     cTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                   <  2:4849460/130‑1 (MQ=34)
                                                                                                                                                      |                                                                                                                                                 
GGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGA  >  CP009273/360008‑360302

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: