Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 360,154 T→G 6.4% intergenic (‑88/+507) 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,1540TG 6.4% 574.3 / 7.4 245intergenic (‑88/+507)lacZ/lacIpseudogene, truncated/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (56/165);  new base G (11/12);  total (67/178)
Fisher's exact test for biased strand distribution p-value = 2.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

AGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGA  >  CP009273/360005‑360302
                                                                                                                                                     |                                                                                                                                                    
tgggccACTTCACCGTTTCCATGCGGCACACCACCCAGCGCCACCATCCGGTCCAGGCGCTCGTTTGATACGTTGTAAAAATTGTGTTTTCCGATCCCCAATTCGTGGCTTTTTTGGGGGGGGGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                                      <  2:1141405/145‑1 (MQ=255)
agagGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                                      >  1:1532374/1‑150 (MQ=255)
agagGCACTCCCCCGCGTGCCAGCGGCTCACCAACCCGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAAGTCCGACCCCCAATTCCTGGAAGTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                                      <  2:3021898/150‑1 (MQ=37)
 tgtcctcctcaggagagcgttcaccgacacccccccgataaaaccacaggcccagtctttcgactgagcctttcgtttaatttgatgccgggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTg                                                                                                                                                     >  2:2453286‑M1/99‑150 (MQ=255)
      tactcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCg                                                                                                                                                <  2:2221220‑M1/57‑1 (MQ=255)
       actcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCATGTCAGTTTATGTCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGc                                                                                                                                               <  2:239786‑M1/58‑1 (MQ=255)
         tcaggagagcgttcaccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                             <  1:2845256‑M1/60‑1 (MQ=255)
         tCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCaa                                                                                                                                             >  1:2144816/1‑150 (MQ=255)
         tCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCaa                                                                                                                                             >  1:4272798/1‑150 (MQ=255)
          aaacGCTTGACAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAATCTGTAAAACATGAGAATTCCGATCCCCAATTCCTCGAAGTTTATGGCGGGCGTCCTGCCAGCCACCCTCCGGGCCGTTTCATAGCaaa                                                                                                                                            >  2:4057193/4‑149 (MQ=38)
           aCCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                                                                                           >  2:4254167/1‑150 (MQ=255)
             cggttGCCAGCGGCTTACCACCCCGCGCCACCATCCAGGGCAGGAGCTCCTATTATAAGCTGTAAAACATGATAATTCCGACCCCCAATTCCTGGCTGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGtt                                                                                                                                         <  1:3862973/147‑1 (MQ=25)
              gCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGTGGGTCGCAACCTTc                                                                                                                                        >  1:3388553/1‑150 (MQ=255)
               cTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCCACCTTcc                                                                                                                                       >  1:3275158/1‑149 (MQ=255)
                  gCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                                    <  2:3983422/150‑1 (MQ=255)
                     aGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                 <  2:2654187/150‑1 (MQ=255)
                     aGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                 <  2:1695711/150‑1 (MQ=255)
                      caccgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGc                                                                                                                                <  2:2979274‑M1/73‑1 (MQ=255)
                        ggCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGGTAAAATACCCTc                                                                                                                              >  1:3483453/1‑150 (MQ=37)
                         cgacaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGGTAAAATCCGCTcc                                                                                                                             <  2:1453350‑M1/76‑1 (MQ=255)
                          cTTACCATCCAGCGCCACCATCCCGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTccc                                                                                                                            <  2:3889060/150‑1 (MQ=255)
                           acaaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCCAATCCGCACCcg                                                                                                                           <  2:2555649‑M1/78‑1 (MQ=255)
                            tACCCTCCAGCGCCACCATCCTGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcgg                                                                                                                          <  1:3629873/150‑1 (MQ=255)
                            caaacaacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTGAGCAACGTTCAaatacgcagcccg                                                                                                                          <  1:1181588‑M1/79‑7 (MQ=255)
                                caacagataaaacgaaaggccccgtctttcgactgagcctttcgtgttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGGGGCGTCGCACCGTTCAACTCCGCTCccggcggc                                                                                                                      <  2:3032383‑M1/83‑2 (MQ=255)
                                aTCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGa                                                                                                                      <  2:4134265/150‑1 (MQ=255)
                                 tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGATTTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                     <  1:4304120/150‑1 (MQ=255)
                                 tCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                     <  2:4211919/150‑1 (MQ=255)
                                 aacagataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAt                                                                                                                     <  2:1766945‑M1/84‑1 (MQ=255)
                                   cAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTTTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTTTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                                                                                                   <  1:2466043/150‑1 (MQ=255)
                                      ataaaacgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                                                                                                <  2:4261715‑M1/89‑1 (MQ=255)
                                        ccaccaTCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCt                                                                                                              >  1:3585621/1‑150 (MQ=255)
                                            cgaaaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTGAGCCACGTTCAAAACCGGTCCCGGCGGGTGTGTCCtacac                                                                                                          <  1:2717108‑M1/95‑3 (MQ=255)
                                            caTCCAGTGCAGGCGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCAACGTTAAAAAACGCGCCCCGCGGCTTTGTCCTAcac                                                                                                          >  1:2393267/1‑148 (MQ=21)
                                             tccccgTGCTGGTGCTCCGTTGATAAGCTGACAAACTTGATATTTCCGAACCCAAATTCCGGGCTTTTTTTGGGGGGGGTCCTGCCCCCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCa                                                                                                         <  1:3360574/145‑1 (MQ=17)
                                               ccAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTACCGGCCGATTTGTCCTACTCacg                                                                                                       >  1:2453286/1‑148 (MQ=255)
                                               aaggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGTGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCGGCGGATTTGTACTACTCAgg                                                                                                       <  2:3337770‑M1/98‑1 (MQ=255)
                                                 ggcccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                                                                                                 <  1:2499876‑M1/72‑1 (MQ=255)
                                                 cgTGCCGGAGCTCCGATGATACGCTGTCAAACATGAGAATTCCGATCCCCAAGTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGgag                                                                                                     <  2:3611716/149‑1 (MQ=255)
                                                   tGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCCACCTTCAAATACCCAACCCGCGGATTTGTCCGACACACgagat                                                                                                   >  1:4178292/1‑149 (MQ=21)
                                                   cccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCGGCGGCGTTGTCCTACTCAGGagag                                                                                                   <  2:451555‑M1/102‑1 (MQ=255)
                                                    ccagtctttcgactgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGGTCGCAACCATAAAAAACGCTCCCCGCCGCTTTGTGCAACACAGGCGAGc                                                                                                  <  1:4155751‑M1/103‑1 (MQ=255)
                                                     cAGGAGCTCGTATGATACTCTGTCATAAATGCTACGACCGTTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCg                                                                                                 <  2:599122/150‑1 (MQ=14)
                                                      aGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGt                                                                                                <  2:3544064/150‑1 (MQ=255)
                                                       ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                               <  2:4162673/150‑1 (MQ=255)
                                                       ggAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGtt                                                                                               <  2:1972622/150‑1 (MQ=255)
                                                         aGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCa                                                                                             <  2:3813323/150‑1 (MQ=255)
                                                          gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCTTTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                            <  1:2672244/150‑1 (MQ=255)
                                                          gCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCATTTTCTGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCAc                                                                                            <  1:1322095/150‑1 (MQ=255)
                                                                ctgagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCt                                                                                                                                                   >  2:4268958‑M1/36‑89 (MQ=255)
                                                                 tgatAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCAACCTTAAAATACGGACCCGGCGGCTTTTTGGTAAACAGGTGGGCGTTAACACacaaa                                                                                     >  2:370003/1‑150 (MQ=18)
                                                                  tatAAGCTGTAAAACAGGACAATTCCTATCCCAAATACGTGGCATTTTATGGCGGGCGTCCGCCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaac                                                                                    <  1:456309/149‑1 (MQ=14)
                                                                  gagcctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAACTTCAAAACCGCGCCCCGCGGAATTGTCCAACTCAGGAGAGCGTTCACCCacaaac                                                                                    <  1:2862227‑M1/117‑1 (MQ=255)
                                                                   agcctttcgttttatttgatgcctggggaagtAATTCCTGGCCGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTGCGCCACGTTCAAATCCGCTCCCGGCGGCTGTGTCCTACTCAGGGGTGCggtcgccggccatca                                                                                   <  2:2504090‑M1/118‑8 (MQ=255)
                                                                    tAAGCTCACAAACATGAGAATACCGACCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGacaaacaa                                                                                  <  1:1951313/150‑1 (MQ=255)
                                                                     cctttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAAGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGCAAaacaac                                                                                 <  2:1570209‑M1/120‑1 (MQ=255)
                                                                     aaGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGGGTCGCAACCTTCAAATACGCGCCCCGCGGATTTGTGCAACAAAGGAGAGCGGTAAACCacaaacaaa                                                                                 >  1:3521524/1‑149 (MQ=17)
                                                                        ttcgttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAACCGCTCCCGGCGGATTTGTCCTACACAGGAGTGCGTTCACCGACAAACCACAGa                                                                              <  1:3995563‑M1/123‑1 (MQ=255)
                                                                           tCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                           <  2:3915038/150‑1 (MQ=255)
                                                                           gttttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATaa                                                                           <  2:2472886‑M1/126‑1 (MQ=255)
                                                                            cAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAACGGTGAAAAACCCTCCCCCCCGGGTTGTAGGACGAAGGAGTGCGTTGACCCACAAAAAACAaaaaaa                                                                          >  1:906227/1‑150 (MQ=12)
                                                                              ttatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCATCGCAACCTTCAAATAACCTCCCCCCGGATTTGTCCAAATCATGAGAGCGTTCACCGACCAACAACAGATAAAAc                                                                        <  1:2323946‑M1/129‑1 (MQ=255)
                                                                              aaCATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTAGCAACCTTCAAAAACCCACCCGGCGGATTTGTGCTACGAAGGTGAGCGGTAACCCGCAAACAACAGCTAAAAc                                                                        >  1:2254706/1‑150 (MQ=17)
                                                                               tatttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                       >  1:4260219‑M1/21‑150 (MQ=255)
                                                                               aCATGAGAATACCGATCCCCAATTCCTGGCAGTTTAGGGCGGGCGTCCGGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACg                                                                       <  1:3846884/150‑1 (MQ=255)
                                                                                cATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGCGCGTTCACCGACAAACAACACATAAAACGa                                                                      >  2:1523453/1‑150 (MQ=255)
                                                                                atttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGTATTGTCCTACTCAGGGGAGCGTTGACCGACACACAACAGTTAAAACGa                                                                      <  2:3736905‑M1/131‑1 (MQ=255)
                                                                                 tttgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTGGCAACGTTCAAATCCGCTCCCGGGGGGTTTGTCCGACTCAGGAGGGGGGTGCCCGCCAAACAACAGATAAAacgta                                                                     <  1:1675401‑M1/132‑3 (MQ=255)
                                                                                   tgatgcctgggtaattAATTCGTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAg                                                                   >  2:3657642‑M1/17‑150 (MQ=255)
                                                                                   tgatgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTCGCAACCTTCAAAAACGCTCCCGGCGGATTTGTCCTACACAGGAGAGCGGTCAACCCCCAACCACAGATAAaaagaaag                                                                   <  2:764703‑M1/134‑1 (MQ=255)
                                                                                     atgcctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATACGCTCCCCGCGGCTTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAAAGAAAGGc                                                                 <  2:3372534‑M1/136‑1 (MQ=255)
                                                                                         ctggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAg                                                             <  2:311023‑M1/140‑1 (MQ=255)
                                                                                          tggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                                                                                                       <  2:3506708‑M1/66‑1 (MQ=255)
                                                                                          tggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAACGTTCAAATCCGCACCCGGCGGTTTTGTCCCTCTCCGGAGAGCGTTCACCCCCAAACAAAAGATAAATCCAATGGCCCAGt                                                            <  1:969356‑M1/141‑1 (MQ=255)
                                                                                           ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTAGCCACCTTCAAAACCGCGCCCCGCGGATTTGTCCTAAACAGGGGAGCGGTCACCGACAAACAACACATAAAAAGAaaggccccctc                                                           <  1:3109259‑M1/142‑5 (MQ=255)
                                                                                           ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                            >  1:1051097‑M1/9‑133 (MQ=255)
                                                                                           ggggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATa                                                                            <  2:1051458‑M1/125‑1 (MQ=255)
                                                                                            gttaattAATTCCTGGCATTTTAGGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                          >  2:846358‑M1/8‑150 (MQ=255)
                                                                                            ggtaattAATTCGTGGCTGTTTTTGGGGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGTTTTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                          >  1:3337770‑M1/8‑150 (MQ=255)
                                                                                            gggaattAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCt                                                          >  1:1766945‑M1/8‑150 (MQ=255)
                                                                                             aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                         <  2:2155885/150‑1 (MQ=38)
                                                                                             aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCtt                                                         <  2:2512589/150‑1 (MQ=38)
                                                                                             aTCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCTTCGCCACGTTCAAATCCGCTCCCGGCGGCTTTGTGCTAATCAGGGGGGGGGTCCCCGACCAACAACAGATTAAACGGAAGGCCCCGTCtt                                                         >  2:3679632/1‑150 (MQ=11)
                                                                                               ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                       <  1:628891/150‑1 (MQ=35)
                                                                                               ccccAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTc                                                       <  2:2715504/150‑1 (MQ=35)
                                                                                                  tAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTGGCGTCGCAACGGTCAAATCCGCGCCCCGCGGCTTTGGCCTACACAGGAGAGCGTTCACCCACCAACAACAGATAAAACGAAAGGCCCAGTCATTCGAc                                                    <  2:418177‑M1/149‑1 (MQ=255)
                                                                                                   atttCCTGGCTTTTTCTGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                   <  1:843478/148‑1 (MQ=9)
                                                                                                   aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACt                                                   <  2:1819434/150‑1 (MQ=32)
                                                                                                   aaTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGTGGCAAAGGTCAAAACCGCGCCCGGCGGGTTTGTCCTACACAGGAGGGCGTTCACCGCCCAACCACAGATAAAACGCAAGGGCCAGTCTTTCGTCt                                                   >  1:1999055/1‑150 (MQ=11)
                                                                                                    aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                  >  2:532174/1‑150 (MQ=32)
                                                                                                    aTTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTg                                                  >  2:2412019/1‑150 (MQ=32)
                                                                                                                                gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTgggg                      <  2:2932360/150‑1 (MQ=31)
                                                                                                                                  ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                    <  1:2718167/150‑1 (MQ=32)
                                                                                                                                  ccGCCACCCTCCGGGCCGTAGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAACACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaa                    <  2:2909998/150‑1 (MQ=17)
                                                                                                                                   cGCCACCCGCCGGGCCGTTGCATCGCAACGTACAACTCCGCGCCCGGCGGATTGGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGATAACACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaat                   <  2:1903363/150‑1 (MQ=11)
                                                                                                                                    gCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAACTCCGCTCCCGGCGGACTGGTCCGACTCAGGAGAGCGTTCACCGACAAACAACAGATACAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                  <  2:2094991/150‑1 (MQ=11)
                                                                                                                                    gCCACCCGCCGGGCCGTTGCTTCGCAACGTTCACAGCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACACCAGATAACACGAAAGGCCCAGGCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaatt                  <  2:4216966/150‑1 (MQ=11)
                                                                                                                                        cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattggga              <  1:2979274/150‑1 (MQ=38)
                                                                                                                                                    ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATGGGGAATTAATTCTTGa  >  1:2884406/1‑150 (MQ=255)
                                                                                                                                                     tGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGgggaattg                 >  1:3932847/1‑134 (MQ=35)
                                                                                                                                                     |                                                                                                                                                    
AGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCCGATGATAAGCTGTCAAACATGAGAATTCCGATCCCCAATTCCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCAGGAGAGCGTTCACCGACAAACAACAGATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCTTTCGTTTTATTTGATGCCTGGGGAATTGGGAATTAATTCTTGA  >  CP009273/360005‑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: