Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,512,642	T→C	V256V (GTA→GTG)	fecC ←	ferric citrate ABC transporter membrane subunit FecC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,512,642	0	T	C	97.9%	142.8 / NA	47	V256V (GTA→GTG)	fecC	ferric citrate ABC transporter membrane subunit FecC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (25/21);  minor base A (1/0);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCA  >  NC_000913/4512401‑4512856                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                       aGCGCCATTTTCCTCCTCGCCTCCTCTCAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCGCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCATCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCCCCAGCAGACCGATAACCGCTTCCGGACCTGCCTCGCTGacg                                                                                                                                                                                                               <  2:214652/251‑1 (MQ=255)               ccTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCaccaacca                                                                                                                                                                                                  <  1:87561/250‑1 (MQ=255)                  cGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGa                                                                                                                                                                                                >  1:224572/1‑249 (MQ=255)                    cctcctCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAagc                                                                                                                                                                                             >  1:94214/1‑250 (MQ=255)                       cctcACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAagcagcc                                                                                                                                                                                         >  1:101190/1‑250 (MQ=255)                         ccacaAGCCAGAAAAATCAAGGGCCTCCAAGCAGCGCCAGCACGGCGCCTGCGGGCAGATCCCCGGGGAAGGCCAGTGCGCG‑GCGGAGTCCATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCGCCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACt                                                                                                                                                                                       <  2:118847/250‑1 (MQ=255)                          cccaAGCCAGACAAAGCAAGGGCTGCCAACCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACGGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCAc                                                                                                                                                                                        <  2:180205/247‑1 (MQ=255)                                   gACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCGCTAACAGATtga                                                                                                                                                                             <  1:226145/251‑1 (MQ=255)                                    aCAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtgat                                                                                                                                                                            <  1:133205/251‑1 (MQ=255)                                     aaaaGCAAGGGCTGCCCATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTCACAAATtgatg                                                                                                                                                                           >  2:809/2‑251 (MQ=255)                                      aaaGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtgatgc                                                                                                                                                                          >  1:139916/1‑250 (MQ=255)                                      aaaGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtgatg                                                                                                                                                                           >  1:83656/1‑250 (MQ=255)                                      aaaGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACAAGAAGCAGCACTAACATATtgatg                                                                                                                                                                           >  2:90797/1‑250 (MQ=255)                                      aaaGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGAGAAACGCCACCGGTCCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtgatga                                                                                                                                                                          <  1:265093/251‑1 (MQ=255)                                           aaGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACCTCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCCCTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCCGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATCTTGATGACCaa                                                                                                                                                                      >  1:20579/1‑250 (MQ=255)                                               gCTGCCAAGCAGCGCCAGCCCTGGGCCGGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCAGCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTa                                                                                                                                                                 <  2:79024/251‑1 (MQ=255)                                                      aTCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACCTATTGATGACCAAACGTAGCCt                                                                                                                                                             >  2:169856/1‑248 (MQ=255)                                                      aTCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCt                                                                                                                                                             <  1:169856/248‑1 (MQ=255)                                                        cAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTc                                                                                                                                                         >  1:63846/1‑250 (MQ=255)                                                        cAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTc                                                                                                                                                         >  1:14165/1‑250 (MQ=255)                                                             gcAGCACTGCGCCGGAGGGCAGATCTCCGGGGAAGGCCCGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCCAAATGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGAGGCACCAGCAGACCGATAAACGCCACCTGACCTGCAAAGCTGACGCACGCACAAACCAGAAGCGGCACTAACATATTGATGACCAAACGAAGCCACGTCCGGttt                                                                                                                                                   >  2:164866/2‑250 (MQ=255)                                                                    tGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTccc                                                                                                                                             >  1:28557/1‑250 (MQ=255)                                                                    ggcgcCTGCGGGCAGATCTCCGGGGAAGGCCGGTGCGCGGGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACGGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCa                                                                                                                                            <  2:14165/250‑1 (MQ=255)                                                                    ggcgcCAGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCa                                                                                                                                            <  2:66426/250‑1 (MQ=255)                                                                       gcCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCg                                                                                                                                         <  2:210118/251‑1 (MQ=255)                                                                         cTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGt                                                                                                                                        >  1:104308/1‑250 (MQ=255)                                                                                           ggAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTg                                                                                                                     >  1:192851/1‑251 (MQ=255)                                                                                            gAAGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACTCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGa                                                                                                                    <  1:182766/251‑1 (MQ=255)                                                                                             aaGGCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGAGCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCTCTGAg                                                                                                                   <  1:214170/251‑1 (MQ=255)                                                                                               ggCCAGTGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGt                                                                                                                 >  1:246877/1‑251 (MQ=255)                                                                                                            gCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTc                                                                                                      >  1:256104/1‑249 (MQ=255)                                                                                                                         gCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCCGGTTCAGTTGATTCGcc                                                                                          >  1:236456/1‑249 (MQ=255)                                                                                                                          ccAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCcagc                                                                                       >  2:67479/1‑251 (MQ=255)                                                                                                                                    cAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTTCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCCCGCACCAACCAGAAGCCGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAgcagcacaacg                                                                              >  2:215412/1‑245 (MQ=255)                                                                                                                                                      caTGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGACGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCCACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGacc                                                            >  2:101409/1‑250 (MQ=255)                                                                                                                                                          cccacgggCAGGACGTTGCGCTGAGCGAAGCCCGCCCAGAAGCGCGCCAGATGTGGCACCAACAGAACGAAAAACGCCCCCGGGCCTTCCACGCTGACGGCCGCACCAACCAGAAGCAGCACCAACAAATTGATGACCAAACGTAGGCCCGGCAGGGTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGaccaccac                                                       <  2:101190/245‑1 (MQ=255)                                                                                                                                                                             gcTGATCGAAGCCTGCCCAGAAGCGCGCCAGCTGTGGCCCCAGCAGACCGATAAGCGCAACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACGGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGAc                                     >  2:110341/1‑250 (MQ=255)                                                                                                                                                                              cTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACAt                                   >  1:171925/1‑251 (MQ=255)                                                                                                                                                                                 aTCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGTCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCt                                <  1:125286/251‑1 (MQ=255)                                                                                                                                                                                 aTCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCCACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATccc                                >  2:213400/1‑250 (MQ=255)                                                                                                                                                                                    gAAGCCAGCCCAGAAACGCGCCAGATGTGGCACCAGCAGACCGATAACCGCCACCGGCCCTGCCACGCTGACGCACGCACCCACCAGAAGCAGCACTAACATACTGATGACCAAAAGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCCGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGcc                             >  2:39410/1‑251 (MQ=255)                                                                                                                                                                                        ccTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAAACAGAAGCAGCACTAACATCTTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGCGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCg                         >  1:31632/1‑251 (MQ=255)                                                                                                                                                                                           gCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGGCAGCggg                       <  1:196482/250‑1 (MQ=255)                                                                                                                                                                                                          aGATGGGGCACCGGCACGCCGAAAAACCCCACCGACCCTGCCTCGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGAGGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGGCCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCAGcc       <  2:204312/251‑1 (MQ=255)                                                                                                                                                                                                              gtgGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTgcc   <  1:265133/251‑1 (MQ=255)                                                                                                                                                                                                               gggCACCAGCGGACCGATAAACGCCACCGGCCCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATAATGATGACCAAACGTAGCCTCGTCAGGGTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTTCCAGCGGGCTTGGGACACTCCGCCTgcca  <  2:240478/250‑1 (MQ=255)                                                                                                                                                                                                                ggCACCAGCAGACCGATAAACCCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTgcca  <  1:125528/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                       AGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCG‑AGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCA  >  NC_000913/4512401‑4512856

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