Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,794,649	T→C	N221N (AAT→AAC)	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,649	0	T	C	100.0%	257.3 / NA	74	N221N (AAT→AAC)	rfaD	ADP‑L‑glycero‑D‑mannoheptose 6‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/34);  total (40/34)

CTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCTGATGTGAATCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCGGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAAT  >  NC_000913/3794403‑3794890                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  cTACGGATACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCtgt                                                                                                                                                                                                                                               >  1:127770/1‑251 (MQ=255)   aCGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTTTTTGGGCGACGTGGCAGATGTAAACCTGTgg                                                                                                                                                                                                                                             >  2:113690/1‑251 (MQ=255)      gTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTc                                                                                                                                                                                                                                          >  1:400511/1‑251 (MQ=255)             aaaaTTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGaaa                                                                                                                                                                                                                                   <  2:380609/251‑1 (MQ=255)              aaaTTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTTAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCCGATGTAAACCTGTGGTTCCTGGaaa                                                                                                                                                                                                                                   >  1:142176/1‑250 (MQ=255)                      gTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGttt                                                                                                                                                                                                                          <  2:313521/251‑1 (MQ=255)                           atgaATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGc                                                                                                                                                                                                                     <  2:234911/251‑1 (MQ=255)                               atatGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGACGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCt                                                                                                                                                                                                                 >  1:295898/1‑251 (MQ=255)                                 atGTTCGTCAAATCCTGCCGGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTGCCGGCATCTTc                                                                                                                                                                                                               <  2:71251/251‑1 (MQ=255)                                     tCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAAt                                                                                                                                                                                                            <  2:253481/250‑1 (MQ=255)                                          aaaTCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGt                                                                                                                                                                                                      >  1:272539/1‑251 (MQ=255)                                                 gcctgAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTc                                                                                                                                                                                               >  1:35656/1‑251 (MQ=255)                                                  cctgAAGCGATCTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCg                                                                                                                                                                                              >  1:6890/1‑251 (MQ=255)                                                           aaCTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATc                                                                                                                                                                                     <  1:371635/251‑1 (MQ=255)                                                             cTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATCGCGAGCGTCGCTTTCCATCTCAACAACCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACATGTGGATCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCt                                                                                                                                                                                   >  2:271855/1‑251 (MQ=255)                                                                   gATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAgg                                                                                                                                                                             >  1:6719/1‑251 (MQ=255)                                                                          ggCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGc                                                                                                                                                                      <  1:178375/251‑1 (MQ=255)                                                                          ggCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGc                                                                                                                                                                      <  1:330023/251‑1 (MQ=255)                                                                           gCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCt                                                                                                                                                                     <  1:60772/251‑1 (MQ=255)                                                                                cGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGc                                                                                                                                                                >  2:373849/1‑251 (MQ=255)                                                                                    aTTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACGTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTTGAAAATGGCGTTTCCGGCATCGTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTAc                                                                                                                                                             <  1:273416/250‑1 (MQ=255)                                                                                          aCGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTGGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCt                                                                                                                                                      <  2:34230/251‑1 (MQ=255)                                                                                            gTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCtt                                                                                                                                                     >  2:261699/1‑250 (MQ=255)                                                                                             tttATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGATAACTTCAAACGCGTCTTCGTCTATGTGGGCGCCGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTa                                                                                                                                                    <  2:109625/250‑1 (MQ=255)                                                                                              ttATGGACCGAGTGAAGGCCATAAAGGCAGTATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTAt                                                                                                                                                   <  2:142176/250‑1 (MQ=255)                                                                                                aTGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGTCTTATcac                                                                                                                                                <  1:113690/251‑1 (MQ=255)                                                                                                aTGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATcac                                                                                                                                                >  2:378443/1‑251 (MQ=255)                                                                                                   gACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCATAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAAGGGCGTTTCCGGCATCTTCAATCTCGGTACTTGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTGCGCTGGCGTATCACAAg                                                                                                                                             <  1:130714/251‑1 (MQ=255)                                                                                                           gAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGcc                                                                                                                                      <  1:353819/250‑1 (MQ=255)                                                                                                            aaGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTCGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGcc                                                                                                                                      >  1:405211/1‑249 (MQ=255)                                                                                                               gCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGcc                                                                                                                                      >  1:296273/1‑246 (MQ=255)                                                                                                                ccATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGcc                                                                                                                                      >  1:167150/1‑245 (MQ=255)                                                                                                                ccATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCg                                                                                                                                <  1:256821/251‑1 (MQ=255)                                                                                                                ccATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGCTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACACGaaacacc                                                                                                                                      >  2:302839/1‑241 (MQ=255)                                                                                                                  aTAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCCCCAACCAGCTTAACACCGGAGAATCACCGAAGCTGTTTGAAGGTAGCGAGAAATTAAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGAAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTCCTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGa                                                                                                                               >  2:226983/1‑250 (MQ=255)                                                                                                                  aTAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGACGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGCCGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGa                                                                                                                               >  2:311053/1‑250 (MQ=255)                                                                                                                   tAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAAt                                                                                                                             >  2:191285/1‑251 (MQ=255)                                                                                                                          agcaTGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTcc                                                                                                                      >  2:433806/1‑251 (MQ=255)                                                                                                                               ggCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTcc                                                                                                                 >  1:4531/1‑251 (MQ=255)                                                                                                                                 cGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCg                                                                                                               >  2:318977/1‑251 (MQ=255)                                                                                                                                     cGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGCAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGAt                                                                                                            >  2:371928/1‑250 (MQ=255)                                                                                                                                      gTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTGCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTTCGCTGGCTTATCACAAGAAAGGCCAGATCGAATGCTTTCCGTTCCGGGATaa                                                                                                          <  2:205347/251‑1 (MQ=255)                                                                                                                                        cGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAAc                                                                                                        <  1:373849/251‑1 (MQ=255)                                                                                                                                         gCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACt                                                                                                       <  1:354300/251‑1 (MQ=255)                                                                                                                                           tttCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTg                                                                                                      >  1:174683/1‑250 (MQ=255)                                                                                                                                              ccATCTCAACACCCAGTTTAACAACGGTGAATCCCCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGaaa                                                                                                   <  2:108466/250‑1 (MQ=255)                                                                                                                                              ccATCTCAACACCCAGCTTACCAACGGTGAATCGCCGAAGCTGTTTGAAAGTATCGAGAACTTCAAACGCGACTTCGTCTATGTGTGCGACGTGGCAGATGTAAACCTGGGTTTCCTGGAAAATGGCGTTACCGGCATCTTCAACCTCGGTACTGGTCGTGCGGAATCATTCCATGCTGTAGCTGTTGCTACGCTGGCTTATCACAAGAAAGGCCAGATCAAATGCTTTCCGTCCCCGGATAAACTGAAAg                                                                                                  <  2:282498/251‑1 (MQ=255)                                                                                                                                               cATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAgg                                                                                                 >  2:34385/1‑251 (MQ=255)                                                                                                                                                         cccAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTGCCTGGAAAATGGCGTTTCCGGCATGTGCAATCGCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCGACGCTGGCTTATCACAAGAAAGGGCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAg                                                                                       <  1:226983/251‑1 (MQ=255)                                                                                                                                                            aGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCt                                                                                                                                                          >  1:369121/1‑181 (MQ=255)                                                                                                                                                            aGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGGAGCTGAGGCTACGCt                                                                                                                                                          <  2:369121/181‑1 (MQ=255)                                                                                                                                                              cTTAACACCGGTGAATCACCGAAGATGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCAATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGCTCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTCTCAGGCGtt                                                                                  >  2:21798/1‑251 (MQ=255)                                                                                                                                                                      cGGTGAATCACCGCAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAACCCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAgg                                                                          >  1:305212/1‑251 (MQ=255)                                                                                                                                                                      cGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAgg                                                                          >  1:6636/1‑251 (MQ=255)                                                                                                                                                                            aTCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATc                                                                    <  2:167150/251‑1 (MQ=255)                                                                                                                                                                             tCACCGAAGCTGTTTGACGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCt                                                                   >  2:261329/1‑251 (MQ=255)                                                                                                                                                                                 cGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGAcg                                                               >  1:58297/1‑250 (MQ=255)                                                                                                                                                                                    aGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGCCTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAAt                                                            >  1:145751/1‑251 (MQ=255)                                                                                                                                                                                     gCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATc                                                           <  1:415910/251‑1 (MQ=255)                                                                                                                                                                                        gTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTgc                                                        <  1:313833/251‑1 (MQ=255)                                                                                                                                                                                          ttGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTg                                                                                                                                                                                            >  1:213290/1‑117 (MQ=255)                                                                                                                                                                                          ttGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTg                                                                                                                                                                                            <  2:213290/117‑1 (MQ=255)                                                                                                                                                                                            gAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGTGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGa                                                                                                                               <  1:59408/176‑1 (MQ=255)                                                                                                                                                                                            gAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCAGATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGa                                                                                                                               >  2:59408/1‑176 (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)                                                                                                                                                                                                                                              agATGTAAACCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCTGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACGAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAAt  >  1:76121/2‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  CTACGGTTACTCAAAATTCCTGTTTGATGAATATGTTCGTCAAATCCTGCCTGAAGCGAACTCGCAGATTGTTGGCTTCCGTTATTTCAACGTTTATGGACCGCGTGAAGGCCATAAAGGCAGCATGGCGAGCGTCGCTTTCCATCTCAACACCCAGCTTAACAACGGTGAATCACCGAAGCTGTTTGAAGGTAGCGAGAACTTCAAACGCGACTTCGTCTATGTGGGCGACGTGGCTGATGTGAATCTGTGGTTCCTGGAAAATGGCGTTTCCGGCATCTTCAATCTCGGTACTGGTCGTGCGGAATCCTTCCAGGCGGTAGCTGATGCTACGCTGGCTTATCACAAGAAAGGCCAGATCGAATACATTCCGTTCCCGGATAAACTGAAAGGCCGCTATCAGGCGTTCACTCAGGCAGATCTGACAAATCTGCGCGCGGCGGGTTACGACAAACCGTTCAAAACCGTTGCTGAAGGTGTAACGGAAT  >  NC_000913/3794403‑3794890

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