Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,395,566	C→T	S119S (TCG→TCA)	ynaI ←	small conductance mechanosensitive channel YnaI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,395,566	0	C	T	91.7%	57.8 / ‑1.2	25	S119S (TCG→TCA)	ynaI	small conductance mechanosensitive channel YnaI
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (11/11);  total (13/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 = 4.69e-01

GCCAGCCAATTTCCGCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACATAATTAATAATGCTACCAACAGCACTGATAATTATGGTGTAGATTAAAAATTGCACCGCATGGAATAAAAATCC  >  NC_000913/1395319‑1395808                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   cccagccaATTTCCGCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTAGCCAATACTGAAAGGACGGTCAAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCAGACCGACAGCCATACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGGTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGAc                                                                                                                                                                                                                                                 <  2:94207/250‑1 (MQ=255)                                                cgaaTCCAGTCGCAAAGACTGAAAGGACGTTCGAAATAGAGCATATGCCCGAAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGttt                                                                                                                                                                                                  <  2:11683/248‑1 (MQ=255)                                                   atacAGTCGCTAATACTGAAAGAAAGGACGAAATAGAGCATAATCCCGGAAAAGAAGATACTCAGATTATCTTTACCGACCATACCGACAGCCAGACAACCACAACAACCAACGATCAGCAAGCCAGAAAGGCTAATGCCGACATGATAGCCATAAAGCCGCACTAACACCACAGTAATGGTGATTTAGATGCTACGTGACATACTCAGCGACCTGGTGAAATCGCGACCTTTTTTACACTGCTGTTTTTc                                                                                                                                                                                               >  2:345214/4‑251 (MQ=255)                                                            ccAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGACTATCTTTACCGGCCATACCGACAGCCAGAACACCAAGACAACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGAGACATAATACGCGCACTGGGGATATCGCGACCTTTTTTAATctgc                                                                                                                                                                                                       >  2:533402/1‑234 (MQ=255)                                                            ccAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCCCCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATctgc                                                                                                                                                                                                       <  1:533402/234‑1 (MQ=255)                                                                               cGAAATAGAGCATACTCCCGGAAAAGCAGTTACTCCGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAACGGTCAGCAAGCCAGAAAGGCTAATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTACTCTGCTGGTTTTCAAACTGATTAATCAGCAGACATAGCtt                                                                                                                                                                    >  1:435922/1‑250 (MQ=255)                                                                                gAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCTCTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGAAGAAATAGCTTa                                                                                                                                                                   >  2:429897/1‑250 (MQ=255)                                                                                               cccGGAAAAGAAGTTACTCAGGATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCg                                                                                                                                                   <  1:435331/251‑1 (MQ=255)                                                                                                       agaagTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAAGACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAgg                                                                                                                                           <  2:115466/251‑1 (MQ=255)                                                                                                            ttACTCAGACTATCTTTAACGGCCATACCGACAGCCAGACCACCAATACCACCAAATTTCCGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAACAGCACATGCACCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCAATGGTGCTATCACGACCTTTTTTCATCTGCTGTTTTTCAAACTGGTTACTACGCAGAAATAGCTTAATTGTCAGAATAACAGCAATCAGGGActt                                                                                                                                      >  2:550994/1‑249 (MQ=255)                                                                                                                        cactgaaccgCCATACCGACAGCCAGACCACCAATACCACCAAATGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAAt                                                                                                                           <  2:159306/241‑1 (MQ=255)                                                                                                                         cTTTACCG‑CCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAAt                                                                                                                           >  1:159306/1‑248 (MQ=255)                                                                                                                          tttACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCCGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGACATCATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTa                                                                                                                                                                                     >  2:458112/1‑190 (MQ=255)                                                                                                                          tttACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTa                                                                                                                                                                                     <  1:458112/190‑1 (MQ=255)                                                                                                                                  ccATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCt                                                                                                                >  1:526131/1‑251 (MQ=255)                                                                                                                                         caCAGCCAGACCACCAATACCTCCAAAGGTCAGCAAGCCAGAAAGGCTCATGCTGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGa                                                                                                         <  2:445423/250‑1 (MQ=255)                                                                                                                                              ccAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAAt                                                                                                    >  1:462709/1‑251 (MQ=255)                                                                                                                                                     accaATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATGTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTGTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGtt                                                                                             <  1:294047/251‑1 (MQ=255)                                                                                                                                                                         gCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACAtaa                                                                         <  1:359751/251‑1 (MQ=255)                                                                                                                                                                                            gCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGAGATTTTGATGATACGTGACATAATCCCCGCACTGGTGATATAGCGTCCTTTTTTAATCTGCTGTCTTTCAAACTGATTAATCAGAAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACATAATTAATAATGCTACCAAc                                                        >  1:65096/1‑249 (MQ=255)                                                                                                                                                                                                    gTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACATAATTAATAATGCTACCAACAGCACTGATa                                              <  1:393657/251‑1 (MQ=255)                                                                                                                                                                                                                                     tGGTGATTTTGATGAAAAGTGACATAATCCGCGCACTGGAGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCCGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGGACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACATAATTAATAATGCTACCAACAGCACTGATAATTATGGTGTAGATTAAAAATTGCACCGCATgg             >  2:331887/1‑251 (MQ=255)                                                                                                                                                                                                                                       gTGATTTTGTTGATACGTGACATAATCAGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGACTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACATAATTAATAATGCTACCATCAGCCCTGATAATTATGGTGTAGATTAAAAATTGCACCGCATGGa            >  2:353599/1‑250 (MQ=255)                                                                                                                                                                                                                                         gATTTTGATGATACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACATAATTAATAATGCTACCAACAGCACTGATAATTATGGTGTAGATTAAAAATTGCACCGCATGGAATa         <  1:94832/251‑1 (MQ=255)                                                                                                                                                                                                                                                atgatACGTGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACATAATTAATAATGCTACCAACAGCACTGATAATTATGGTGTAGATTAAAAATTGCACCGCATGGAATAAAAATcc  >  1:60927/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   GCCAGCCAATTTCCGCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGATAACGCCTGGGGTGATAAATTTGAGTTTATAGTTTTCTATAACATAATTAATAATGCTACCAACAGCACTGATAATTATGGTGTAGATTAAAAATTGCACCGCATGGAATAAAAATCC  >  NC_000913/1395319‑1395808

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