Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,794,829	A→G	T281T (ACA→ACG)	rfaD →	ADP‑L‑glycero‑D‑mannoheptose 6‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,794,829	0	A	G	100.0%	199.9 / NA	57	T281T (ACA→ACG)	rfaD	ADP‑L‑glycero‑D‑mannoheptose 6‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/28);  total (29/28)

AGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCTGATGTGAATCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCGGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATG  >  NC_000913/3794582‑3795072                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    aGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGCCTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAAt                                                                                                                                                                                                                                                  >  1:145751/1‑251 (MQ=255)  gCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATc                                                                                                                                                                                                                                                 <  1:415910/251‑1 (MQ=255)     gTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTgc                                                                                                                                                                                                                                              <  1:313833/251‑1 (MQ=255)         gAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTgcgcg                                                                                                                                                                                                                                           <  1:14614/250‑1 (MQ=255)               aGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCggg                                                                                                                                                                                                                                    <  1:261329/251‑1 (MQ=255)                          cAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACaaa                                                                                                                                                                                                                          <  1:138722/250‑1 (MQ=255)                          cAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAAc                                                                                                                                                                                                                         <  1:311053/251‑1 (MQ=255)                            aaCGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCg                                                                                                                                                                                                                       >  1:349455/1‑251 (MQ=255)                                   cTTCGTCTAAGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGATCaaa                                                                                                                                                                                                                 >  1:97619/1‑250 (MQ=255)                                     tCGTCTATGTGGGCGACGTGGCCGATGTAAACCTGTGGTTCCAGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCTTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGACAGGCCGCTATCAGGCGTTCACTCAGGCCGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAAcc                                                                                                                                                                                                              >  2:384868/1‑251 (MQ=255)                                       gTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAATGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGCCAAACCGTTCAAAACCGt                                                                                                                                                                                                            >  1:354342/1‑251 (MQ=255)                                           aTGTGGGCGACGTGGTAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCt                                                                                                                                                                                                        <  2:127770/251‑1 (MQ=255)                                                  cgacggggcagaggtaaacccgtgGTTCCTGGAAAATGGCGTTTCCTGCATCTTCAATCTCGGTACTGGTCGTGCGGAAGCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCTCTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGTTGAAGgtg                                                                                                                                                                                                 <  1:363662/230‑1 (MQ=255)                                                           agATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAAt                                                                                                                                                                                        >  1:76121/2‑251 (MQ=255)                                                             atgtaaaccTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGACACGGAATAc                                                                                                                                                                                      >  2:98043/9‑251 (MQ=255)                                                                 aaaccTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATgg                                                                                                                                                                                  <  1:334961/247‑1 (MQ=255)                                                                  aaccTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGTTGTAACGGAATACATGGc                                                                                                                                                                                 <  2:272539/248‑1 (MQ=255)                                                                  aaccTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGc                                                                                                                                                                                 <  1:144672/248‑1 (MQ=255)                                                                       gtgGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCTCTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGc                                                                                                                                                                            <  1:353709/251‑1 (MQ=255)                                                                        tgGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGCATACATGGCCTGGCt                                                                                                                                                                           >  2:320391/1‑251 (MQ=255)                                                                                        gCGTTTCCGGCAGCGTCAAGCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATaa                                                                                                                                                           <  1:371928/251‑1 (MQ=255)                                                                                         cGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATa                                                                                                                                                            >  2:81209/1‑249 (MQ=255)                                                                                         cGTTCCCGGCATCTTCAATCTCGGTTCTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAACCCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAag                                                                                                                                                          <  2:296273/251‑1 (MQ=255)                                                                                              ccGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGctc                                                                                                                                                     <  1:398932/251‑1 (MQ=255)                                                                                              ccGGCATCTTCAATCTCGGTACTGGGCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGct                                                                                                                                                      <  1:9691/250‑1 (MQ=255)                                                                                                ggCATCTTCAATCGCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTg                                                                                                                                                   <  2:295898/251‑1 (MQ=255)                                                                                                     cTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATg                                                                                                                                               >  1:117406/1‑250 (MQ=255)                                                                                                        cAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGaaaa                                                                                                                                           <  1:191285/251‑1 (MQ=255)                                                                                                         aaTCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAAACGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAAt                                                                                                                                          >  1:306783/1‑251 (MQ=255)                                                                                                                    tGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCg                                                                                                                               >  2:381021/1‑251 (MQ=255)                                                                                                                         gagCGGAATCCTTCCAGGCTGAAGCTGATGCTACGCTGGCTTATCACACGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGATACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCg                                                                                                                          >  2:17416/3‑251 (MQ=255)                                                                                                                           gCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTc                                                                                                                        >  1:49400/1‑251 (MQ=255)                                                                                                                            cGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACTACAAACCGTTCAAAACCGTTGCTGAAGGGGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCt                                                                                                                       <  1:20724/251‑1 (MQ=255)                                                                                                                                  ccTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGtt                                                                                                                 >  1:342040/1‑251 (MQ=255)                                                                                                                                                    aTGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACAGGGCCTGGCTGAATCGCGACGCATAAGAGCTCGGCAGGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCg                                                                                               <  2:400511/251‑1 (MQ=255)                                                                                                                                                        tACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaa                                                                                            >  1:111762/1‑250 (MQ=255)                                                                                                                                                         aCGCTGGCTTATCACAAGAAAGGCCAGATCGAATACAGTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCATATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaaa                                                                                           <  2:139639/250‑1 (MQ=255)                                                                                                                                                              ggCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTctct                                                                                     <  2:145751/251‑1 (MQ=255)                                                                                                                                                                      acaAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAATGCCGCTATAAGGCGTTCACTCAGGCAGATCTGACGACTCTGCGCGCGGCGGGTTACGCCAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGCCTCTATCACACg                                                                             >  2:73252/1‑251 (MQ=255)                                                                                                                                                                        aagaaagaCCAGATCGAATACATTCCGTTACCGGATAAACTAACAGGCCGCTATCAGGCGTTCACTCAGACAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGCAAATACTGGTGATCGGCCCGTCTTGGGTTg                                                                                                                >  2:285303/1‑214 (MQ=255)                                                                                                                                                                        aaGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACTACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTg                                                                                                                <  1:285303/214‑1 (MQ=255)                                                                                                                                                                           aaaGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCCCGCTcc                                                                         >  2:427740/1‑250 (MQ=255)                                                                                                                                                                                      cGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACGCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGTCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGGCCGTCTTGGGTTGGCGACATGATTATGTCGCTAAGTCTCTAGCGCTCGCTCCAGGCGCGCTAt                                                              <  2:182240/250‑1 (MQ=255)                                                                                                                                                                                           aCATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATcccc                                                          >  1:373381/1‑249 (MQ=255)                                                                                                                                                                                           aCATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAg                                                        <  2:342040/251‑1 (MQ=255)                                                                                                                                                                                           aCATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCAATCAGGCAGATCTGAAGAATCTGCGCGCGGCGGGTTACGACAAACCGTCCAAAACCGTTGCTGCAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGCATTGGGTTGGCGACATGTTGATGTCGCAAAGACTCTATCGCACGCGCCAGGCGCGCTATcccc                                                          >  2:241844/1‑249 (MQ=255)                                                                                                                                                                                                 cGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGAGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGAGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGAGATATACCCAGGCGAt                                                   >  1:255248/1‑250 (MQ=255)                                                                                                                                                                                                  gTTCCCGGATAATCTGACAGGCCGATATCAGGCGTTCACTCCGGCCGATCTGACGAATCTGCGCGCGGCGGCATACGACAAACCGTTCAAAACCGTTGCTGAAGGCGTACCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGACGATGTCGCAAAGTCTCTATCGCACGCTCCCGGCGCGCTATCCCAAGGCGATa                                                  >  2:21026/1‑250 (MQ=255)                                                                                                                                                                                                  gTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCCGGCAGATCTGACGAATCTGCGAGCGGCGGGTTACGACAACCCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCCCCAGGCGCGCTATCCCCAGGCGATa                                                  >  2:182948/1‑250 (MQ=255)                                                                                                                                                                                                     cccGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCg                                              <  2:6719/251‑1 (MQ=255)                                                                                                                                                                                                        ggATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAtg                                           >  2:297507/1‑251 (MQ=255)                                                                                                                                                                                                                  aaaGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCAcc                                 <  2:6636/251‑1 (MQ=255)                                                                                                                                                                                                                    aGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCgg                               >  1:4935/1‑251 (MQ=255)                                                                                                                                                                                                                         gCTATCAGGCTTGCACTCAGGCAGAGCTGACGACGCTGCGCGCGGCTGGTTACGACAAACCGTTCAAATCCGATGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCTCGACGCATAAGAGCTCTGCATGAAAATACTTGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACTCTCCAGGCGCTCTATCCCCAGGCGTTCATCGATGTGATGGCACTGGcctgt                          <  2:165330/251‑5 (MQ=255)                                                                                                                                                                                                                                   gTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGGTGATGTCGCAAAGACTCTATCGCACGCTCCAGGCGCGCTATCCCCCGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat                >  2:135663/1‑251 (MQ=255)                                                                                                                                                                                                                                    ttCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCCGAAGGGGTAACGGAAGACATGGACTGGATGAATCGCGCCGCATAATAGCTCTGCATCAAAATACTGGTGATCGGCCCGACTTAGGTAGGCGACATGATGATGTCGCAATGGCTCTATCGCACGCTCGCGGCGCGCGATCCCCACGCGCTAATCGGTGTGATGGCACCGGCATGGTGCCGTCCatt               >  1:197916/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aTCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATg  >  1:220569/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    AGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCTGATGTGAATCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCGGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATG  >  NC_000913/3794582‑3795072

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑