Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,053,561 A→T 6.1% intergenic (‑233/‑82) yeeL ← / → shiA pseudogene, glycosyltransferase homology/shikimate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,053,5610AT 6.1% 222.0 / 3.4 66intergenic (‑233/‑82)yeeL/shiApseudogene, glycosyltransferase homology/shikimate transporter
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base T (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

ACTAAATGCAATACGTCTCAACAATGATGCTAAAAACATACCTAACCTCGCCTCCCTACTGGTTATAATGCAATGCAGTCTATCAGACTCATCAGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGCAAGTAAGCCCGGCGATGGGAACGCTCGC  >  NC_000913/2053290‑2053836
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                   
aCTAAATGCAATACGTCTCAAAAAAGATGGTAAAAACAGACCTAACCTCGCCCCCCTGCGTGTTATAATGCAATGCAGTCTATCAGACTCATCAGTGTGCCATTTTGTGCATATGCGGACTGTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATaa                                                                                                                                                                                                                                                                                  <  4:326887/275‑1 (MQ=255)
          aTACGTCTCAACAATGATGCTAAAAACATACCTAACCTCGCCTCCCTACTGGTTATAATGCAATGCAGTCTATCAGACTCATCAGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTta                                                                                                                                                                                                                                                                        >  2:520137/1‑275 (MQ=255)
                ctcAACAATGATGCTAAAAACATACCTAACCTCGACTCCCTACTGGTTATAATGCAATGCAGTCTATCAGACTCATCAGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCCCTAACCTGTGGGTCTTCTGCTTATTACTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCCTGTTGCCGATAATCAGCCATCAACTCACTGGCCGATGCGCTCTGTCCAGGTACTCTCATAAGGTCCGTTAATACATc                                                                                                                                                                                                                                                                  >  4:543390/1‑275 (MQ=255)
                                                      ccTACTGGTTATAATGCAATGCAGTCTATCAGACTCATCAAGGTGCCATTTTGTGCATATGCGTACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTACTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATaa                                                                                                                                                                                                                            <  2:543391/275‑1 (MQ=255)
                                                               tataATGCAATGCAGTCTATCAGACTCATCAGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCNTCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTCCAGATTCGCTCTGTCGCGTTTGTCACATAAGGTTCTTNGATAAATAGACGCCAATGTGAGTTANCTTCTACTTACTAATTATaaaaaaaaaataaaa                                                                                                                                                                                                               >  3:130296/1‑269 (MQ=255)
                                                                 taATGCAATGCAGTCTATCAGACTCATCAGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTa                                                                                                                                                                                                                  <  2:449222/274‑1 (MQ=255)
                                                                         tgcaGTTTATCGAACTCATCATGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGa                                                                                                                                                                                                         <  4:520136/275‑1 (MQ=255)
                                                                         tgcaGTCTATCAGACTCATCAGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAAATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGa                                                                                                                                                                                                         >  4:74020/1‑275 (MQ=255)
                                                                         tgcaGTCTATCAGACTCATCAGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATCTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAAtt                                                                                                                                                                                                                                >  4:38869/1‑252 (MQ=255)
                                                                             gTCTATCAGACTCATCAGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATCTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATaa                                                                                                                                                                                                                            <  2:38869/252‑1 (MQ=255)
                                                                              tCTATCAGACTCATCAGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCAt                                                                                                                                                                                                     >  2:640704/1‑274 (MQ=255)
                                                                                  tCAGACTCATCAGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCAt                                                                                                                                                                                                                                                                                    <  2:18137/191‑1 (MQ=255)
                                                                                          atcaGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCTTAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACg                                                                                                                                                                                        >  2:104666/1‑275 (MQ=255)
                                                                                          atcaGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAAtg                                                                                                                                                                                                                                                        >  4:551995/1‑211 (MQ=255)
                                                                                              gggTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAg                                                                                                                                                                                                                                                    <  2:551996/211‑1 (MQ=255)
                                                                                                    cATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTa                                                                                                                                                                                                                  >  4:574448/1‑239 (MQ=255)
                                                                                                        ttGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACaaa                                                                                                                                                                                                              <  2:574449/239‑1 (MQ=255)
                                                                                                           tgCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGttt                                                                                                                                                                                                                                                                         <  3:185285/177‑1 (MQ=255)
                                                                                                           tgCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGttt                                                                                                                                                                                                                                                                         >  1:185285/1‑177 (MQ=255)
                                                                                                                                    tatCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATaa                                                                                                                                                                                                                            >  2:46054/1‑197 (MQ=255)
                                                                                                                                        tctAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAg                                                                                                                                                                                                                        <  4:46054/197‑1 (MQ=255)
                                                                                                                                         ctAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATCCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAACCAGGCTCATTATGTAATCCACGCTCATCTCCACTCGTCTCGATGACTGGACGCTTTCGTTAAGTCgcgc                                                                                                                                         >  2:214960/1‑275 (MQ=255)
                                                                                                                                                 tgGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGCATTTTTCATGTTGCCGCTAATCAGCCATAAACTCATTTGCAGATGCGCTCTGTCGAGGNAGTCTCANAAGNTTCGTTNATAGATCGACGGCAATGATAGTTACTTTTTCCATACTAATTATAAAAAGACAGTACaca                                                                                                                                                                                                              >  3:29060/1‑196 (MQ=255)
                                                                                                                                                 tgGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCAt                                                                                                                                                                                                                                       >  4:705830/1‑173 (MQ=255)
                                                                                                                                                 tgGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAGAGACAGTACAAACAGGATCAt                                                                                                                                                                                                     >  2:390295/1‑207 (MQ=255)
                                                                                                                                                 tgGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACaaa                                                                                                                                                                                                              <  1:29060/198‑1 (MQ=255)
                                                                                                                                                     tCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATaa                                                                                                                                                                                                                                                                                  <  4:719954/126‑1 (MQ=255)
                                                                                                                                                     tCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTa                                                                                                                                                                                                                                   <  2:705831/173‑1 (MQ=255)
                                                                                                                                                     tCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAGAGACAGTACAAACAGGATCATTATg                                                                                                                                                                                                 <  4:390294/207‑1 (MQ=255)
                                                                                                                                                        tctGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGtt                                                                                                                          >  2:22286/1‑275 (MQ=255)
                                                                                                                                                               aaTCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAg                                                                                                                    >  2:455946/1‑274 (MQ=255)
                                                                                                                                                                         caacaCCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTg                                                                                                         <  2:74020/275‑1 (MQ=255)
                                                                                                                                                                            cacCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCTTAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATaa                                                                                                                                                                                                                            >  2:715056/1‑157 (MQ=255)
                                                                                                                                                                            cacCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGNCAATGTGAGTTACCTTTTCCATANTAATNATaaa                                                                                                                                                                                                                           >  4:668018/1‑158 (MQ=255)
                                                                                                                                                                            cacCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATAc                                                                                                                                                                                                                                     >  2:409882/1‑148 (MQ=255)
                                                                                                                                                                            cacCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTgg                                                                                                           >  4:584605/1‑270 (MQ=255)
                                                                                                                                                                            cacCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCcgtcgt                                                                                                  >  1:169033/1‑279 (MQ=255)
                                                                                                                                                                            cacCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCcg                                                                                                      >  2:268195/1‑275 (MQ=255)
                                                                                                                                                                             acCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATaaa                                                                                                                                                                                                                           >  2:661438/1‑157 (MQ=255)
                                                                                                                                                                                aGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCTTAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAg                                                                                                                                                                                                                        <  4:715055/157‑1 (MQ=255)
                                                                                                                                                                                aGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAAt                                                                                                                                                                                                                                 <  4:409881/148‑1 (MQ=255)
                                                                                                                                                                                aGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGa                                                                                                                                                                                                                       <  2:668019/158‑1 (MQ=255)
                                                                                                                                                                                aGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGcc                                                                                                       <  2:584606/270‑1 (MQ=255)
                                                                                                                                                                                 gCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGa                                                                                                                                                                                                                       <  4:661437/157‑1 (MQ=255)
                                                                                                                                                                                      tCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAAt                                                                                                                                                                                                                                 <  3:5875/142‑1 (MQ=255)
                                                                                                                                                                                      tCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAAt                                                                                                                                                                                                                                 >  1:5875/1‑142 (MQ=255)
                                                                                                                                                                                         tGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAg                                                                                                                                                                                                                                                    >  4:578809/1‑120 (MQ=255)
                                                                                                                                                                                             cTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTAc                                                                                                                                                                                                                                                <  2:578810/120‑1 (MQ=255)
                                                                                                                                                                                                tCATCTTCCATCGAATTTTTCATGTTGCCGTTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACATTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCTTTAAGTCGCCCCCGCCGCGCTCCGTTGGGCAGCTTCCCTTGTTCCGTCGCCGACTGGTATGtttt                                                                                   >  4:134280/1‑274 (MQ=255)
                                                                                                                                                                                                tCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTCATCCCATAAGGTTCGTTTATAGATCGACGGCAATGTCAGTTACCATTTCCATACTAATTATAAAAAGACAGTCCAAACAGCATCCTAATGCACTCCCCGCTCATCTCTCCCACTCCCCATGAAGGGCCGCTTTCGTTAAGTCGCTCCCGACGCGCTGCGTTCGGCCGCCGCGCTGGTGCCGTTGTCGCCCGGttgtatttt                                                                                  >  4:649113/1‑267 (MQ=255)
                                                                                                                                                                                                tCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATaaa                                                                                                                                                                                                                           >  2:126483/1‑138 (MQ=255)
                                                                                                                                                                                                    cTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGa                                                                                                                                                                                                                       <  4:126483/138‑1 (MQ=255)
                                                                                                                                                                                                          tCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCCTAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTAt                                                                         >  2:118011/1‑274 (MQ=255)
                                                                                                                                                                                                                 tttCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATtat                                                                                                                                                                                                                              <  3:110919/118‑1 (MQ=255)
                                                                                                                                                                                                                 tttCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATtat                                                                                                                                                                                                                              >  1:110919/1‑118 (MQ=255)
                                                                                                                                                                                                                  ttCATGTTGCAGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCAccgccg                                                            <  1:130296/279‑1 (MQ=255)
                                                                                                                                                                                                                         tgcgctaaaccaaccATAAAATCATTTAACGGTGCGCATTGATGGGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCGTACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCAc                                                         <  4:22286/262‑1 (MQ=255)
                                                                                                                                                                                                                         tGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGGGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCAc                                                         <  2:134280/275‑1 (MQ=255)
                                                                                                                                                                                                                         tGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGGTCGTTGATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGTGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCAc                                                         <  2:467501/275‑1 (MQ=255)
                                                                                                                                                                                                                         tGCCGCTAATCAGCAATAAAAGACTTTGCAGATGCCGCCTGTCGAGGTAGTTTCATAAGGTTCGGTTAAAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGAAc                                                         <  4:455945/275‑1 (MQ=255)
                                                                                                                                                                                                                                                     cagtgtcgcTCTGTTGGGGTAGTCTCTTAAGGTTCGTTTATAGATCGACGGCAATTGGAGTTACCGTTTCCATACTAATTATAAAAAGACAGTACAACCAGGATCCTTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCg                               <  4:104666/267‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGCAAGTAAGCCCGGc                 <  4:118011/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tgttgctcAAAAGGTTGGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAATACAGTACAAACAGGATCATTATGGACTCCACGCCTATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGCAAGTAAGCCCGGCGATggg           <  4:458849/270‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCAACGCTCATCTCCACTCGTCCCGATGAAGGGCCGCTTTCGTTAAGTCGCGCCCGAAGAGCAGCGTTAGGCAGCTGCGCTGGTGCCGTCGGCTACTGGCATGCGTTTTTACGCTATGGAAACACCGCCGCCCTGGGGTTTAATTTCGAGNNNTTCCAGCAGGTAAACCCGGCGAGGGGAACGcc     >  3:217921/1‑277 (MQ=255)
                                                                                                                                                                                                                                                                            ctcATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTa                                                                                                                                                                                                                  >  2:154636/1‑71 (MQ=255)
                                                                                                                                                                                                                                                                            ctcATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGCAAGTAAGCCCGGCGATGGGAACGCTcgc  >  1:169224/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                   
ACTAAATGCAATACGTCTCAACAATGATGCTAAAAACATACCTAACCTCGCCTCCCTACTGGTTATAATGCAATGCAGTCTATCAGACTCATCAGGGTGCCATTTTGTGCATATGCGGACTTTTATGTTTCATATCTCTAACCTGTGGGTCCTCTGCTTAATCCTTAAACAACACCAGCAACTCCTGCGCTTTCATCTTCCATCGAATTTTTCATGTTGCCGCTAATCAGCCATAAAATCATTTGCAGATGCGCTCTGTCGAGGTAGTCTCATAAGGTTCGTTTATAGATCGACGGCAATGTGAGTTACCTTTTCCATACTAATTATAAAAAGACAGTACAAACAGGATCATTATGGACTCCACGCTCATCTCCACTCGTCCCGATGAAGGGACGCTTTCGTTAAGTCGCGCCCGACGAGCTGCGTTAGGCAGCTTCGCTGGTGCCGTCGTCGACTGGTATGATTTTTTACTCTATGGCATCACCGCCGCACTGGTGTTTAATCGCGAGTTTTTCCCGCAAGTAAGCCCGGCGATGGGAACGCTCGC  >  NC_000913/2053290‑2053836

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: