Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,157 T→G 8.6% 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.6% 431.1 / 5.6 188intergenic (‑91/+504)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (40/124);  new base G (10/13);  total (51/137)
Fisher's exact test for biased strand distribution p-value = 7.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 193 total aligned reads displayed.

CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGA  >  CP009273/360010‑360305
                                                                                                                                                   |                                                                                                                                                    
ctactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCtt                                                                                                                                                     <  2:134591‑M1/55‑1 (MQ=255)
     cACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCaa                                                                                                                                                >  1:1936796/1‑149 (MQ=255)
     cACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                                <  2:1672719/149‑1 (MQ=255)
      aggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTGCGCAAc                                                                                                                                               <  1:2270254‑M1/61‑1 (MQ=255)
      aCCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAAc                                                                                                                                               >  1:2727904/1‑149 (MQ=255)
       ccGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCGGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCATGCCCGCCACCCTCCGGGCCGGGGGTTCGCagcc                                                                                                                                              >  1:1839055/1‑146 (MQ=38)
         gCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATACGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                            <  2:257075/149‑1 (MQ=255)
         agagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                            <  2:4221431‑M1/64‑1 (MQ=255)
           ttGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATACGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                                          <  2:3589252/149‑1 (MQ=255)
              gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCGTGGCAGTTTAAGGCGGGCGTCCGGCCCGCCCCCCGCCGGGGCGTTGGTGCGCGACGTTCGAAt                                                                                                                                       <  2:930567‑M1/69‑1 (MQ=255)
              gttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAAt                                                                                                                                       <  1:1041727‑M1/69‑1 (MQ=255)
              cccgCGGCTTACCATCCAGCGCCACCCTCCCGTGCAGGAGCTCCGATGATACGCTGTCAAACAGGAGAATTCCGAACCCCACTTCCTGGCTGTTTTTGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCaact                                                                                                                                       <  1:1771873/146‑3 (MQ=25)
              ccAGCGGCTTACCATCCAGCGCCACCATCCCGTGCAGGAGCTCCTATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                       <  2:1122329/149‑1 (MQ=255)
                 gCGGCTTACCATCCAGCGCCACCATCCAGTGCCGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                    >  1:3849675/1‑149 (MQ=255)
                  cGGCTCACCATCCCGCGCCACCAGCCAGTGCTGGCGACCCGTTGATACGCTGTAAAACATCATAAGTCCGACCCCAATTTCCTGGCTGTTTATGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                                   <  1:3797355/149‑1 (MQ=17)
                  accgacaaacaacagatacaacgaaaggcccagtctttcgcatgagcctgtcgttttatttgatgcctggggaatcAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                                   >  2:3865899‑M1/77‑149 (MQ=255)
                   ccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgccgggggaattATTTCCTGGCAGGTTATGGCGGGCGTCCTGCGCGCCACCCTCCGGGCCGTGGGGTGGCACCGTTAAAATCCgcg                                                                                                                                  <  2:1496339‑M1/74‑2 (MQ=255)
                      ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                               <  1:3964964/149‑1 (MQ=255)
                      ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                               >  2:873772/1‑149 (MQ=255)
                      ttACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTccc                                                                                                                               >  2:708575/1‑149 (MQ=255)
                      acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                               >  2:1325085‑M1/73‑149 (MQ=255)
                       caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCCACGTTCAAAACCGCTCCcg                                                                                                                              <  2:2539442‑M1/78‑1 (MQ=255)
                            aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCACCCCGCGGa                                                                                                                         <  1:2470500‑M1/83‑1 (MQ=255)
                             ccAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                        <  2:918094/149‑1 (MQ=255)
                             acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACCTTCAAATACGCTcccggcgggt                                                                                                                        <  1:651916‑M1/84‑3 (MQ=255)
                             acagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGgag                                                                                                                        <  2:3228446‑M1/84‑2 (MQ=255)
                                tataaaacgacagccccagtctttcgcctgagcctttcgttttatttgatgccggggtaattAATTCCTGGCAGTTTAGGGCGGGGGTCCTGCCCGCCACCCTCCGCGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                     >  1:2059880‑M1/63‑149 (MQ=255)
                                gcgcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCCCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                     >  2:1275968/1‑149 (MQ=255)
                                gataacccgaacggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTg                                                                                                                     >  2:3125557‑M1/63‑149 (MQ=255)
                                  gcCACCATCCAGTGCAGGAGCTCCGATTATAAGCTTTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                   <  1:320171/149‑1 (MQ=255)
                                  gcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGATAATTCCGACCCCCAAGTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                   <  1:2406318/149‑1 (MQ=255)
                                  gcCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGAGTTGta                                                                                                                   >  2:3344147/1‑148 (MQ=255)
                                   aaaacgaaaggcccagtctttcgactgagcctttcgttgtatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTACTGCCCGCCACCCTCCGGGCCGTTGGTTAGCCACGTTGAAATACGCTAccggcggaagtttac                                                                                                                  <  1:2044293‑M1/90‑8 (MQ=255)
                                     aacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTACCGGCGGAGTTGTCCTa                                                                                                                <  2:321232‑M1/92‑1 (MQ=255)
                                          ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAg                                                                                                           <  2:1994886/149‑1 (MQ=255)
                                            ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGga                                                                                                         <  2:2790113‑M1/99‑1 (MQ=255)
                                                cAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGc                                                                                                     <  2:2863054/149‑1 (MQ=255)
                                                  ggAGCTCGGATGATATGCTGTCAAACATGATAATTCCGATCCCCAATTCGTGGCAGTTTTTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                   <  1:1678907/149‑1 (MQ=18)
                                                   gAGCTCCGATGATATGCTGTCAAACATGAGAATTCCGATCCCAAATTCCTGGCTGTTTATGGGGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                                  <  1:3060433/149‑1 (MQ=18)
                                                    ctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                 <  2:2115056‑M1/107‑1 (MQ=255)
                                                    ctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTc                                                                                                 <  2:2115819‑M1/107‑1 (MQ=255)
                                                       tcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTAGCAAACTTCAAAAACGCACCCCGCCGATTTGTCCAACACAGgagagcggtgaac                                                                                              <  1:4179033‑M1/110‑7 (MQ=255)
                                                       tCCGCTGATAAGCTGTCACACATGAGAAGTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                              <  2:998227/149‑1 (MQ=255)
                                                       tCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAcc                                                                                              <  2:1756187/149‑1 (MQ=255)
                                                         cGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGGTCACCGa                                                                                            >  2:3595410/1‑149 (MQ=255)
                                                           atgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAAACCCACCCCGCGGCATTGTACTAAACAGGGGTGCGTTAACCGaaa                                                                                          >  1:1605051/1‑147 (MQ=18)
                                                            tgagcctttcgttttatttgatgcctggggaattAATTCATGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaa                                                                                         <  2:1773922‑M1/115‑1 (MQ=255)
                                                              atAAGCTGTCAAACATGTTATGTCCGAGCCCCAATTCCTGGCAGTTTTGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                       <  2:815417/149‑1 (MQ=14)
                                                               tatgCTGTCAAACATGAGATTTCCGATCCCCTTTTCCTGGCTGTTTTGGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCATCGacaaaca                                                                                      <  1:163368/146‑1 (MQ=17)
                                                               tAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaaca                                                                                      <  2:4150710/149‑1 (MQ=255)
                                                                  tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAacaaca                                                                                   <  2:3804047‑M1/121‑1 (MQ=255)
                                                                  tttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCAACGTTAAAATCCGCACCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACcgaccaacaaaa                                                                                   <  2:1138657‑M1/121‑3 (MQ=255)
                                                                    tGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCCACCTTCAAATACCCACCCGGCGGGTTTGTACTACGCAGGAGAGCGGTCACCCACCAacaacaca                                                                                 >  2:2527781/1‑147 (MQ=17)
                                                                     gTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGAt                                                                                <  2:3529584/149‑1 (MQ=255)
                                                                      tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                               <  2:3451664/149‑1 (MQ=255)
                                                                      catagttttgatgcccggtgatttAATTCCGGGCTGTGTATGGCGGGCGTCCGCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGatc                                                                               >  1:992957‑M1/25‑148 (MQ=255)
                                                                       cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                              <  2:3371646/149‑1 (MQ=255)
                                                                       cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                              >  2:342918/1‑149 (MQ=255)
                                                                         aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaaaa                                                                            <  1:276446/149‑1 (MQ=255)
                                                                          aCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTCCAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCCAACCACCAGATAAAAc                                                                           >  2:3014651/1‑149 (MQ=17)
                                                                             tGAGAATACCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGAGTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGaa                                                                        <  2:3956465/149‑1 (MQ=255)
                                                                               gatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCGCCCCGCGGCTTTTTGCTACTACGGAGAGCGTTCACCCACAAACAACACATAAAACGAAAg                                                                      <  2:1590508‑M1/134‑1 (MQ=255)
                                                                                gaATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTGGCAACGTTCAAAAACGCGCCCGGCGGAGTTGTGCTAAGCAGGGGAGCGGTCACCGACCAACAACAGATATAACGCAAgg                                                                     >  1:4212518/1‑149 (MQ=14)
                                                                                 tgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACCTTCCAATCCCCGCCCGGCGGATTTGTCCTACTCAGGAGAGCGGTCACCCACAAACAACAGAAAAAaccaaacgc                                                                    <  2:1248486‑M1/136‑4 (MQ=255)
                                                                                   ttCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACCAACAACAGATAAAACGAAAGGccc                                                                  <  2:2153542/149‑1 (MQ=255)
                                                                                   tcccGATCAACAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACGCAGGAGAGCGTTCACCGACAAACCACAGATAAAACGAAAGGccc                                                                  <  2:1616878/147‑1 (MQ=255)
                                                                                   gtCCGATCACAAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGccc                                                                  <  2:1936796/148‑1 (MQ=255)
                                                                                    ctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCa                                                                 <  2:322532‑M1/139‑1 (MQ=255)
                                                                                     tggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                                <  2:1664825‑M1/140‑1 (MQ=255)
                                                                                     tggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAAACTTCAAATACCCACCCGGCGGATTTGTCCTAAACACGAGCGCGTTCACCCACAAACAACACAAAAAAAGAAAGGCCCAg                                                                <  2:1933245‑M1/140‑1 (MQ=255)
                                                                                     ccgtgcaccAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                                <  2:1652869/142‑1 (MQ=12)
                                                                                       gATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTc                                                              >  2:2693081/1‑149 (MQ=38)
                                                                                        ggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAACTTCAAATCCCCTCCCGGCGGCTTTGTCCTACACAGGCGAGCGTTCACCGCCAAACAACAGAAAAAACGAAAGGCCCAGTCt                                                             <  2:1208690‑M1/143‑1 (MQ=255)
                                                                                         tCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                            <  2:1168469/149‑1 (MQ=37)
                                                                                           cccATTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                          <  2:1882193/149‑1 (MQ=17)
                                                                                            ttAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCg                                                         <  2:3224896‑M1/147‑1 (MQ=255)
                                                                                              aaTTCCTGGATGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGAc                                                       <  1:942665/149‑1 (MQ=17)
                                                                                               aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTTTCCTACTCAGGAGAGCGTTAACCGCCAAACAACAGATAAAAAGAAATGCCCAGTCTTTCGACt                                                      >  2:537936/1‑149 (MQ=11)
                                                                                                ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                            >  1:2862652/1‑62 (MQ=33)
                                                                                                ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                            <  2:2862652/62‑1 (MQ=33)
                                                                                                ttCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                     <  2:2727966/149‑1 (MQ=31)
                                                                                                ttCCGGGCTGTTTTGGGCGGGCGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                     <  1:4050746/149‑1 (MQ=11)
                                                                                                                     cgcccTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATGTGTCCCACTCACGAGAGCGTTCACCGACAAACAACAGATACACCGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTAATTGATg                                <  2:2905856/146‑1 (MQ=1)
                                                                                                                     cgaccTGCCCGCCACCCTCCGGGCCGTTGCTGCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATg                                <  2:2620105/146‑1 (MQ=1)
                                                                                                                        ccggCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCt                             <  1:472003/146‑1 (MQ=1)
                                                                                                                            cccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                         <  1:3806454/149‑1 (MQ=31)
                                                                                                                                ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAACGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAACCGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                     <  1:1299132/149‑1 (MQ=255)
                                                                                                                                ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCCGGAGAGCGTTCACCGACAAACAACAGATAAAACGACAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                     <  1:1275968/149‑1 (MQ=21)
                                                                                                                                ccACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                     <  1:3224896/149‑1 (MQ=34)
                                                                                                                                 cACCCTCCGGGCCGTTGCTTCGCAACGGTAAAAAACGCTCCCCGGGGATTTGTCCTACACACGAGAGCGTTAACCGACAAACAAAAAATAAAACGAAATGCCCAGTAATTCGACTGTGCCTTTCGTTTTATTTTTTTCCTGgggaaatt                    >  2:1795884/1‑146 (MQ=11)
                                                                                                                                    ccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggga                 >  2:1043612/1‑149 (MQ=38)
                                                                                                                                     cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaa                <  1:1051962/149‑1 (MQ=38)
                                                                                                                                      tCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaat               <  2:2428601/149‑1 (MQ=39)
                                                                                                                                       ccGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattgggaatt              <  1:2109832/149‑1 (MQ=40)
                                                                                                                                         gggCCGTTGCTTCTCAACGTTCAAATCCGCTCCCGGCTGATTTTTCCTACTCATGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaa            <  1:2900042/149‑1 (MQ=255)
                                                                                                                                           gCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGaattaatt          <  1:1664825/149‑1 (MQ=255)
                                                                                                                                            ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTc         >  2:662511/1‑149 (MQ=255)
                                                                                                                                             cGTTGCTTCGCAACGTACAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTACACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt        <  1:13575/149‑1 (MQ=255)
                                                                                                                                             cGTTGCTACGCAACGTACCAATCCGCTCCCGGCGGATTTGTCCTACCCAGGAGAGCGTTCCCCGACACACCACAGATAACACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCt        <  1:1496339/149‑1 (MQ=14)
                                                                                                                                                   ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGa  <  1:1334339/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGAAGA  >  CP009273/360010‑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: