Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,445,315	C→T	P402S (CCC→TCC)	tamB →	translocation and assembly module subunit TamB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,445,315	0	C	T	92.6%	66.7 / ‑6.5	28	P402S (CCC→TCC)	tamB	translocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/13);  minor base G (1/1);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATG  >  NC_000913/4445068‑4445560                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      cgcgcTGCGCGAACAGCTGGAGATTGGCGGTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATcccgtccg                                                                                                                                                                                                                                                    <  1:140044/251‑1 (MQ=255)                     agcttGGCGTTACTCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTTGCGGAAGACGGATTTCCGCTCAACGTGGAAGTGCACAGAAAACCGCTTTACTGGCCGTTCACTGGTGCGACGCAGTATCAGGCGGCTGATCTGTAACTGAAACTTACCTGCAAAATGACAGATTACACGCTCTCACTGCGTACGGCAGTGAAAGGACAGGAGATCCCCTCAGAGACTATTGCCCTCGAAGcc                                                                                                                                                                                                                                >  2:129657/4‑251 (MQ=255)                         tGGCGGTAATCTGTCCGGTCCGGTGGAGATGGATGTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTGTACTGGCCGTTCACTGGTGAGAAGCAGTAGCAGGCGGATGATCTGAAACGGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGGGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCaaa                                                                                                                                                                                                                             <  1:162508/250‑1 (MQ=255)                            cGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTa                                                                                                                                                                                                                         <  1:31961/251‑1 (MQ=255)                               tAATCTTTGCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTAGTACCCTCGACGCAAAAGGTAAt                                                                                                                                                                                                                       <  2:47451/250‑1 (MQ=255)                                             cggtGGATATGGATTTACGCGCCCAGACGCGACGGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAAt                                                                                                                                                                                                        <  2:159649/251‑1 (MQ=255)                                                     agggagttaagcgaccataCGCGACTGGCGGAAGCCGGAGTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGTTGAGATGCAGTATCAGGCGGATGATCTGAAACGGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGTTCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACATCAGATCATTCTCGagac                                                                                                                                                                                                <  2:224479/233‑4 (MQ=255)                                                                  cccAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGCCAATCTCGACAAACTCACCGTcgcg                                                                                                                                                                                   <  1:43753/251‑1 (MQ=255)                                                                        cgcgACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTg                                                                                                                                                                             >  1:8680/1‑251 (MQ=255)                                                                                  ggAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCCCTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGGAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGaaa                                                                                                                                                                    >  2:170473/1‑250 (MQ=255)                                                                                                 gCTCAACGTGGCAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGAGAGAAGCAGTATCAGGCGGGAGCTATGAGACTGACACTTACCAGCAAAATGACCGATTACACGATCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCAAGACCGCGACAATTACCCTCGACGCCAACGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTAACGGCGt                                                                                                                                                    >  1:266038/1‑251 (MQ=255)                                                                                                        gTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCg                                                                                                                                              >  2:96796/1‑250 (MQ=255)                                                                                                          ggAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGCTTACACGCTCTCTATGCGTACGGCAGTGCAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAACGGTAATGAACAGCAGGTCAATCTCGCCAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACt                                                                                                                                           >  1:135369/1‑251 (MQ=255)                                                                                                                    aGCAAACAGCTTTACTGGCCGGTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCAATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGGCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGAAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGATGCACGACTGGCAACAGGc                                                                                                                                 >  2:229852/1‑251 (MQ=255)                                                                                                                       aaaCAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACACACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCAt                                                                                                                              >  2:77918/1‑251 (MQ=255)                                                                                                                                  atggacgtcaaCTGGTGAGAAGAAGTATCAGGCGGATGATCTGGCCATGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGAAATGAACAGCAGGGCAATCTCGACAAACTCACCGTAGCGGCGCTGGAAGGGAAAACTGAACTCACGGCGATGCTGGGCTGGCGACAGGACATTAGTTGGcgcg                                                                                                                   >  2:167841/11‑251 (MQ=255)                                                                                                                                        gTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCCCCGTCGCGGCGCTGGAAGGGAAAAATGAAATCAAGGCGTTGCTCGAATGGAAACAGGCCATTAGTTGGCGCGGTGAGc                                                                                                             >  2:189650/1‑251 (MQ=255)                                                                                                                                            aCTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAAc                                                                                                         >  1:101402/1‑251 (MQ=255)                                                                                                                                                            aTCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACGCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAac                                                                                         <  2:161688/251‑1 (MQ=255)                                                                                                                                                                  cGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGcc                                                                                   <  2:66202/251‑1 (MQ=255)                                                                                                                                                                             aaactgaaactTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGACGCGGCGAGCTAACGCTTAACGGCATTAACACCGCCACAGAGTTccc                                                                        >  2:182920/1‑251 (MQ=255)                                                                                                                                                                                  gaaactTTCCGGCAAAATGACCGATTACACGCGCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAAGCTCGACAAACTCACCGTCGCGTCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGGGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGAc                                                                    <  2:143224/250‑1 (MQ=255)                                                                                                                                                                                         aCCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGgtcgtc                                                            <  2:48718/251‑1 (MQ=255)                                                                                                                                                                                          ccGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGCCGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAACGAGTTCCCGGACTGGCCGTCg                                                           >  1:248111/1‑251 (MQ=255)                                                                                                                                                                                               aaaaTGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATAAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTACCGCTTAACGGAATTAACAACGCCAAAGAGTTCCCGGACTGGCCGTCGaacct                                                      >  1:239217/1‑248 (MQ=255)                                                                                                                                                                                                                       aTGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTCCCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCcg                                 <  1:47717/248‑1 (MQ=255)                                                                                                                                                                                                                       aTGCGTACGGCAGTGAAGGGACAGGAGATCCCGTCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCgg                              <  1:54091/251‑1 (MQ=255)                                                                                                                                                                                                                         gCGTACGGCGGGGAAGGGACAGAAGATCCCGGACGCGACTATTACCCGCGACGCCAAAGGTAATGAACAGCAGATAAATCTCGACAAGCTCACCGGCGCGGGGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCGCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGTCATTAGCACCGCCAAAGAGTTCCCGGACTGGCCGTGGAAGCTCAATGGCTTGATTAAAACCCGCggga                            <  2:137597/251‑3 (MQ=255)                                                                                                                                                                                                                                           aCAGGAGATCCCGTCCGCGACTATTGCCCGCGACGCAAAAGGTAATGAACAGCAGGTCAATCCCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCCATGGCTTGATTAAACCCCGCGGTAGCCTGTACGGCGGCACCt          <  2:5291/251‑1 (MQ=255)                                                                                                                                                                                                                                                    cccgtccgCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACGCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATg  <  2:8680/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      CGCGCTGCGCGAACAGCTGGAGATTGGCGTTAATCTTTCCGGTCCGGTGGATATGGATTTACGCGCCCAGACGCGACTGGCGGAAGCCGGATTGCCGCTCAACGTGGAAGTGAACAGCAAACAGCTTTACTGGCCGTTCACTGGTGAGAAGCAGTATCAGGCGGATGATCTGAAACTGAAACTTACCGGCAAAATGACCGATTACACGCTCTCTATGCGTACGGCAGTGAAGGGACAGGAGATCCCGCCCGCGACTATTACCCTCGACGCCAAAGGTAATGAACAGCAGGTCAATCTCGACAAACTCACCGTCGCGGCGCTGGAAGGGAAAACTGAACTCAAGGCGTTGCTCGACTGGCAACAGGCCATTAGTTGGCGCGGTGAGCTAACGCTTAACGGCATTAACACCGCCAAAGAGTTCCCGGACTGGCCGTCGAAACTCAATGGCTTGATTAAAACCCGCGGTAGCCTGTACGGCGGCACCTGGCAGATG  >  NC_000913/4445068‑4445560

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