Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,360,645	T→A	8.9%	intergenic (+29/‑531)	gltD → / → gltF	glutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,360,645	0	T	A	8.9%	149.7 / 4.2	45	intergenic (+29/‑531)	gltD/gltF	glutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (18/23);  new base A (2/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GGCGTTTGGTTTCCGTCCACACAACATGGAATGGCTGGCAAAACACAGCGTCGAGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAA  >  NC_000913/3360388‑3360876                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        ggCGTTTGGTTTCCGTCCACACAACATGGAATGGCTGGCAAAACACAGCGTCGAGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAAtt                                                                                                                                                                                                                        >  2:259860/1‑275 (MQ=255)   cGTGTGGTTTCCGTCAACACAACATGGAATGGCTGGCAAAACCCAGCGTCGGGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCGGCAACCCGAAAATCTGTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGAACGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTa                                                                                                                                                                                                                       <  4:36608/274‑1 (MQ=255)                tCCAAACAACAGGGAATGGCTGGCAAAAAACAGCGTCGAGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTgcg                                                                                                                                                                                                          <  1:86773/274‑1 (MQ=255)                    cacaACATGGAAGGGCTGGCAAAACACAGCTTCGGGCTGGATTCACAAGGCCGTATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAAGCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCgtt                                                                                                                                                                                                      <  2:341171/274‑1 (MQ=255)                                   ctggcAAAACACAGCGTCGAGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGCCGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGAAGAAGTTCGCGTTGTTTTATGGCacaagc                                                                                                                                                                                      >  4:519706/1‑270 (MQ=255)                                                  gtctggCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCGTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTg                                                                                                                                                                        <  2:332436/269‑1 (MQ=255)                                                     gAGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGAcc                                                                                                                                                                     <  1:34159/274‑1 (MQ=255)                                                           gATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGttt                                                                                                                                                                         >  2:337297/1‑264 (MQ=255)                                                               cacaAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGAcc                                                                                                                                                                     <  4:337296/264‑1 (MQ=255)                                                                acaAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataat                                                                                                                                                          <  4:259859/274‑1 (MQ=255)                                                                                                                    aaCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCtc                                                                                                                                                                                                                                    <  8:781445/148‑2 (MQ=255)                                                                                                                    aaCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTg                                                                                                                                                                                                                                    <  6:813909/148‑1 (MQ=255)                                                                                                                                             atatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCaa                                                                                                                                                                                                                           <  8:89625/132‑1 (MQ=255)                                                                                                                                             atatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCaa                                                                                                                                                                                                                           >  7:89625/1‑132 (MQ=255)                                                                                                                                             atatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAg                                                                                                                                                                                                              >  6:923303/1‑145 (MQ=255)                                                                                                                                             atatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAg                                                                                                                                                                                                              <  5:923303/145‑1 (MQ=255)                                                                                                                                                               aTCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGt                                                          >  4:253415/1‑275 (MQ=255)                                                                                                                                                               aTCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACACAGTCTGGCGa                                                                                                                                                                                                                                >  6:940360/1‑109 (MQ=255)                                                                                                                                                               aTCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACACAGTCTGGCGa                                                                                                                                                                                                                                <  5:940360/109‑1 (MQ=255)                                                                                                                                                                      ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtgtg                                                                                                                                                                                                                                        <  7:291268/94‑1 (MQ=255)                                                                                                                                                                      ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtgtg                                                                                                                                                                                                                                        >  8:291268/1‑94 (MQ=255)                                                                                                                                                                      ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTgg                                                                                                                                                                                                                                   >  7:450405/1‑99 (MQ=255)                                                                                                                                                                      ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTgg                                                                                                                                                                                                                                   <  8:450405/99‑1 (MQ=255)                                                                                                                                                                         gACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACctc                                                 <  2:253416/274‑1 (MQ=255)                                                                                                                                                                          aCCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGa                                                                                                                                                                       >  2:436884/1‑155 (MQ=255)                                                                                                                                                                                     cGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCt                                                                                                                                                                                               >  5:129991/1‑120 (MQ=255)                                                                                                                                                                                     cGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCt                                                                                                                                                                                               <  6:129991/120‑1 (MQ=255)                                                                                                                                                                                         ggTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCta                                                                                                                                                              >  2:365869/1‑149 (MQ=255)                                                                                                                                                                                             gTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataat                                                                                                                                                          <  4:365868/149‑1 (MQ=255)                                                                                                                                                                                             gTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCaa                             <  2:519707/274‑1 (MQ=255)                                                                                                                                                                                              tAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGa                                                                                                                                                                       <  5:640810/135‑1 (MQ=255)                                                                                                                                                                                              tAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGa                                                                                                                                                                       >  6:640810/1‑135 (MQ=255)                                                                                                                                                                                                                tatGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCgttg                                                                                                                                                                                                     <  5:212446/87‑1 (MQ=255)                                                                                                                                                                                                                tatGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCgttg                                                                                                                                                                                                     >  6:212446/1‑87 (MQ=255)                                                                                                                                                                                                                     aCTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAAtttt                                                                                                                                              >  5:486211/1‑137 (MQ=255)                                                                                                                                                                                                                     aCTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAAtttt                                                                                                                                              <  6:486211/137‑1 (MQ=255)                                                                                                                                                                                                                       tggctggAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTaa  >  2:161052/1‑275 (MQ=255)                                                                                                                                                                                                                                  ttAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataat                                                                                                                                                          <  8:386492/112‑1 (MQ=255)                                                                                                                                                                                                                                  ttAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataat                                                                                                                                                          >  7:386492/1‑112 (MQ=255)                                                                                                                                                                                                                                                    aaTTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGCCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTTTTCACCAAGTCTGATGAGTGAAGCACAGAAATCAATTTCTCATTTACAATCCGTCGCa                                                                   >  3:11206/1‑181 (MQ=255)                                                                                                                                                                                                                                                     aTTCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTg                                                                                                                                                                        >  7:507159/1‑79 (MQ=38)                                                                                                                                                                                                                                                     aTTCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTg                                                                                                                                                                        <  8:507159/79‑1 (MQ=38)                                                                                                                                                                                                                                                      ttCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTc                                                                                                                                                                                                >  6:110761/1‑54 (MQ=255)                                                                                                                                                                                                                                                      ttCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTc                                                                                                                                                                                                <  5:110761/54‑1 (MQ=255)                                                                                                                                                                                                                                                        cGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGAAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTa                                                               <  1:11206/181‑1 (MQ=255)                                                                                                                                                                                                                                                              aCTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtga                                                                                                               <  6:301873/127‑1 (MQ=255)                                                                                                                                                                                                                                                              aCTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtga                                                                                                               >  5:301873/1‑127 (MQ=255)                                                                                                                                                                                                                                                                 gtgtCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAAttt                                                                                                                                               >  7:266258/1‑92 (MQ=255)                                                                                                                                                                                                                                                                 gtgtCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAAttt                                                                                                                                               <  8:266258/92‑1 (MQ=255)                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        GGCGTTTGGTTTCCGTCCACACAACATGGAATGGCTGGCAAAACACAGCGTCGAGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAA  >  NC_000913/3360388‑3360876

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