Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,609,208	(A)5→4	intergenic (‑212/‑9)	tusA ← / → yhhQ	sulfur transfer protein TusA/putative queuosine precursor transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,609,204	0	A	.	100.0%	173.6 / NA	41	intergenic (‑208/‑13)	tusA/yhhQ	sulfur transfer protein TusA/putative queuosine precursor transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (21/20);  total (21/20)

GAAAAAAGGGTAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATCAAAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCGTTATTTCATCT  >  NC_000913/3609124‑3609278                                                                                 |                                                                           gAAAAAAGGGTAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATCa                                                                            <  2:631663/81‑1 (MQ=255)  aaaaaaGGGTAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATCa                                                                            >  2:631737/1‑80 (MQ=255)   aaaaaGGGTAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATCaaa                                                                          <  2:631664/81‑1 (MQ=255)      aaGGGTAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCTCTCACCCCAATGGTTAATCaaaagg                                                                       >  2:631732/1‑79 (MQ=255)       aGGGTAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATCaaaaggt                                                                      <  1:631667/81‑4 (MQ=255)       aGGGTAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATCaaaaggt                                                                      <  2:631667/81‑4 (MQ=255)       aGGGTAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATCaaaaggt                                                                      <  1:631666/81‑4 (MQ=255)       aGGGTAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATCaaaagg                                                                       >  1:631748/1‑78 (MQ=255)       aGGGTAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATCaaaagg                                                                       >  2:631713/1‑78 (MQ=255)          gTAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTAc                                                                   >  1:631713/1‑80 (MQ=255)           tAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACaa                                                                 >  1:631719/1‑81 (MQ=255)             aGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAat                                                                <  2:631668/80‑1 (MQ=255)                  gTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAAc                                                          >  2:631720/1‑81 (MQ=255)                  gTTATTACGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAAc                                                          >  1:631722/1‑81 (MQ=255)                     attatGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGt                                                        <  2:631669/80‑1 (MQ=255)                      ttatGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGtt                                                       <  2:631670/80‑1 (MQ=255)                       aatGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGtttt                                                     >  2:631702/2‑81 (MQ=255)                         tGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTtc                                                    <  2:631671/80‑1 (MQ=255)                           ttAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTtctc                                                  <  2:631672/80‑1 (MQ=255)                             aaGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGc                                                <  1:631674/80‑1 (MQ=255)                              aGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCa                                               <  1:631675/80‑1 (MQ=255)                               ggTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCaaa                                             >  1:631711/1‑81 (MQ=255)                               ggTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCaaa                                             >  1:631746/1‑81 (MQ=255)                                gTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAAc                                            <  2:631675/81‑1 (MQ=255)                                          tgtCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCta                                   >  1:631729/1‑80 (MQ=255)                                              gtacgtaTTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAgg                              >  1:631730/1‑81 (MQ=255)                                              gtacgtaTTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAgg                              >  1:631741/1‑81 (MQ=255)                                              gtacgtaTTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAgg                              >  2:631725/1‑81 (MQ=255)                                               tacgtaTTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGc                             <  1:631677/81‑1 (MQ=255)                                                    aTTGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCgttgt                        >  2:631764/1‑81 (MQ=255)                                                     ttGGGTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCgttgtt                       <  1:631678/81‑1 (MQ=255)                                                       gggTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCt                     >  2:631743/1‑81 (MQ=255)                                                       gggTTCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCt                     >  2:631723/1‑81 (MQ=255)                                                          ttCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGt                  >  1:631749/1‑81 (MQ=255)                                                          ttCCCTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGt                  <  1:631679/81‑1 (MQ=255)                                                            cccTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGtt                 <  2:631679/80‑1 (MQ=255)                                                            cccTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGtt                 <  1:631681/80‑1 (MQ=255)                                                            cccTCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTa                >  1:631728/1‑81 (MQ=255)                                                            cccTCACCCCAATGGTTAATC‑AAAAGATACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTTTTCTGGTTa                >  1:631723/1‑81 (MQ=255)                                                               tCACCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCg             >  1:631737/1‑81 (MQ=255)                                                                 aCCCCAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCGtt           <  1:631682/81‑1 (MQ=255)                                                                   cccAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCGTTAt         <  1:631683/81‑1 (MQ=255)                                                                    ccAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCGTTAtt        <  2:631683/81‑1 (MQ=255)                                                                     cAATGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCGTTAttt       <  2:631684/81‑1 (MQ=255)                                                                      aaTGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCGTTATTTc      <  2:631686/81‑1 (MQ=255)                                                                      aaTGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCGTTATTTc      <  1:631686/81‑1 (MQ=255)                                                                       aTGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCGTTATTTCa     >  1:631727/1‑81 (MQ=255)                                                                       aTGGTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCGTTATTTCa     >  2:631724/1‑81 (MQ=255)                                                                         ggTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCGTTATTTCATc   >  2:631742/1‑81 (MQ=255)                                                                          gTTAATC‑AAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCGTTATTTCATct  <  2:631687/81‑1 (MQ=255)                                                                                 |                                                                           GAAAAAAGGGTAAGCACGTTATTATGTTAAGGTAACAGACGTGTCGTACGTATTGGGTTCCCTCACCCCAATGGTTAATCAAAAAGGTACAATATGAACGTTTTCTCGCAAACTCAACGCTATAAGGCGTTGTTCTGGTTATCGTTATTTCATCT  >  NC_000913/3609124‑3609278

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