Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 314,258 G→A 6.8% intergenic (‑453/‑99) ykgP ← / → eaeH pseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913314,2580GA 6.8% 267.6 / 6.4 85intergenic (‑453/‑99)ykgP/eaeHpseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology;factor; Not classified
Reads supporting (aligned to +/- strand):  ref base G (44/35);  new base A (3/3);  total (47/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

TGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCA‑GCATGCGGTTCAG  >  NC_000913/314081‑314512
                                                                                                                                                                                 |                                                                                                                                                                                                                                                                 
tGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGTTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAATTCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTTTCTACCTATGTTCGTTAATTCGTTCATCGGGAAAtttccta                                                                                                                                                                  >  4:286318/1‑270 (MQ=255)
                  tataACTTAGGTATCTATTTAATTACTTGCACAAAAAGCCAAATTTTCCCCCATAAATAAAAATATAAGCCCTAGCCCAACCAACTTATGAGTGGCTATAGGCACGGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGTAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCCTCAGGAAATTATCTCAATGTCACATTATAAAAc                                                                                                                                                <  4:128132/275‑1 (MQ=255)
                         tAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACACGCACACGACATTCTTTTTGTCTGCCTATGTTCGTTAATTCGTTCTTCAGGCAATTATCCCCTTGTTCCATTATAAAAATGGTcct                                                                                                                                         >  4:120172/1‑273 (MQ=255)
                                  aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                                                                                                                                              <  7:259931/149‑1 (MQ=255)
                                  aTTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCGGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGAAAAGCAGAAGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑Aaa‑catcac                                                                                                                              >  2:369466/1‑271 (MQ=255)
                                   tttAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAg                                                                                                                                                                                                                                                             <  6:303363/149‑1 (MQ=255)
                                      aaTGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATTCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTTCCTATGTTCGTTAATTCGTTCATCCGGAAAtt                                                                                                                                                                       >  1:8709/1‑232 (MQ=255)
                                             tGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGACATCCCCGGTTTTTCTGACAAGCAAACGCCATTATTTTTGTCTGCCTATGTTCGTTAATTCGTTCACCTGGAAATTATTCCAATGGCACATTATAAAAAAGGTCTT‑AAA‑CAACCAAGTTTTCGtt                                                                                                                    >  4:24251/1‑274 (MQ=255)
                                             tGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGc                                                                                                                                                                                                                            >  4:519761/1‑172 (MQ=255)
                                                acaAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAAt                                                                                                                                                                                                                                                >  8:58835/1‑149 (MQ=255)
                                                 caAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGAc                                                                                                                                                                                                                        <  2:519766/172‑1 (MQ=255)
                                                  aaaaaGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATcc                                                                                                                                                                                                                                              >  5:207951/1‑149 (MQ=255)
                                                         tAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGttt                                                                                                                                                                                                                                       >  5:227662/1‑149 (MQ=255)
                                                         tAAACTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGATTTCGTTATTCAGTTCTGGc                                                                                                       <  2:24251/275‑1 (MQ=255)
                                                                    ccATAAATAAATATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                                                                                                                                              >  7:313751/1‑115 (MQ=255)
                                                                    ccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGaaaa                                                                                                                                                                                                                                                              <  8:313751/115‑1 (MQ=255)
                                                                    ccATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGc                                                                                                                                                                                                                            <  6:207951/149‑1 (MQ=255)
                                                                      ataaataaAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAg                                                                                                                                                                                                                          >  5:248831/1‑149 (MQ=255)
                                                                               aaTATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCatta                                                                                                                                                                                                                 <  7:58835/149‑1 (MQ=255)
                                                                                 tataATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGaa                                                                                                                                                                                                                                                          <  6:434930/106‑1 (MQ=255)
                                                                                 tataATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGaa                                                                                                                                                                                                                                                          >  5:434930/1‑106 (MQ=255)
                                                                                  ataATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTAttt                                                                                                                                                                                                              <  5:376245/149‑1 (MQ=255)
                                                                                      tCCCGCGCCCAACCGCCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGa                                                                                                                                                                                                                                                           <  8:248660/100‑1 (MQ=255)
                                                                                      tCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGa                                                                                                                                                                                                                                                           >  7:248660/1‑100 (MQ=255)
                                                                                         cgcgCCCAACCACTTGATGAGTGGCTATAGGCACTGGATATACCAGGTCGCGGCGCCCTTTCTTACATAAAGCTATTTCCTTTTCTGCGGAAACGGCAATAGTGAAATCCCCGGATTTTCTGACAAGCATACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTTATCAGGTAATGATCTCAATGTCACTTTATAAAACAGGTCTTCTAC‑CAACCCCGATTTCGTTTTTCTTTTCTCGCCGTCCGCGTTGCGTTGGCCAatctgttctt                                                                         >  4:195986/1‑266 (MQ=255)
                                                                                           cgcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGAc                                                                                                                                                                                                                        >  3:38850/1‑130 (MQ=255)
                                                                                           cgcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGATTTCGTTATTCAGTTCt                                                                                                          >  4:127730/1‑238 (MQ=255)
                                                                                           cgcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGATTTCGTTATTCAGTTCt                                                                                                          <  2:127731/238‑1 (MQ=255)
                                                                                            gcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCt                                                                                                                                                                                                    <  6:227662/149‑1 (MQ=255)
                                                                                            gcCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCt                                                                                                                                                                                                    <  6:54788/149‑1 (MQ=255)
                                                                                               cAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCa                                                                                                                                                                                                                    <  1:38850/130‑1 (MQ=255)
                                                                                                aaCCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGt                                                                                                                                                                                                >  8:215200/1‑149 (MQ=255)
                                                                                                          tgaGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGACATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTc                                                                                                                                                                                      >  5:236935/1‑149 (MQ=255)
                                                                                                             gTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACaa                                                                                                                                                                                                                              >  3:33333/1‑106 (MQ=255)
                                                                                                              tGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCg                                                                                                                                                                                     <  5:71249/146‑1 (MQ=255)
                                                                                                              tGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCg                                                                                                                                                                                     >  6:71249/1‑146 (MQ=255)
                                                                                                              tGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGa                                                                                                                                                                                                                                                           >  5:45814/1‑76 (MQ=255)
                                                                                                              tGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCCGCAGAAAAGGa                                                                                                                                                                                                                                                           <  6:45814/76‑1 (MQ=39)
                                                                                                                 cTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAg                                                                                                                                                                                                                          <  1:33333/106‑1 (MQ=255)
                                                                                                                        cACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTAc                                        <  2:286322/275‑1 (MQ=255)
                                                                                                                         aCTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAAtt                                                                                                                                                                       >  5:474385/1‑149 (MQ=255)
                                                                                                                          cTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTCTTTGTGTCTGCCTATGTTCGTTNATTCGTTCATCAGGAAATTTTCTCAATGTCACAttct                                                                                                                                                     >  3:50659/1‑164 (MQ=255)
                                                                                                                          cTGGAGATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAACCAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTAc                                        <  2:195988/274‑1 (MQ=255)
                                                                                                                           tGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                               >  4:181889/1‑155 (MQ=255)
                                                                                                                           tGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAACACAGGTCAT‑AAC‑AAACCACGATTTCGTTAGTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTTTTTTTCCACTCCCTGTCCCCTTTCcccc                                     >  4:84722/1‑275 (MQ=255)
                                                                                                                              atatatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaaaa                                                                                                                                                 <  1:50659/166‑1 (MQ=255)
                                                                                                                               tatatTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTca                                                                                                                                                           <  2:181891/155‑1 (MQ=255)
                                                                                                                                  atTAGGTGGCGGTGCACTTTCTTACATAACGCTATTTCCTTTTCTGCGGAAAAGGATATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATg                                                                                                                                                              >  8:246996/1‑149 (MQ=255)
                                                                                                                                   ttAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGt                                                                                                                                                             >  8:297774/1‑149 (MQ=255)
                                                                                                                                       gTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTcaca                                                                                                                                                         >  7:110635/1‑149 (MQ=255)
                                                                                                                                       gTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCTTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATggcgg                          >  2:491236/1‑274 (MQ=255)
                                                                                                                                        tGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCAGGTTTTTCTGCCAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTNAATGTCACATTATCAAACATGTCAT‑AAG‑AAACCACGAGTTCGTTATTAAGTTCTTGCCCGCTGAGTGGCGGGGGNCAATCGCGCTGTTCAGTTTCTTTTTCCATTNNATTTCCCCTTTCCCCCCGTAATTGCGGCa                        >  3:72970/1‑275 (MQ=255)
                                                                                                                                           cGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACa                                                                                                                                                                                                                               >  1:135197/1‑75 (MQ=255)
                                                                                                                                           cGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTc                                                                                                                                                                                              >  4:72946/1‑108 (MQ=255)
                                                                                                                                           cGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACTTTTACCCCAGTAATGGCGGCACGt                     >  2:337099/1‑275 (MQ=255)
                                                                                                                                           cGGTGCACTTTC‑TACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGc                                                                                                                                                                                                                            >  2:204076/1‑77 (MQ=255)
                                                                                                                                               gCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCa                                                                                                                                                                                                                           <  3:135197/75‑1 (MQ=255)
                                                                                                                                               gCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTa                                                                                                                                                                                          <  2:72947/108‑1 (MQ=255)
                                                                                                                                               gCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATaaaa                                                                                                                                                 <  6:236935/149‑1 (MQ=255)
                                                                                                                                               gCACTTTC‑TACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGAc                                                                                                                                                                                                                        <  4:204073/77‑1 (MQ=255)
                                                                                                                                                   tttCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAgg                                                                                                                                             <  8:110635/149‑1 (MQ=255)
                                                                                                                                                   tttCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAgg                                                                                                                                             >  6:27316/1‑149 (MQ=255)
                                                                                                                                                        tACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                               >  2:135352/1‑126 (MQ=255)
                                                                                                                                                          cATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑aaa                                                                                                                                     <  5:27316/149‑1 (MQ=255)
                                                                                                                                                            tAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTca                                                                                                                                                           <  4:135351/126‑1 (MQ=255)
                                                                                                                                                             aaaGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑Aaa‑caa                                                                                                                                 <  6:248831/149‑1 (MQ=255)
                                                                                                                                                             aaaGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑Aaa‑caa                                                                                                                                 >  6:327537/1‑149 (MQ=255)
                                                                                                                                                             aaaGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCACCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTCTAAAACAGGTCAT‑AAA‑CAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGTTTCTTTTTCCACTCGCTGTACCCTTTACCCAAGTAATGGGTGCACGTGCACACGCAT‑CGGTTCac  >  4:396521/1‑274 (MQ=255)
                                                                                                                                                              aaGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAg                                                                                                                                              <  6:135461/137‑1 (MQ=255)
                                                                                                                                                              aaGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAg                                                                                                                                              >  5:135461/1‑137 (MQ=255)
                                                                                                                                                               aGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                               >  4:163087/1‑119 (MQ=255)
                                                                                                                                                               aGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATtata                                                                                                                                                    >  2:317113/1‑130 (MQ=255)
                                                                                                                                                                 gTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCAc                                                                                                                             <  6:474385/149‑1 (MQ=255)
                                                                                                                                                                   aTTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTca                                                                                                                                                           <  2:163089/119‑1 (MQ=255)
                                                                                                                                                                   aTTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAAc                                                                                                                                                <  4:317109/130‑1 (MQ=255)
                                                                                                                                                                    tttCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGAt                                                                                                                          >  2:79519/1‑149 (MQ=255)
                                                                                                                                                                        cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGATTTCg                                                                                                                      >  5:25146/1‑149 (MQ=255)
                                                                                                                                                                        cTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGATTTCg                                                                                                                      <  4:79518/149‑1 (MQ=255)
                                                                                                                                                                         ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAgg                                                                                                                                             >  6:5836/1‑127 (MQ=255)
                                                                                                                                                                         ttttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAgg                                                                                                                                             <  5:5836/127‑1 (MQ=255)
                                                                                                                                                                          tttCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTAtc                                                                                                                                                                    >  4:32096/1‑103 (MQ=255)
                                                                                                                                                                             cTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGATTTCGTTAtt                                                                                                                 <  6:25146/149‑1 (MQ=255)
                                                                                                                                                                              tGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCaa                                                                                                                                                                <  2:32096/103‑1 (MQ=255)
                                                                                                                                                                                cGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAAt                                                                                                                                                               >  2:311434/1‑102 (MQ=255)
                                                                                                                                                                                cGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCt                                            >  4:224846/1‑215 (MQ=255)
                                                                                                                                                                                 |                                                                                                                                                                                                                                                                 
TGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTGGATATATTAGGTGGCGGTGCACTTTCTTACATAAAGGTATTTCCTTTTCTGCGGAAAAGGAAATCGGGAAATCCCCGGTTTTTCTGACAAGCAGACGCCATTATTTGTGTCTGCCTATGTTCGTTAATTCGTTCATCAGGAAATTATCTCAATGTCACATTATAAAACAGGTCAT‑AAA‑CAACCACGATTTCGTTATTCAGTTCTGGCCCGCTGCGTGGCGTGGGCAAATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCA‑GCATGCGGTTCAG  >  NC_000913/314081‑314512

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑