Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,515,348	A→C	L444V (TTG→GTG)	fecA ←	ferric citrate outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,515,348	0	A	C	100.0%	170.7 / NA	53	L444V (TTG→GTG)	fecA	ferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (22/31);  total (22/31)

CGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATG‑ACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATAC  >  NC_000913/4515106‑4515592                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     gggttcGTTAGAGCTCTCGTCTTCGTTCGGGCCTGTTATGGCGTTGTTCGTGTATGCAC‑CGATATGATCGAAACTCATACCCGGCGTGAGGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAATTGCCCGCTGCTGGTTGTGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCCCATAGCGAt                                                                                                                                                                                                                                              <  2:380553/246‑1 (MQ=255)     gCGTTATAGCTCACTTCTTCGTTCGTGCCTGTGATGGCGTTGTTCTGGTATG‑ACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAgcc                                                                                                                                                                                                                                          <  1:78510/251‑1 (MQ=255)          attgCTCGCTTCGTCGTGCGTGCCTGTGATGGCGTTGTTCATGTATT‑ACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTAGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGTAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTTGCGGTGTAGTAACTC‑ATTTCATGCGTTGATTCATTCACATAGCGATAgcccacgc                                                                                                                                                                                                                                     <  1:71283/248‑1 (MQ=255)          attgCTCCCTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTTGGATG‑ACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGGCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAgcccacgc                                                                                                                                                                                                                                     <  2:81874/248‑1 (MQ=255)              cTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATG‑ACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAgcccacgcccac                                                                                                                                                                                                                                 >  2:63833/1‑251 (MQ=255)                         gtgcgtacCTGTGATGGCGTTGTTCTAGTATG‑ACTCGATATGTTAGAAACGCCTACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAAATGCCCGCTGCTGGTGGCGGGGGAGTAAAGC‑ATTTCATGCGTTGATTCATTCACGTAACGATAGCACACGCACACTTCgacggcgg                                                                                                                                                                                                                      >  2:60236/1‑244 (MQ=255)                          tgcgtgcCTGTGCTGGCGTTGTTCTGGCATG‑ACTCGATATGTTCGCAACGCATACCCGGCGTGATGGTACAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGCGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACCGCCACTGCCCGATGCTGGTGGCGGAGTAGTCACGA‑AATTCATTCGTTCATTCATTCCCATAGCGATAGACCACGCCCACTTCGTGCGCgg                                                                                                                                                                                                                      >  2:21858/1‑250 (MQ=255)                               gcCTGTGATGGCGTTGTTCTGGTATG‑ACTCGATCTGTTCGAAACGCCTACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGCTACCAGGCGTGCGCATCGGTGCCGGAACGCGTATCGAGGTCGTAAGGGCTTGAGCCGGACGGCACCTGCCCGCTGCTGGTGGCGGGGTAGTAACGC‑ATTTAATGCGTTGATTCATTCACAAAACGATAAACCACCACCCCTTCGTGCGCGGAAGGGc                                                                                                                                                                                                                >  2:258794/1‑251 (MQ=255)                                ccTGTGATGGCGTTGTTCTGGTATG‑ACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGcc                                                                                                                                                                                                               <  1:276605/251‑1 (MQ=255)                                   gtgATGGCGTTGTTCTGGTATG‑ACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGa                                                                                                                                                                                                             <  1:211589/250‑1 (MQ=255)                                   gtgATGGCGTTGTCCTTGTATG‑ACTCGATATGTTCGAAACGCATCCCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCTCGCCCTCTTCGTGCGCGGAAGGGCCg                                                                                                                                                                                                              <  2:291424/249‑1 (MQ=255)                                         gCGTTGTTCTGGTATG‑ACTCGATATGATCGAATCGCATACCCGGCGTGATGGTCCAGTTGCCGATGTTGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGGGTATCGCGGGCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATaaa                                                                                                                                                                                                      <  2:272795/251‑1 (MQ=255)                                             tgttCTGGTATG‑ACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGGTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATaa                                                                                                                                                                                                       >  1:212331/1‑246 (MQ=255)                                              gttCTGGTATG‑ACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCt                                                                                                                                                                                                 <  1:63833/251‑1 (MQ=255)                                                       atgaaccGATATTTACGAAACGCATACCCGGCGTGATGGTCCAGTGGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCCCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTa                                                                                                                                                                                          <  2:134644/244‑1 (MQ=255)                                                             cgctatGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGGGGTGTAAGAACGC‑GTTTCGTGCGTTTATTCATTCACATAGCGGTATCCACCTCCCACTTCGGTCGCGGGAGTGCCAGTAATCAAGGTCTGGGTGTAGCGTGGc                                                                                                                                                                                   >  2:297838/4‑251 (MQ=255)                                                                 atGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCaa                                                                                                                                                                               <  1:292051/251‑1 (MQ=255)                                                                  tGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGATGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGGTTCATTCACCTAGCGATAGCCCACACCCACTTCGTGCGCGGTAGGGCCGATCATAAAAAACTGGCAGTAGCGTGGCTccat                                                                                                                                                                              >  2:94872/1‑248 (MQ=255)                                                                                   cGGCTTGAGGGTCCTGTTGCCTAGGTCGATTTTGTCATCCTGATACCTGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTTCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGTTAGCCCACGCCCTCTTCGTGCTCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATCCCGCGCTCCCagtagt                                                                                                                                                             <  2:257477/251‑1 (MQ=255)                                                                                       gTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTAcg                                                                                                                                                         <  1:373322/251‑1 (MQ=255)                                                                                           gggTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGTTGCCGGAACGCGTATCGCGGGCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGTTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACgcggc                                                                                                                                                     <  1:316453/250‑1 (MQ=255)                                                                                                     gCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGACCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTgg                                                                                                                                           >  2:278778/1‑250 (MQ=255)                                                                                                             cGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTg                                                                                                                                   <  1:60236/251‑1 (MQ=255)                                                                                                               aTTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAGCGCGTATCGCGGTCGTAAGGACTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTCGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCCCGCCCACTTCGTGCGCGGACGGGCCGATCATAAAGATCTGACTGTAGCGTGGCTCAATACCGCGCACCCCGTAGTTACGCGGCGCGAGGGTGATGCGTTTGcc                                                                                                                                 >  1:199746/1‑251 (MQ=255)                                                                                                                   tGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTg                                                                                                                              <  1:379484/250‑1 (MQ=255)                                                                                                                      cATCCAGATACCAGGGGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGTGCGTGATCCGGACGGCAACGGCCCGCGGCTTGTGGCGGTGGAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACTCCCACTTCGTGCGCGGAAGGGTCGATCATAAAGCTCTTGTGGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTcc                                                                                                                          <  2:74890/251‑1 (MQ=255)                                                                                                                            gATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAAAGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCACACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGCGAGGGGGCTGCGTTTGCCTTGCTCCAGGt                                                                                                                      >  1:196593/1‑249 (MQ=255)                                                                                                                             aTACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGc                                                                                                                   >  1:280290/1‑251 (MQ=255)                                                                                                                                       tGCGCCTCGGTGCCGGAACGCGTATCTCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGGGGCAGGGCCGATCATAAAGATCTGGCTTTAGCGTGGCTCAAtcccccggac                                                                                                                                                                     >  2:256695/1‑182 (MQ=255)                                                                                                                                       tGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCAc                                                                                                                                                                     <  1:256695/191‑1 (MQ=255)                                                                                                                                          gccccggggCCGGAACGCGTTTCGCGGTCGTAAGGTCTTGAGCCGGACGGCACCTGCCCGCTGCTGGGGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATGCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGGTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGt                                                                                                      <  2:146471/243‑1 (MQ=255)                                                                                                                                            cTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATc                                                                                                                                                                                                           >  1:279179/1‑148 (MQ=255)                                                                                                                                            cTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACGGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTGCGTGCGCGGAAGGGCCGATc                                                                                                                                                                                                           <  2:279179/148‑1 (MQ=255)                                                                                                                                                           ggtttcgcgggctgaaGGTCTTGAGCTGGACGGCTCCTGCCAGTTACTGCTGGCGATGAAGTAACGCTATTT‑AGGCGTGGATTCAGTCACACACCGATAGGGCCCGCTCACTTTGTGTGCGGAGTGGCCGATCATAAAGATCTGGCTGTTGCGTGTCTCAATAGCGCGCAGCCAGTAGATACGCGGCGATAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCTGGTTTTGGGTGGAGAActcat                                                                                     <  2:30234/237‑5 (MQ=255)                                                                                                                                                           cgTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCt                                                                                     >  1:230329/1‑251 (MQ=255)                                                                                                                                                            gtaacgcgTTCGGAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTTGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACGTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTg                                                                                    <  2:195884/247‑1 (MQ=255)                                                                                                                                                              aTCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGGAGGGCCGATCATAAAGATCTTGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGGTACGCGGGGAGAGGGTGATGCGTTTGCCTTTCTCCAGGGAGCCGATGCGGAGGGTTTGGGTGTAGGACCCCTGa                                                                                   >  2:144120/1‑250 (MQ=255)                                                                                                                                                                      cGTAAGGGCTTGAGCCGGACGGCAATTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGTCCGATCATAAATATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAAt                                                                          <  1:324711/251‑1 (MQ=255)                                                                                                                                                                       gTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGGAGAACCCCTGAATGTTGAAtt                                                                         >  1:361090/1‑251 (MQ=255)                                                                                                                                                                          aGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTAt                                                                      >  1:197728/1‑251 (MQ=255)                                                                                                                                                                             gCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCAGAATGTTGAATTTATgct                                                                   <  1:96246/251‑1 (MQ=255)                                                                                                                                                                                 gAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATgctggc                                                                <  1:364873/250‑1 (MQ=255)                                                                                                                                                                                    cctgACGGCAACTGCCCTCGGCTGGTGGCGGTGTAGTACCGC‑ACTTCATGCGTGGATTCATTCACATAGCGATAGCCCACGCCCTCTGCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCCCCCAGTAGGTACGCGGCGAGAGGGTGATGCGGTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTc                                                            <  2:61425/248‑1 (MQ=255)                                                                                                                                                                                          ggCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGCATTTATGCTGGCTGTctggctg                                                      >  1:206185/1‑251 (MQ=255)                                                                                                                                                                                                        gctgGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCTCTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGccc                                         <  1:199691/250‑1 (MQ=255)                                                                                                                                                                                                        gctgGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGccc                                         <  1:345868/250‑1 (MQ=255)                                                                                                                                                                                                              ggCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTc                                   <  1:338779/250‑1 (MQ=255)                                                                                                                                                                                                                         aaCGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGTCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCGCCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGcga                       <  2:60640/251‑1 (MQ=255)                                                                                                                                                                                                                          aCGC‑ATTTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGcgac                      <  1:278778/251‑1 (MQ=255)                                                                                                                                                                                                                           cGC‑ATGTCATGCGTTGATTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGAGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGCCCCGCCGACTTGcgacg                     >  2:168327/1‑251 (MQ=255)                                                                                                                                                                                                                            gC‑ATTTAATGCGTTGATTCATTAACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAACGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGCGAGGGTGATGAGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTTCATTTATGCTGGCAGTCTGGCTGGAACTGGCAGATCAGGCTCGCCAGcgtggaacg                     >  2:235678/1‑242 (MQ=255)                                                                                                                                                                                                                               aTTTCATGCGTTGCTTCATTCACATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGAcc                  >  1:82511/1‑251 (MQ=255)                                                                                                                                                                                                                                               cattcaCATAGCGATAGACCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCACTACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGCGTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACGGGTAGCCCAGGCTCGCCAGCGTGCGACGACCCCAGAAGCGGTCCTAc  >  2:104835/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     CGGTGCGTTATAGCTCACTTCTTCGTGCGTGCCTGTGATGGCGTTGTTCTGGTATG‑ACTCGATATGTTCGAAACGCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGC‑ATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATAC  >  NC_000913/4515106‑4515592

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