Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	22,048	G→A	L214L (CTG→CTA)	ribF →	bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	22,048	0	G	A	91.7%	93.2 / ‑4.1	38	L214L (CTG→CTA)	ribF	bifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (9/24);  total (10/28)
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

AAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATT  >  NC_000913/21807‑22294                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        aaGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTCTAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTc                                                                                                                                                                                                                                                  >  1:192566/1‑249 (MQ=255)    gCGATTTCTTGTTATTACAGAAAGCTGGTATTGAATACGTCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGtt                                                                                                                                                                                                                                             <  1:343947/251‑1 (MQ=255)                    aCAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGGGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCTCCCGTGTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAggg                                                                                                                                                                                                                              <  1:193991/250‑1 (MQ=255)                       gaaatcttggaTGGAATACGGCTTCGATATCTCCTGTTCGCAATCTTTATGCGAAGGGGGCGTGCGTATCCGTAGCGCCGCTGTGCGACAGGCGCGGGCGGATGACAATCTGGCGCTGGCAGAAAGTTTACTGGGGCTCCCGTTTACCATCGCCGGGCGTGCAGTCCACGGTGATGAATTAGGGCGCACTATAGGTGTCCCGACGGCGAATGGTCCGGTAAGCCGTCGTGTTTCCCCGGAGAAAggggctt                                                                                                                                                                                                                          <  2:44272/241‑4 (MQ=255)                          aGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGGGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGCAAGGGGTTTATg                                                                                                                                                                                                                       >  1:12320/1‑251 (MQ=255)                                cATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGt                                                                                                                                                                                                                   >  1:221959/1‑249 (MQ=255)                                           gCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGc                                                                                                                                                                                                      >  1:341982/1‑251 (MQ=255)                                              tCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGcc                                                                                                                                                                                                     <  1:25073/249‑1 (MQ=255)                                                 atatCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTCGAAGTGCTGGGCCttggg                                                                                                                                                                                                >  2:108154/1‑247 (MQ=255)                                                                     gtttGCGAAGGTGGCGTGCGTATCAGCAGCCCCGCTTTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTCCCGCTTCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCg                                                                                                                                                                            <  2:303349/249‑1 (MQ=255)                                                                         gCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCa                                                                                                                                                                       <  1:206823/251‑1 (MQ=255)                                                                                    gTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAac                                                                                                                                                            <  1:172292/251‑1 (MQ=255)                                                                                       cgtaTCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTACGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACg                                                                                                                                                         >  1:240059/4‑251 (MQ=255)                                                                                          aTCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGccc                                                                                                                                                      <  1:271038/251‑1 (MQ=255)                                                                                                   aCCGCTGTGCGTCAGGCGCTGGCGGATGACAATATGGCTCTGGCAGAATGTTTACTGGGGCACCCGTTGGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGc                                                                                                                                             <  1:298135/251‑1 (MQ=255)                                                                                                            cggcaggcgatggcagATGACAATCTGGCTCTGGCAGAAAGTTTATTGGGGTCTCCGGTTGCCATCTCCGGGCGTGTATTCCTCGGTGATGAATTAGGGCGCACTAGAGGTTTCCCGACGGCGAAGGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTc                                                                                                                                     <  2:285546/235‑1 (MQ=255)                                                                                                             gTCAGGCGCTGGCGGATGTCAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGTGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAATCATCTGAACACGCCCAACGGTTGCCGGTATTCGc                                                                                                                                   <  1:122254/251‑1 (MQ=255)                                                                                                                       ggCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTg                                                                                                                          <  1:234497/249‑1 (MQ=255)                                                                                                                        gCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGccc                                                                                                                                                      >  2:274079/1‑221 (MQ=255)                                                                                                                        gCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGccc                                                                                                                                                      <  1:274079/221‑1 (MQ=255)                                                                                                                                              cAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGAAGCTGGAAGTGCATTTGTTCGATGTTGc                                                                                                   >  2:64411/1‑250 (MQ=255)                                                                                                                                                gaaagTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAAt                                                                                                <  1:87546/248‑1 (MQ=255)                                                                                                                                                     tgtacgggggCACCTGTTTGCCCTCGCCGGGGGTGTAGTCCACGGTGATGAATGAGGTCGTGGTATAGGGTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGGTTCCCCGGTGAAAGGGGTTTATGCGGTAGAATTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGGGGCAAACATCGGAACACGCCCAACGGTTGCCGGTAGTCGCCAGCAGCTTGAAGTGCATTTGTTAGATGTTGCAATGGAc                                                                                            <  2:378035/244‑1 (MQ=255)                                                                                                                                                           ggggCACCCGTGTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTTAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATTGATCTTTACg                                                                                     <  2:240059/251‑1 (MQ=255)                                                                                                                                                                        gCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACa                                                                        <  1:271258/251‑1 (MQ=255)                                                                                                                                                                               ccGGGTGTGTATACCACTGTGTTGAAGTATGGCGCACTATTGTTTTCCCGATGTCGAATGTACAGCTACTGCGTCAGGTTTCCCCGGTGCAAGGGGGTTATTCGGTAGAAGTTCTGGTCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCTGAACACGCCCAACGGTGGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATAccattagt                                                                  <  2:172048/250‑8 (MQ=255)                                                                                                                                                                                            ccACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCACGGTGAAAGGGGGTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTaaaaaa                                                     <  1:208674/250‑1 (MQ=255)                                                                                                                                                                                             cACGGTGATGAATTAGGGCGCACGATAGGTTTCCCGACGGCGAATGGACCGTTACGCCGTTAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATTCAAGTAGTGCGGCGTAAAAAAAt                                                   <  2:341982/251‑1 (MQ=255)                                                                                                                                                                                                          aaGGGCGCTCTATAGGTTTCCCGACTGCGAATGTTCCTCTCAGCCGTCTGGGTGCCCCGGTGCATGGGGTTTATGCTGTAGAAGTGCGTGGCCG‑CGTTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGTTTTCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTCGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGagcag                                      <  2:326115/250‑1 (MQ=255)                                                                                                                                                                                                              gcgcACTATAGGTTTCCCGACGGAGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGAt                                  <  2:221959/251‑1 (MQ=255)                                                                                                                                                                                                                gcaatataGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGAttt                                <  1:303053/247‑1 (MQ=255)                                                                                                                                                                                                                 cACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTg                               >  1:151963/1‑251 (MQ=255)                                                                                                                                                                                                                     ataGGTTTCCCGACGTCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGTAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGGCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGt                            <  1:234338/250‑1 (MQ=255)                                                                                                                                                                                                                         gTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGtt                                                                                                     >  1:225561/1‑173 (MQ=255)                                                                                                                                                                                                                         gTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACAGCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGtt                                                                                                     <  2:225561/173‑1 (MQ=255)                                                                                                                                                                                                                                gACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAAc                <  1:334458/251‑1 (MQ=255)                                                                                                                                                                                                                                       aaTGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAgc         <  1:64411/251‑1 (MQ=255)                                                                                                                                                                                                                                               tgtGCGCCGTCAGGTGTCCCCGGTGAAAGGGGTTTATGCGGGAGAAGCG‑TGGGGCTGCGGTGAAATGGCGGTATCCGGCGTGGCAAACAGTGGAACACGCCCAAAGGTTGACGTTATTCGCCAGCAGCTGGAAGTTCATTTGTTAGATGTTGCAATGGACCTGTACGGTCGGTATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGGCGCTGGACGAACTGAAAGCGCAGAtt  <  2:141388/248‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        AAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCT‑CGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATT  >  NC_000913/21807‑22294

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