Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,302,355	G→T	T133T (ACC→ACA)	napA ←	periplasmic nitrate reductase subunit NapA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,302,355	0	G	T	95.6%	126.7 / ‑2.2	45	T133T (ACC→ACA)	napA	periplasmic nitrate reductase subunit NapA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (23/20);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

TGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAA  >  NC_000913/2302108‑2302590                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            tGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGtt                                                                                                                                                                                                                                          >  1:289552/1‑251 (MQ=255)  gTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGtt                                                                                                                                                                                                                                          <  1:153238/250‑1 (MQ=255)  gTTTGCGCCCCACAGCACAAACGCGGCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGttt                                                                                                                                                                                                                                         <  1:356732/251‑1 (MQ=255)           ccACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTATTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTtct                                                                                                                                                                                                                                <  2:25360/251‑1 (MQ=255)             acaGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTctct                                                                                                                                                                                                                              >  2:336389/1‑251 (MQ=255)                 cacaAACGCGTCAGCCTGCTCGCTGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTAACGCGCGTTCGGGTCGCTGTTGTTCGAACGGAAGCCCGCTTTAACCACCTTGGACGCGGCATAACCTTCCCAGATAGTCCCCTGACCAGAACCGAACATACCCACCAATTCCGGCCCTTTTTATTCCAGGACTGTATTGAATTTTTCTTccc                                                                                                                                                                                                                          >  2:2761/1‑250 (MQ=255)                 cacaAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCa                                                                                                                                                                                                                          <  1:349273/251‑1 (MQ=255)                  acaAACGCGTCAGCCTTCTCTATGTCGTTATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACTCGCTTTCGGGTCGATGTTGTTCGAACTGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCTAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCAt                                                                                                                                                                                                                         <  2:426556/251‑1 (MQ=255)                                 tGCTCGAGGACGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGAGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGGCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGATGATTAGAATTTCTCTTCCAACACATCGCAGGCCTg                                                                                                                                                                                                          >  2:375013/1‑251 (MQ=255)                                         gtcgtcATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAgg                                                                                                                                                                                                  <  1:298389/251‑1 (MQ=255)                                             tcATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGGTCCCAGGTGAt                                                                                                                                                                                              >  2:76054/1‑251 (MQ=255)                                                   cagcCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGt                                                                                                                                                                                        <  1:257457/251‑1 (MQ=255)                                                                            aGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATACCCTTCCCAGATAGTCCACTGACCAGAACCGAACATAACGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGTTCCCAGGGGATTAGGGTAAATTCGCTTTCTTTGTCATAttt                                                                                                                                                                >  2:275839/1‑250 (MQ=255)                                                                                          ccAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTACGGGCTGTTTTGAATTTCTCTTCCAGCACATCGAAGGCCTGATCCCAGGGGATTGGGGAAAATTCGCCTTCTTTGTCTGATTTACCGtttttattagc                                                                                                                                                 >  2:71907/1‑245 (MQ=255)                                                                                                actGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACgcagca                                                                                                                                            <  1:336389/250‑1 (MQ=255)                                                                                                     cGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCCGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCt                                                                                                                                      >  2:346704/1‑251 (MQ=255)                                                                                                      gACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTg                                                                                                                                     <  1:385610/251‑1 (MQ=255)                                                                                                           cATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGATCATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCa                                                                                                                                <  2:324367/251‑1 (MQ=255)                                                                                                            aTACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCa                                                                                                                                >  1:37181/1‑250 (MQ=255)                                                                                                              aCAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCTGCTGCGTCAAAc                                                                                                                             <  2:450572/251‑1 (MQ=255)                                                                                                                             cGGGACGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATCCGTAGCATCGGCTGCGTCAAACGGTATTTAccttcc                                                                                                               >  2:341413/1‑246 (MQ=255)                                                                                                                                      gttgttCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTgg                                                                                                     >  1:445915/1‑251 (MQ=255)                                                                                                                                      gttgttCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTgg                                                                                                     <  1:80140/251‑1 (MQ=255)                                                                                                                                      gttgttCGAACGGAAGCCCGCGTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTgg                                                                                                     <  1:341413/251‑1 (MQ=255)                                                                                                                                       ttgttCGAACGGAAGCCCGCTTTCACCAGCTTGGACGCGGCCTAACCTTCCCAGATCGTCCCCTGACCCGAACCGCCCATACCGATCGATTCCAGCCCTTTTTCTTTCAGGTCTGTTTTGAATTTCTCTTCCATCACATCGTAGGCCTGATACCAGGTGGTTGGGTTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACTCCGCAGCGGCTGCGTCAAACGGTCTTTACCGAACATGATTTTgtg                                                                                                    >  2:157679/1‑249 (MQ=255)                                                                                                                                         gttCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATTATTTTGGGCa                                                                                                  <  1:226065/251‑1 (MQ=255)                                                                                                                                         gttCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGACTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGc                                                                                                   >  1:387463/1‑250 (MQ=255)                                                                                                                                              aaCGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTTAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAACAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGaa                                                                                              >  1:137456/1‑250 (MQ=255)                                                                                                                                               aCGGAAGCCCGCTTTAAACAGCATGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGCGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAAt                                                                                            >  1:399048/1‑251 (MQ=255)                                                                                                                                                            ttAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGTAAATAGCCCTTAATGc                                                                                <  1:346704/250‑1 (MQ=255)                                                                                                                                                            ttAAACAGCTTGGACGCGGCATAACCTCCCCAGATAGTCCACTGACCAGAACCGAACATACCTATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTTAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAACAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATg                                                                                 <  2:137456/249‑1 (MQ=255)                                                                                                                                                              aaaCAGCTTGGACGCGGCATAACCTTCCCAGATTGTTCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTGTTTCTTTCAGGGCTGGTTTGAATTTCGCTTCCATCACATCGAAGGCCTGATCCCAGGGGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACAGTACATGATGTTGGGCAGGAAATAGCCCTTAATGCAGt                                                                             <  2:388472/251‑1 (MQ=255)                                                                                                                                                                 cAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCCTATTTACCGTTTTTCATACGCAGCCGCGGCTGCGTCCAACGGTCTTTACCGTACCTGATTTTAGGCAGGAAATAGCCCTTAATGCAGTTc                                                                           >  2:347108/1‑250 (MQ=255)                                                                                                                                                                       ggACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGCTCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGCAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTAAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCAc                                                                    >  1:284169/1‑251 (MQ=255)                                                                                                                                                                       ggACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCAc                                                                    <  1:76054/251‑1 (MQ=255)                                                                                                                                                                             ggCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACgcagca                                                                                                                                            <  2:334736/173‑1 (MQ=255)                                                                                                                                                                             ggCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACgcagca                                                                                                                                            >  1:334736/1‑173 (MQ=255)                                                                                                                                                                                      ttCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGt                                                     >  1:75795/1‑251 (MQ=255)                                                                                                                                                                                                 tCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGCAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGccc                                          >  1:392857/1‑251 (MQ=255)                                                                                                                                                                                                  ccACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCt                                         <  2:385786/251‑1 (MQ=255)                                                                                                                                                                                                           aGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGccc                                >  2:212605/1‑250 (MQ=255)                                                                                                                                                                                                                 gAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGACGGCCTGATCCCCGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGACCGGGTCGCCCTGACAGGccacaccac                          >  2:317944/1‑246 (MQ=255)                                                                                                                                                                                                                               atTCAGGCCCTTTTTCTTTAAGGGCTGTTTTGAATTTCTCTTCCCTAACATCGAAGGCCTGGTCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGCAATCGCCCTTAATGCAGTTCAGACCACGGTTACCCGGTGCGTCCGGGTCGCACTGAAAGGCCACCACCCGTCCCTGCTGGGtt            >  2:311427/1‑251 (MQ=255)                                                                                                                                                                                                                                    ggCCCTTTTTCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCTGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGAc       >  1:194399/1‑251 (MQ=255)                                                                                                                                                                                                                                          ttttCTTTCAGGGCTGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATGCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGaa  <  1:459683/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            TGTTTGCGCCCCACAGCACAAACGCGTCAGCCTGCTCGATGTCGTCATAGCAGCCCATCGGCTCATCCATACCAAAGGTACGCATAAAGCCAACTACTGCCGACGCCATACAGTGACGCGCGTTCGGGTCGATGTTGTTCGAACGGAAGCCCGCTTTAAACAGCTTGGACGCGGCATAACCTTCCCAGATAGTCCACTGACCAGAACCGAACATACCGATCGATTCCGGCCCTTTTTCTTTCAGGGCGGTTTTGAATTTCTCTTCCATCACATCGAAGGCCTGATCCCAGGTGATTGGGGTAAATTCGCCTTCTTTGTCATATTTACCGTTTTTCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAA  >  NC_000913/2302108‑2302590

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