Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,392,997	G→A	V357V (GTC→GTT)	queG ←	epoxyqueuosine reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,997	0	G	A	96.3%	80.7 / ‑5.4	28	V357V (GTC→GTT)	queG	epoxyqueuosine reductase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (14/12);  minor base C (1/0);  total (16/12)
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.55e-01

CAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGAGTT  >  NC_000913/4392757‑4393238                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  cAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACg                                                                                                                                                                                                                                         >  1:225548/1‑251 (MQ=255)       ttAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAAGCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCAttt                                                                                                                                                                                                                                    <  2:218328/249‑1 (MQ=255)                  ttCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGATGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTAt                                                                                                                                                                                                                         <  1:61598/250‑1 (MQ=255)                                          tAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGCTCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTAtg                                                                                                                                                                                                >  1:203898/1‑251 (MQ=255)                                                      aaCCTATTGAAAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATGCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCAGTTCGTCTCTCTATATGCTGCGCAATCGCCCACGCTATGTGCACATCgag                                                                                                                                                                                     <  2:203898/250‑1 (MQ=255)                                                                tagataAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGCCGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGAGCTCATCGGGATGTGGGTGCt                                                                                                                                                                          >  2:332409/1‑251 (MQ=255)                                                                                tAAGATTCCGTTTTGTCTAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGCATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGCGAAGTGGGTGCTCACCTTTCCGACTTTc                                                                                                                                                          >  2:264291/1‑251 (MQ=255)                                                                                     ttCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTAAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCCTCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAgcg                                                                                                                                                     >  2:17696/1‑251 (MQ=255)                                                                                              gtcgagtcgaTCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCaaa                                                                                                                                             <  2:289460/250‑1 (MQ=255)                                                                                              gtcgaggcgaTCACTTGACGATTTTATTCGTCTTGAATTGCTATGCTTTTTGATTTTTATTCACAAGCTGGGGATTAAGCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGTCAAGTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGTGCTGTCaaa                                                                                                                                             <  2:257874/250‑1 (MQ=255)                                                                                                            ttGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCa                                                                                                                               <  2:178846/250‑1 (MQ=255)                                                                                                               aCGATTTTATTCGTCTTTAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAggg                                                                                                                           <  1:264291/251‑1 (MQ=255)                                                                                                                                       aaTGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGat                                                                                                    >  1:235841/1‑250 (MQ=255)                                                                                                                                                ttatttttattCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGc                                                                                          >  1:239916/1‑251 (MQ=255)                                                                                                                                                  atttttattCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCa                                                                                        <  2:225548/251‑1 (MQ=255)                                                                                                                                                          tCACATGCTGTGGATGAATCAAGCGTTACGCTGTAACCCTTTTTCAATCACCCTAACCAGACGCTTTTTTTTCGGCAATTGCACTTCAACTATTCACTCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTTCTCACCTTTACTACTTTCCAGCTCTTTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCTATATTACGCAGCCAACGCAgg                                                                                 <  1:333114/250‑1 (MQ=255)                                                                                                                                                           cacaAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCACCCAACGCAGGTg                                                                               >  1:183474/1‑251 (MQ=255)                                                                                                                                                                  tgtgGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGAACTCCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGCCAAACTCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAAt                                                                         >  1:190938/1‑250 (MQ=255)                                                                                                                                                                                gCGTCAGGCGGTACACCTTTTTCACTCACCCGAACCAGACGCTGTTTTTTCGGCACTTGCACTTCCACTATGCACGCCTTTCGTCTCTCTATTTGCTGCGCACTCGCCCAAGCTATGTGGTCATCGAGAAGCGGGTGCTCACCTTTCCGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCCACATTAAGCAACCCAAGCAGGAGACCAATACGACGAATCGCCgt                                                          >  2:196651/1‑250 (MQ=255)                                                                                                                                                                                  gTCACGCGATAACCCGTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTTCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGACGTGGGTGCTCCCCTTTACGACTTTCCCCCGCTGTCACACTCACTTCCTCCCCCGCCGCATTGCCTACTGCTACGGCGCTATTCCCCCGCCAACGCACGTCCCCAATaccacgcctccccgccc                                                        >  2:367295/1‑240 (MQ=255)                                                                                                                                                                                    cACGCGGTAACCCTTTTTCCATCACCCGAGCCAGACGCTGTTTTTTCGGCAATTGCCCTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGACGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTCATACTACGGCGATATTACGCAGCCACCGCCGGTTACCAATCCGACGAATCGCCGATCCt                                                      >  2:173460/1‑251 (MQ=255)                                                                                                                                                                                    cACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCt                                                      >  1:174250/1‑251 (MQ=255)                                                                                                                                                                                        cgGTAACCCTTTTTCAATCACCCGAACCAGACGCGTTTTTTTCGGCAATGGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTGCCg                                                  <  2:259713/251‑1 (MQ=255)                                                                                                                                                                                         ggTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGt                                                 >  1:266646/1‑251 (MQ=255)                                                                                                                                                                                          gTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTg                                                >  1:213269/1‑251 (MQ=255)                                                                                                                                                                                                           aCCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTtctc                               >  2:348470/1‑251 (MQ=255)                                                                                                                                                                                                                                 tCGGCAAGTGCACTTCACCTAGGCACGCATTTCGTCTCTCTAGTTTCTGCGCACTCGCCCCCGCTATGTTCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACTAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCATGCGAATAACTc          <  2:120974/250‑1 (MQ=255)                                                                                                                                                                                                                                         tgctcTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGGGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGAGtt  <  1:165197/246‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  CAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGAGTT  >  NC_000913/4392757‑4393238

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