Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,084,932	G→T	F297L (TTC→TTA)	fdoG ←	formate dehydrogenase O subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,084,932	0	G	T	100.0%	212.9 / NA	61	F297L (TTC→TTA)	fdoG	formate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (26/35);  total (26/35)

CAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACC  >  NC_000913/4084694‑4085172                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 cAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTg                                                                                                                                                                                                                                      >  1:234331/1‑251 (MQ=255) cAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGAGCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTg                                                                                                                                                                                                                                      <  1:150500/251‑1 (MQ=255) cAACCACATCAGGCGAGTAGCGGGAAACGAGCTGTTTCAGCACGTTCCACACGCCGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGACCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGCACAGGCCATCTTCAAAGCCGTAATCCTCACGAACGTTCCGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTg                                                                                                                                                                                                                                      >  2:291165/1‑251 (MQ=255)    ccACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTc                                                                                                                                                                                                                                   >  2:154084/1‑251 (MQ=255)          cTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAgg                                                                                                                                                                                                                             >  1:420505/1‑251 (MQ=255)            ggCGTGTAGGGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTa                                                                                                                                                                                                                           <  1:118980/251‑1 (MQ=255)            ggCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTGGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTGCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGTTATAGGCTTCGGAGTATTCGCGGTTTAATTTTTCATTGTGCAGCAgg                                                                                                                                                                                                                             <  2:67323/249‑1 (MQ=255)                 gtAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCa                                                                                                                                                                                                                      <  1:280475/251‑1 (MQ=255)                     cgtgAAACGTGCTGTGTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGGCCATTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCAGCTGCAAAGCCGTAAGCCTCACGCACGATCAGGCTGGCGTTGGGATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTTTTCAGCAGTTACAGCAATAc                                                                                                                                                                                                                  <  2:250864/248‑1 (MQ=255)                        ggaaCGTGCTGTTTCAGCAAGTTCCACACGCAGTGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTGTCGTCCAGTTCATAAGTCCGGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTGCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTGGGTATAGGCGTCGGTGTATTCGCGGTGTAATTTTTCATTGTGCAGGAGGGACAGCAATACGcc                                                                                                                                                                                                               <  2:209766/249‑1 (MQ=255)                        gAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGcc                                                                                                                                                                                                               >  1:303779/1‑251 (MQ=255)                         aaaCGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGTATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCTTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCt                                                                                                                                                                                                              <  1:211798/251‑1 (MQ=255)                         aaaCGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGA‑CAGGCCCTCTTCCAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAACAATAAGCCTg                                                                                                                                                                                                             >  1:47125/1‑251 (MQ=255)                             gTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGACGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGGTCAGCAGGTACAGCAATACGCCTGAc                                                                                                                                                                                                           >  2:232086/1‑250 (MQ=255)                             gTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGGAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACa                                                                                                                                                                                                          <  2:379652/251‑1 (MQ=255)                               gCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAgg                                                                                                                                                                                                        <  2:98782/251‑1 (MQ=255)                                                cacaCGCCGCGCGGTTGTTGCAGCGTGGTATCGCGTTTGGAAAAGCCGTTTTCGTCCAGTTCATACGTCCAGGTGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCACAACCGTAATCATCACGCACGCTCTGGCTGGCGTTGGCATCGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTAcc                                                                                                                                                                                       >  1:331538/1‑251 (MQ=255)                                                 acacGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCg                                                                                                                                                                                      <  1:87544/251‑1 (MQ=255)                                                     gCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAAc                                                                                                                                                                                  >  1:238901/1‑251 (MQ=255)                                                      cAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACg                                                                                                                                                                                 >  1:398068/1‑251 (MQ=255)                                                         cgcgGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAAt                                                                                                                                                                              >  1:95213/1‑251 (MQ=255)                                                                   gCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTGCATAAGTCCAGGAGGATTTATCGGACTGGCGTTTTTCCGCTTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTCCCGGAACGAATAGGGGCAtag                                                                                                                                                                    <  2:37673/251‑1 (MQ=255)                                                                             aTCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTa                                                                                                                                                          <  1:359787/251‑1 (MQ=255)                                                                             aTCGCGTTTGGCGAAGCCGTTTTCGTCCAGTACATAAGTCCAGGAGGCTTTCTCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAAGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGATGGCGTTGGTATAGGCTTCGGTGAATTCGCGATTTAATTTTTCATTGTTCAGCAGGTACAAAAATACGCCTGACAGGCAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCt                                                                                                                                                           >  1:224889/1‑250 (MQ=255)                                                                               cgcgTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAcc                                                                                                                                                        <  2:169716/251‑1 (MQ=255)                                                                                         gAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAc                                                                                                                                              >  1:118802/1‑251 (MQ=255)                                                                                         gAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAc                                                                                                                                              <  2:64125/251‑1 (MQ=255)                                                                                         gAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAc                                                                                                                                              >  1:87756/1‑251 (MQ=255)                                                                                         gAAGCCGTTTTCGTCCAGTTCATAAGGCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAc                                                                                                                                              <  2:240448/251‑1 (MQ=255)                                                                                               gTTTTCGTCCAGTTGATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCATTTCCGGATCGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTaa                                                                                                                                        <  1:368012/251‑1 (MQ=255)                                                                                                  ttCGTCCAGTTCATAAGTCCAGGAGTATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAACGCCGTAATCCGCTCGCACGATCAGGCTGGCGTTGGTATAGGCTGAGGTGTATTCGCGGGTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATTGCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTaaat                                                                                                                                      <  2:230752/250‑2 (MQ=255)                                                                                                    cGTCCAGTTCAGAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCTTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAAGACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCAAAGTAGGCAGTTACCGCCGCCGTACGCGTAAAGCg                                                                                                                                    <  2:50746/250‑1 (MQ=255)                                                                                                     gTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAg                                                                                                                                  <  2:337212/251‑1 (MQ=255)                                                                                                       ccAGTTCATAAGTCCAGGAGGATTTATCGTTCTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGga                                                                                                                                <  1:178072/251‑1 (MQ=255)                                                                                                       ccAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATCCGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTCCCGCCGCCGTACGCGTAAAGCGAgg                                                                                                                                 <  2:77841/250‑1 (MQ=255)                                                                                                        cAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgat                                                                                                                               >  1:256840/1‑251 (MQ=255)                                                                                                           gtcaaaaGTCCAGGAGGATTTATCGTACTTGCGTTTTGTCGCGTCGTCGCCGGTGAGCAGGCCATCTTCAACGCCGTAATCCTCACGCGCGTTCAGTCTGGCGTTTGTATAGGCTTCGGTGTATGCGCGGGTTAATTTTTCATTGTTCAGCAGGTTCAGCAATACGCCTGACAGGAAAGCTATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTCCCGCCGCCGTACGCGTAAATCGAGgatcga                                                                                                                            <  1:41588/246‑1 (MQ=255)                                                                                                              aTAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgatcgatc                                                                                                                          >  1:352055/1‑250 (MQ=255)                                                                                                                    ccAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATc                                                                                                                    >  2:319936/1‑250 (MQ=255)                                                                                                                     cAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGga                                                                                                                                >  1:163581/1‑237 (MQ=255)                                                                                                                     cAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTGCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCAGAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGga                                                                                                                                <  2:163581/237‑1 (MQ=255)                                                                                                                           gaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcgc                                                                                                            <  1:343489/251‑1 (MQ=255)                                                                                                                            aTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcgcg                                                                                                           >  2:60479/1‑251 (MQ=255)                                                                                                                               tATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCg                                                                                                        >  1:304132/1‑251 (MQ=255)                                                                                                                                             ttttCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTgg                                                                                           <  1:232086/250‑1 (MQ=255)                                                                                                                                                   gcttcgtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTACAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCAt                                                                                    <  2:36256/248‑1 (MQ=255)                                                                                                                                                             cGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCgg                                                                          >  1:179987/1‑251 (MQ=255)                                                                                                                                                                 gAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAAcc                                                                      <  1:183013/251‑1 (MQ=255)                                                                                                                                                                    cAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCg                                                                    >  1:181142/1‑250 (MQ=255)                                                                                                                                                                      ggCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTAGAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGGCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGAcc                                                                 <  1:41086/251‑1 (MQ=255)                                                                                                                                                                        ccATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCgg                                                               >  1:200975/1‑251 (MQ=255)                                                                                                                                                                             ttCAAAGCCGTAACCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATCCGCCTGACAGGAAAGCAATGTCAGTACCGGACCGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGa                                                           <  2:322431/250‑1 (MQ=255)                                                                                                                                                                                                cgcacgATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGTTTTAATGTTTCATTTTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCTCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGcc                                       <  1:20932/251‑1 (MQ=255)                                                                                                                                                                                                        cAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATcacca                               <  1:175154/251‑1 (MQ=255)                                                                                                                                                                                                                   tGGTATAGGTGTCGGTGTATGCGTGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCGGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTTAGCTACCTCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGTAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCTCCACGACGAGGtt                     <  1:77165/250‑1 (MQ=255)                                                                                                                                                                                                                     gTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGc                  <  1:154084/251‑1 (MQ=255)                                                                                                                                                                                                                       aaaGGCTTCGGTGTATGCGCGGTTTAATTTTTCAGTGTTCAGCAGGTACAGCATTACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAAAAGGGGCATAGTAGTCAGGTACCGCGGCCGTACGCGTAAAGCGAGGATCGATCACAAGCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGACATGCGGAACCCGGCCGTGTGAGCTTCAGCGGCGTGACTGCCCATCCCCAAGACGAGGtttgct                 <  2:7299/248‑5 (MQ=255)                                                                                                                                                                                                                                 gtgtATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTc      <  2:234331/251‑1 (MQ=255)                                                                                                                                                                                                                                 gtgtATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTc      >  1:232805/1‑251 (MQ=255)                                                                                                                                                                                                                                  cgatctCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGCTCGATCACAATCAGCTTCGCGCCGTTGTGATTTTTGGCTTCCATCGCCCAACGGAACCCGACCGGGTGAGCTTCAGCAGCGTTCACGCCCAACACCACGCCGAGGTTGGCGTTCATGATGTCg     >  1:265846/6‑251 (MQ=255)                                                                                                                                                                                                                                     aTTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGAcc  >  1:163834/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 CAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACGAGGTTGGCGTTCTTGATGTCGACC  >  NC_000913/4084694‑4085172

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