Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,496,751	G→A	W67* (TGG→TGA)	nirD →	nitrite reductase (NADH) small subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,496,751	0	G	A	97.7%	149.3 / NA	43	W67* (TGG→TGA)	nirD	nitrite reductase (NADH) small subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/27);  minor base C (1/0);  total (16/27)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGGGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAA  >  NC_000913/3496510‑3496992                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  ccGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGc                                                                                                                                                                                                                                          <  1:39818/251‑1 (MQ=255)                ccAGTAACGCTGGTGGAGGACAACGCATGAGACAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTGCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAg                                                                                                                                                                                                                           <  1:151/251‑1 (MQ=255)                         ctGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGc                                                                                                                                                                                                                  >  2:97820/1‑251 (MQ=255)                            gtggAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTa                                                                                                                                                                                                               <  2:57229/251‑1 (MQ=255)                                     aaCGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGc                                                                                                                                                                                                      >  2:227392/1‑251 (MQ=255)                                     aaCGCATGAGCCAGTGGAAAGACATCGGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTGTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGGGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGATCTGTGAGTCGCCAGCCCGCTGAAGAAACAGCGTTTTCGCTTAAGCGACGGc                                                                                                                                                                                                      <  2:84656/251‑1 (MQ=255)                                       cGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCAATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCtt                                                                                                                                                                                                    >  2:180399/1‑251 (MQ=255)                                             aGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCAt                                                                                                                                                                                              >  2:83128/1‑251 (MQ=255)                                                     gAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACg                                                                                                                                                                                      <  1:120573/251‑1 (MQ=255)                                                      aaaGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTc                                                                                                                                                                                                                    <  2:29467/220‑1 (MQ=255)                                                      aaaGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTc                                                                                                                                                                                                                    >  1:29467/1‑220 (MQ=255)                                                         gACAACTGCAAACTCGATGACATGCTGCCTGAAACCGGCGTATGCGCGCTGTTAGGGGACGAGCCGGTCGCGATTTTCCGCCCGTCTCACAGCGATCAGGAGTTTGCGCACAGCAACCTCGACCCGTTCTTCGCGTCCAGCGTGCTGTCACGCGGACTGATTGCGTATCACCAGGTCGATCTGTGCGGCGTCAGCCCGCGGCACAAACAGCGTCTTCGCTTAAGTGACggtgtgttcaagggaggcgaac                                                                                                                                                                                   >  2:82330/1‑230 (MQ=255)                                                                gCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGTCTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTc                                                                                                                                                                            <  2:19552/250‑1 (MQ=255)                                                                     aTCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCaa                                                                                                                                                                      <  1:227392/251‑1 (MQ=255)                                                                         aTGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGCCGAACAGTTTTCCGTCAAACAt                                                                                                                                                                  >  2:179236/1‑251 (MQ=255)                                                                                           ccGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGaaa                                                                                                                                                <  1:180009/251‑1 (MQ=255)                                                                                             ggCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTACGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGa                                                                                                                                              <  2:169338/251‑1 (MQ=255)                                                                                                    gcgcgcTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCgt                                                                                                                                        <  1:176144/250‑1 (MQ=255)                                                                                                         cTGTTAGGTGACGAGAATGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGGGTTTGCGATCCGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGATTTCCGACAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGc                                                                                                                                   >  2:147727/1‑250 (MQ=255)                                                                                                                              gcgATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTaa                                                                                                             <  2:216630/251‑1 (MQ=255)                                                                                                                                 aTTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACg                                                                                                                                                              <  1:209305/199‑1 (MQ=255)                                                                                                                                 aTTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACg                                                                                                                                                              >  2:209305/1‑199 (MQ=255)                                                                                                                                       cGCCCGTATCACGGCGAGCAGGTGTTTGCGATCAGCAACATCGACCCGGTCTTCGATTCCAGCGGGCTGGTACGCGGACTGATTGCGGGACACGAGGGCGAGCTGTGAGTCGCGAGCCCGCTGAAAAAACATCGTTTTCGCGTAAGCGACGGCTTGTGCAAGGATGACCAACAGTTGTCCGTCAAACAGTACGAAGCGCGAGTGAAAGACGGCATGGTGCAGCTGCGCGGTTAATGGTTTAACGTGAGgcg                                                                                                    <  2:104852/251‑1 (MQ=255)                                                                                                                                       cGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGgcg                                                                                                    <  1:207338/251‑1 (MQ=255)                                                                                                                                                 acaGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTcc                                                                                          >  2:105623/1‑251 (MQ=255)                                                                                                                                                      gATCAGGAGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTCAGCGACGGCTTGTGCATGGCAGACGAACAGTTTTCCGTCAAACATTCCGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGCGGCGCAATGCCTCCCCttt                                                                                     >  2:211611/1‑251 (MQ=255)                                                                                                                                                            gtgtTTGCGATCAGCAACATAGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCAt                                                                               >  2:37445/1‑251 (MQ=255)                                                                                                                                                                ttGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCGTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGTCGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGt                                                                            <  1:105404/250‑1 (MQ=255)                                                                                                                                                                                agACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCgg                                                           <  1:37445/250‑1 (MQ=255)                                                                                                                                                                                     cGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTatat                                                      <  2:99919/251‑1 (MQ=255)                                                                                                                                                                                     cGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTatat                                                      <  2:188279/251‑1 (MQ=255)                                                                                                                                                                                      gTTCTTCGAGTCCCGCGTGCTCTCACGCGGACTGCTTGCGGAACACCACGCCGAGCTCTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGTCGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACCTTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATCCCTCCCCTTTTTGCATGGTCCTGTAATACTCTTCGTtatctt                                                     >  2:237306/1‑248 (MQ=255)                                                                                                                                                                                           tcGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCa                                                  <  1:129448/249‑1 (MQ=255)                                                                                                                                                                                            cGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACTGGAGGCGCAATGCCTCCCCTTTTTTCAGGGTCCTGTAATAGTCGCCGGTATAGCGCAGGa                                               <  2:24316/251‑1 (MQ=255)                                                                                                                                                                                             gAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACTGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGAc                                              <  2:181763/251‑1 (MQ=255)                                                                                                                                                                                                 ccatcGTGCTGTCACGCGGTCTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGTTGCGCGTTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACAttt                                          <  1:224886/247‑1 (MQ=255)                                                                                                                                                                                                     cGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTa                                      <  1:85735/251‑1 (MQ=255)                                                                                                                                                                                                             cACGCGGACTGATTGCGGAACGCCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACttttt                              <  2:194400/251‑1 (MQ=255)                                                                                                                                                                                                             cACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACtttt                               <  1:130909/250‑1 (MQ=255)                                                                                                                                                                                                              aCGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTg                             >  1:4282/1‑251 (MQ=255)                                                                                                                                                                                                                              gAACACCTGAGCGAGCTGTGAGTCGACAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGAAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGtttt             >  2:245158/1‑251 (MQ=255)                                                                                                                                                                                                                                       ggCGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTAtttttt     >  1:178588/1‑250 (MQ=255)                                                                                                                                                                                                                                         cGAGCTGTGAGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCCGCGCGGTTAATGTTTTAACGGTAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAAGCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTGTTATTTTTTGTGTTTATGTTTTaaa  <  1:211611/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  CCGTATGAACGTATCCCAGTAACTCTGGTGGAGGACAACGCATGAGCCAGTGGAAAGACATCTGCAAAATCGATGACATCCTGCCTGAAACCGGCGTCTGCGCGCTGTTAGGTGACGAGCAGGTCGCGATTTTCCGCCCGTATCACAGCGATCAGGTGTTTGCGATCAGCAACATCGACCCGTTCTTCGAGTCCAGCGTGCTGTCACGCGGACTGATTGCGGAACACCAGGGCGAGCTGTGGGTCGCCAGCCCGCTGAAAAAACAGCGTTTTCGCTTAAGCGACGGCTTGTGCATGGAAGACGAACAGTTTTCCGTCAAACATTACGAAGCGCGAGTGAAAGACGGCGTGGTGCAGCTGCGCGGTTAATGTTTTAACGGGAGGCGCAATGCCTCCCCTTTTTGCATGGTCCTGTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTTTATTTTTTGTTTTTATTTTTTAAA  >  NC_000913/3496510‑3496992

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