Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 152,527 T→A 7.4% T96S (ACT→TCT)  yadM ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913152,5270TA 7.4% 186.3 / 5.3 54T96S (ACT→TCT) yadMputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base A (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

CTTATTGATAGATGAAATTAAACGTCGCCAGTGCACGAAAATTCCCCGGTACGCCTTGCCCTGCATCTGTTTCACGTAGCGCAACGGTCATTTCAAGACCCTTATCGGGCTGCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGT  >  NC_000913/152242‑152780
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                             
cttATTGATAGATGAAATTAAACGTCGCCAGTGCACGAAAATTCCCCGGTACGCCTTGCCCTGCATCTGTTTCACGTAGCGCAACGGTCATTTCAAGACCCTTATCGGGCTGCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGt                                                                                                                                                                                                                                                            >  1:379944/1‑289 (MQ=255)
cttATTGATAGATGAAATTAAACGTCGCCAGTGCACGAAAATTCCCCGGTACGCCTTGCCCTGCATCTGTTTCACGTAGCGCAACGGTCATTTCAAGACCCTTATCGGGCTGCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCACTTCCGGTCAGAGAGGt                                                                                                                                                                                                                                                            >  1:529929/1‑289 (MQ=255)
                 ttAAACGTCGCCAGTGCACGAAAATTCCCCGGTACGCCTTGCCCTGCATCTGTTTCACGTAGCGCAACGGTCATTTCAAGACCCTTATCGGGCTGCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGc                                                                                                                                                                                                                                           <  2:552926/289‑1 (MQ=255)
                                        aTTCCCCGGTACGCCTTGCCCTGCATCTGTTTCACGTAGCGCAACGGTCATTTCAAGACCCTTATCGGGCTGCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAAcc                                                                                                                                                                                                                    <  1:279853/289‑1 (MQ=255)
                                                               cATCTGTTTCGCGTAGCGCAACGGTCATTTCAAGACCCTTATCGGGCTGCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGAt                                                                                                                                                                                               >  2:693090/1‑287 (MQ=255)
                                                               cATCTGTTTCACGTAGCGCAACGGTCATTTCAAGACCCTTATCGGGCTGCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAACGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTTGTCGAAGATGAATCACCAGAATGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCAGGTCAGAGTGGTATCAATCCAACTGATGCAACTGCTGCACCCACTGGCAACCAGTTTAAACTCATCCTgagct                                                                                                                                                                                               >  2:529174/1‑285 (MQ=255)
                                                                         aCGTAGCGCAACGGTCATTTCAAGACCCTTATCGGGCTGCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTa                                                                                                                                                                                   >  1:509936/1‑289 (MQ=255)
                                                                                            tCAAGACCCTTATCGGGCTGCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGAc                                                                                                                                                                <  1:693090/289‑1 (MQ=255)
                                                                                                     ttATCGGGCTGCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAg                                                                                                                                                       >  1:27919/1‑289 (MQ=255)
                                                                                                     ttATCGGGCTGCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAg                                                                                                                                                       >  1:40057/1‑289 (MQ=255)
                                                                                                               gCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCaa                                                                                                                                             >  1:540523/1‑289 (MQ=255)
                                                                                                                     cGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTa                                                                                                                                       <  2:27919/289‑1 (MQ=255)
                                                                                                                              tCCAGCGTATCTTTTCCGAACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTtattat                                                                                                                               <  1:519355/288‑1 (MQ=255)
                                                                                                                                cAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGttatta                                                                                                                                >  2:324343/1‑285 (MQ=255)
                                                                                                                                cAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGttatta                                                                                                                                <  1:324343/285‑1 (MQ=255)
                                                                                                                                     tATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCcatca                                                                                                                       >  1:530200/1‑289 (MQ=255)
                                                                                                                                            tCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCa                                                                                                                                              >  2:135771/1‑259 (MQ=255)
                                                                                                                                            tCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCa                                                                                                                                              <  1:135771/259‑1 (MQ=255)
                                                                                                                                                               ggAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTg                                                                                                                 >  1:712455/1‑269 (MQ=255)
                                                                                                                                                               ggAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTg                                                                                                                 <  2:712455/269‑1 (MQ=255)
                                                                                                                                                               ggAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGttatt                                                                                                         >  2:59009/1‑277 (MQ=255)
                                                                                                                                                               ggAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGttatt                                                                                                         <  1:59009/277‑1 (MQ=255)
                                                                                                                                                               ggAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAgtgt                                                                                             >  1:477867/1‑289 (MQ=255)
                                                                                                                                                                 aaaGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCACTTCCGGTCAGAGAGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGtgtg                                                                                            <  2:529929/288‑1 (MQ=255)
                                                                                                                                                                   aGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGttatt                                                                                                         >  1:49572/1‑273 (MQ=255)
                                                                                                                                                                   aGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGttatt                                                                                                         <  2:49572/273‑1 (MQ=255)
                                                                                                                                                                   aGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATg                                                                                         >  1:11754/1‑289 (MQ=255)
                                                                                                                                                                                  agtCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCt                                                                          <  2:540523/289‑1 (MQ=255)
                                                                                                                                                                                    tCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCttt                                                                        >  1:340990/1‑289 (MQ=255)
                                                                                                                                                                                      cGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTTACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTg                                                                       >  2:135263/1‑288 (MQ=255)
                                                                                                                                                                                        ttttttGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGtta                                                                                                           >  2:310670/1‑250 (MQ=255)
                                                                                                                                                                                        ttttttGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGtta                                                                                                           <  1:310670/250‑1 (MQ=255)
                                                                                                                                                                                         tttttGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAg                                                                   <  1:456236/289‑1 (MQ=255)
                                                                                                                                                                                                   ccATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTaaa                                                         <  2:509936/289‑1 (MQ=255)
                                                                                                                                                                                                   ccATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTaaa                                                         <  2:477867/289‑1 (MQ=255)
                                                                                                                                                                                                                   tGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATcc                                         <  2:379944/289‑1 (MQ=255)
                                                                                                                                                                                                                    gtcgtcGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCCTCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTTGCGCTTTGACAGTGGCTGTAAACTATACATTAATATCCg                                        >  2:229882/1‑289 (MQ=255)
                                                                                                                                                                                                                      cgtcgAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATc                                                                                                                                                                          >  1:356208/1‑157 (MQ=255)
                                                                                                                                                                                                                      cgtcgAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATc                                                                                                                                                                          <  2:356208/157‑1 (MQ=255)
                                                                                                                                                                                                                         cgAAGATTAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                   <  1:229882/289‑1 (MQ=255)
                                                                                                                                                                                                                                   tCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAgg                         <  2:11754/289‑1 (MQ=255)
                                                                                                                                                                                                                                   tCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAgg                         <  2:40057/289‑1 (MQ=255)
                                                                                                                                                                                                                                       cAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGa                     <  2:285653/289‑1 (MQ=255)
                                                                                                                                                                                                                                                gTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataata             <  2:530200/288‑1 (MQ=255)
                                                                                                                                                                                                                                                     ataaGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCa       <  1:529174/289‑1 (MQ=255)
                                                                                                                                                                                                                                                          cTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt  >  1:28324/1‑289 (MQ=255)
                                                                                                                                                                                                                                                             aGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCAt      <  1:279509/282‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCAAt      >  2:279509/1‑282 (MQ=255)
                                                                                                                                                                                                                                                             aGGTGAGCTTGATGATGCATTTCCGGTCAGAGAGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATc                                                                                                                                                                          <  1:580423/118‑1 (MQ=255)
                                                                                                                                                                                                                                                             aGGTGAGCTTGATGATGCATTTCCGGTCAGAGAGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATc                                                                                                                                                                          >  2:580423/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                                      gTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATc                                                                                                                                                                          <  1:632886/93‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      gTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATc                                                                                                                                                                          >  2:632886/1‑93 (MQ=255)
                                                                                                                                                                                                                                                                                           agTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCt                                                                                                                                                          >  2:54111/1‑104 (MQ=255)
                                                                                                                                                                                                                                                                                           agTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCt                                                                                                                                                          <  1:54111/104‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                             
CTTATTGATAGATGAAATTAAACGTCGCCAGTGCACGAAAATTCCCCGGTACGCCTTGCCCTGCATCTGTTTCACGTAGCGCAACGGTCATTTCAAGACCCTTATCGGGCTGCATCTCGTCTGTGCTCCAGCGTATCTTTTCCGCACTGTTAGGTTTAAGGAAAGTGGCATCATCAGTAGTCCGTTTTTTGAAACCCATACCGATATTACTTGTCGTCGAAGATGAATCACCAGACTGCGGTATAATAAGCTTAGGTGAGCTTGATGATGCATTTCCGGTCAGAGTGGTATCAATCCAACTGATGCCACTGCTGCACCCACTGGCAACCAGTTTAAAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGT  >  NC_000913/152242‑152780

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