Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,084,903	T→G	N307T (AAC→ACC)	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,903	0	T	G	97.8%	151.9 / ‑6.3	46	N307T (AAC→ACC)	fdoG	formate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (20/25);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCAC  >  NC_000913/4084664‑4085144                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  aCGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCt                                                                                                                                                                                                                                        <  2:206869/251‑1 (MQ=255)       cTTTTAGCGTACCACAGAAGATTTCAACCACATCTGTCGTGTAGCCAGACACGTGCTCTTTCCGCAAGTTCCAAACGCAGAGCGGTTGTTGCAGCGAGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTAATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCCGGCTGGCGGTGGAATAGGCTTCgggg                                                                                                                                                                                                                                  >  2:424908/1‑249 (MQ=255)          ttGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTa                                                                                                                                                                                                                                <  2:108060/250‑1 (MQ=255)               gTACCACAGATGATTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTcgcg                                                                                                                                                                                                                          <  1:4169/251‑1 (MQ=255)                tACCACAGATGTTTTCAACCACATCTGGCGTGAAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGAAAGGCTTCGGGGTTTTCGCgg                                                                                                                                                                                                                         >  2:377611/1‑251 (MQ=255)                       gATGTTTTCAGCCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCTGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGGTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCTGACTTGCGGTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCTCGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATGAGCGGTTGAAtt                                                                                                                                                                                                                  <  2:264891/251‑1 (MQ=255)                       gATGTTTTCAACCACATCTGGAGTGTAGCGGGAAACGTGCTGTTTAAGCAAGTTCCACACGCAGCGCGGAGGGTGCAGCGTGGTATCGCGTTTGGCGAAACCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAACGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGGGTATTcaaggttgcatt                                                                                                                                                                                                                  >  2:305977/1‑240 (MQ=255)                                 aCCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGCATTTTTCATTGtt                                                                                                                                                                                                        >  1:229921/1‑251 (MQ=255)                                  ccACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGGTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTc                                                                                                                                                                                                       <  1:282990/251‑1 (MQ=255)                                  ccACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGGGTATTCGCGGTTGAATTTTTCATTGTTc                                                                                                                                                                                                       >  1:275133/1‑251 (MQ=255)                                         tGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGGGGTATAGGACTCGGTGTATTCGCGGGTGAATTTTTCATTGTTCAGCAGGt                                                                                                                                                                                                >  1:174431/1‑251 (MQ=255)                                              tgtAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGGGGTATAGGCTTCGGGGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTAc                                                                                                                                                                                              >  2:382344/1‑248 (MQ=255)                                                  gCGGGAAACGTGCTGTGTCATAAAGTTCCACAGGCAGCGCGGATGTGGCAGCGTGGTATCGTGTTTGGCGAGGCCTTTTTCGGCCAGTACATAAGTCCAGGAGGATTTATCGTACGTGCGTTTGTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCTCGGGCAGGGTGGCGGGGGTATAGGCTTCGGTGTATTCGCGGTTGAAATTTTCATTGTTCAGCAGGTACAGCAATa                                                                                                                                                                                       <  1:50139/251‑1 (MQ=255)                                                   cGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGAGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTCTAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCGGCAGGTACAGCAATa                                                                                                                                                                                       >  1:47088/1‑250 (MQ=255)                                                   cGGGAAACGTGCTGTCTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATAc                                                                                                                                                                                      <  1:130444/251‑1 (MQ=255)                                                     ggAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGGGTATTCGCGGTTGAATTTTTCATTGTTTAGCAGGTACAGCAATACGc                                                                                                                                                                                    >  2:297838/1‑251 (MQ=255)                                                       aaaCGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGGGGTATAGGCTTCGGTGGATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGccc                                                                                                                                                                                  >  2:138840/1‑250 (MQ=255)                                                       aaaCGTGCTGTTTCAGAAAGTGCCACCCGCAGCGCGGATGTTGCAGCGTGGTAGCGCGTTTGGCGAAGCAGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCt                                                                                                                                                                                  <  1:293372/251‑1 (MQ=255)                                                                gTTTCAGCTAGTTCCACACGCTGCGCGGATGTTGAAGCGTGGTATCTCGTTTGGCGAAGACGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGaaa                                                                                                                                                                         <  2:49555/251‑1 (MQ=255)                                                                  ttCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGc                                                                                                                                                                       >  1:388171/1‑251 (MQ=255)                                                                    cAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCa                                                                                                                                                                      >  1:300484/1‑250 (MQ=255)                                                                       cAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGt                                                                                                                                                                  <  1:224068/251‑1 (MQ=255)                                                                               acacGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCg                                                                                                                                                          <  1:242669/251‑1 (MQ=255)                                                                                              gTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGcc                                                                                                                                           >  2:209429/1‑250 (MQ=255)                                                                                                     cGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATagtagt                                                                                                                                    <  1:146709/251‑1 (MQ=255)                                                                                                      gtggtagagcgTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTGCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGTCGCATagtagt                                                                                                                                    <  2:74110/242‑1 (MQ=255)                                                                                                             cgcgTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAcc                                                                                                                            <  1:34335/251‑1 (MQ=255)                                                                                                                    ggCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGTGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGACAAGGGGCAAAGTAGACAGCTAccgccgccg                                                                                                                     >  2:363901/1‑251 (MQ=255)                                                                                                                     gCGAAGCCTTTTTCGTCCAGTTCATAAGACCAGGAGGAGTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGGGAACAGGCCATCTGCAAAGCCGTAATCCTCACCCCCGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGGACCGGAACGAATAGGGGCATAGAAGTCAGCTACCGCCGCCGt                                                                                                                    <  2:388171/251‑1 (MQ=255)                                                                                                                                   gTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAg                                                                                                      <  1:280826/251‑1 (MQ=255)                                                                                                                                     ccAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGga                                                                                                    <  2:76221/251‑1 (MQ=255)                                                                                                                                      cAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgat                                                                                                   >  1:250166/1‑251 (MQ=255)                                                                                                                                       aGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgatc                                                                                                  <  1:198157/251‑1 (MQ=255)                                                                                                                                         ttCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACCGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCAACCGCCGCCGTACGCGTAAAGCGAGgatcg                                                                                                 >  1:238606/1‑250 (MQ=255)                                                                                                                                                  ccAGGAGGATTTATCGTACTGGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTGCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTCCAGCAGGTACAGCAATACGCCTTACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCa                                                                                       <  2:57860/251‑1 (MQ=255)                                                                                                                                                    aggaggATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGCGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACAGGAACGAATAGGGGCATAGTAGTCAGATACAGCCGCCGAACGCGTAAAGCGAGGGACGATAAAAATCAGc                                                                                     >  1:21426/1‑251 (MQ=255)                                                                                                                                                         gaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcgc                                                                                >  1:24848/1‑251 (MQ=255)                                                                                                                                                             tATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCAAAATCAGCTTCGCGCCg                                                                            >  2:385680/1‑251 (MQ=255)                                                                                                                                                             tATCGTACTGGCGTTTTGCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAAGCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTGCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGTATAGGAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGcc                                                                             <  2:168925/250‑1 (MQ=255)                                                                                                                                                               ccttACTTGCGTTTTTCCGCGCCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACTATCAGGCTGGCGGTGTTATAGGCTTCGGTGTATGCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGtt                                                                          <  1:80118/248‑1 (MQ=255)                                                                                                                                                                                            ggTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTTTAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCCTAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGAACGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCggc                                             >  2:297899/1‑250 (MQ=255)                                                                                                                                                                                                              aaaGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCtt                           <  2:168700/251‑1 (MQ=255)                                                                                                                                                                                                                              cgcacgATCAGTCAGGCGGTTGTAGAGGCGTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTTACAGAAAAGCAATGTCAGTACCTGAACGAATAGGGGCATAGTGTTCAGCTTCCGCCGCCGTACTCGTAAAGCGAGTATCGATCACAATCAGCTTCTTGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGcc           <  1:181778/251‑1 (MQ=255)                                                                                                                                                                                                                                cacgATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCCCCGTTGTGCATTTTGGCTTCCATCGCCCCGCGGAACCCGACCGGGTGCGCTTCAGCGGCGTTACCGcccc         >  2:331033/1‑250 (MQ=255)                                                                                                                                                                                                                                      cAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGTTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATcacca   <  1:305977/251‑1 (MQ=255)                                                                                                                                                                                                                                       aGGCTGGCGGGGGTATAGGCTTCGGTGTATTCGCGGGTGAATTTTTCATTATTCAGCAGGTACACCAATACTCCTGACATGAAAGCACTGTCATTACCGGAACGAATAGGGGCATAGTAGGCAGCTACCGACGCCGTACGCGTACAGCGAGTCTCGATCACAATCAGCTTCGCGCCGTTGTGAATTGTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCTGCGTTACCGCCCATcaccac  <  2:171222/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  ACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCAC  >  NC_000913/4084664‑4085144

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