Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,446,482	C→T	L791L (CTG→TTG)	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,446,482	0	C	T	97.1%	100.2 / ‑4.2	35	L791L (CTG→TTG)	tamB	translocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (21/13);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

GAACTGTGCGTGCCGCAAACTATCGATGCGGGGGCCGAAGG‑GCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGCTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGT  >  NC_000913/4446245‑4446708                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                   gAACTGTGCGTGCCGCAAACTATCGATGCGGGGGCCGAAGG‑GCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCCCCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGtt                                                                                                                                                                                                                       <  1:262151/251‑1 (MQ=255)          gtgcCGCAAACTATCGATGCGGGGGCCGAAGG‑GCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACt                                                                                                                                                                                                              >  1:263975/1‑251 (MQ=255)                            gCGGGGGCCGAAGGAGCGCG‑GGTGGAGAATATCAACCGCTCAGACCTCGCCATGCTTAAACCGTTTATGCCAGAAACCACGCAGGCCCGCGGTATCTACCCGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCATGTAAC‑GCAAACCGTCAACGATGCGGCTTTGCCCGTGGCGGTTCTGACTCTGAATCTGCcggttgaatt                                                                                                                                                                                            <  2:358411/251‑8 (MQ=255)                             cGGGGGCCGAAGG‑GCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGGGAAGGTAAC‑GCAACCCGTCCACGATGCGGCGTTGCCCGAGGCGTTTCAGACACAGAATCTGACggaggactat                                                                                                                                                                                           >  2:269166/1‑247 (MQ=255)                                        gg‑gCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCg                                                                                                                                                                                >  1:359211/1‑251 (MQ=255)                                                       gAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCt                                                                                                                                                                  >  1:107576/1‑251 (MQ=255)                                                                         ccTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGAcc                                                                                                                                                 >  1:129544/1‑250 (MQ=255)                                                                                                    aGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGAccggtgc                                                                                                                     >  1:370566/1‑251 (MQ=255)                                                                                                          cACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGAcaggtgcaggtg                                                                                                                >  1:167802/1‑250 (MQ=255)                                                                                                          cACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAAAGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGCTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTAGTGGAcaagtgcagg                                                                                                                  >  1:185653/1‑248 (MQ=255)                                                                                                                           tCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGgccgc                                                                                              >  1:223863/1‑251 (MQ=255)                                                                                                                            cTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGgccgcc                                                                                             >  1:161309/1‑251 (MQ=255)                                                                                                                                    gTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGGGGCGTTTCAGAAACTGAATCTGACGGCGGCATTGCGTAACAACCGTGCCGCATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGTCAGGTGCAGGTGACCGCTCCGCAAGGCCGCCGTAATCtt                                                                                     >  2:353760/1‑251 (MQ=255)                                                                                                                                        aGCGGACGTTGCCGGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCTCCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGTAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTtggtg                                                                                 <  1:167907/251‑1 (MQ=255)                                                                                                                                              cGTTGGCTGGGACCCCACGAAAGAGGGGCTGCCGCAGGGCATTATCCCCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGAGGCGGCGTTGCCCGGGGCGTTCCAGACACGGAATCTGACGGCGGAATGGCGTAACAACCGGGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACg                                                                           <  2:107576/251‑1 (MQ=255)                                                                                                                                                  gCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGTTGCATGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCaa                                                                       <  1:324878/251‑1 (MQ=255)                                                                                                                                                   ccTGGTACACCAGGAAAGAGGGGCTGCCGCAGGGCAGTATCTCCCTTTCGGGGCGTAACGTGCAGTTAAC‑GCAACCCGTCAACGATGCGGCTTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGGAACAACCGTGCCTAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAATGCCGCCGTAATCTTGGTGGCAACGTCAat                                                                      <  1:228955/251‑1 (MQ=255)                                                                                                                                                     ggggACCCCACGAAAGAGGGGCTGCCGCAGGGCAGTATCCCCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAatat                                                                    <  1:92492/250‑1 (MQ=255)                                                                                                                                                       ggACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCTCCCTTGCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATcc                                                                  <  1:223543/251‑1 (MQ=255)                                                                                                                                                         acacCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGt                                                                >  1:41955/1‑251 (MQ=255)                                                                                                                                                                gAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAgg                                                                                                  <  2:83525/210‑1 (MQ=255)                                                                                                                                                                gAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAgg                                                                                                  >  1:83525/1‑210 (MQ=255)                                                                                                                                                                gAAAGAGGGGCGGCCGCATGGCAGTACCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑CCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTCCAGCCACTGAATCTGACGGCGGATTTGCTTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGTCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCa                                                         <  1:372590/251‑1 (MQ=255)                                                                                                                                                                     gggggttggcGCAGGGCCGTTTCACCCTTGCGGGGCGTAACTTGAAGG‑AACAGCAAACCGTCAATGAGGCGGCGTTGCCCGTGGTGTTTCAGACTCGGAATCGGCCGGCGGAATGGCTTAACAACCGTGCCGAAGTGGGCTGGACCTTCCGCTGGACAAATAACGGCCAGTTTGATGGACCGGTGGTGGTGACCGTTCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCtg                                                    <  2:359211/242‑2 (MQ=255)                                                                                                                                                                       gggCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGt                                                                                                                                                                               <  1:115385/126‑1 (MQ=255)                                                                                                                                                                       gggCTGCCGCAGCGCAGTATCACCCTTTAGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTAAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGAAACAACCGt                                                                                                                                                                               >  2:115385/1‑126 (MQ=255)                                                                                                                                                                          cTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGAGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCgat                                               >  1:166497/1‑251 (MQ=255)                                                                                                                                                                                cAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGCTGATaa                                          >  1:69230/1‑250 (MQ=255)                                                                                                                                                                                cAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGAGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATaa                                          >  1:367184/1‑250 (MQ=255)                                                                                                                                                                                 aGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCGGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGGGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGCCAATATACGTAACTTCAACCTTGCGATGAAAAAc                                        >  1:51159/1‑251 (MQ=255)                                                                                                                                                                                         aTCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCCATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCtt                                >  1:95436/1‑251 (MQ=255)                                                                                                                                                                                                                gCAGGTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCACCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCgg         <  1:269166/251‑1 (MQ=255)                                                                                                                                                                                                                    gTAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGAt     >  1:381185/1‑251 (MQ=255)                                                                                                                                                                                                                     tAAC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGTATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGACGGACAGGTGCAGGTGACCGATCCGCAAGGCACCCGAAATCTTGGTGGCCACGTCGATATCCGTAACTTCAACCTTGCGATGTTAAACCCCATCTTTACCCGCGGGGAAAAAGCAGAGGGGacg    >  2:297058/1‑249 (MQ=255)                                                                                                                                                                                                                       aC‑GCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGCCACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATCACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGCTGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGt  >  1:166559/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                   GAACTGTGCGTGCCGCAAACTATCGATGCGGGGGCCGAAGG‑GCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGTAACGTGCAGGTAAC‑GCAAACCGTCAACGATGCGGCGCTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGT  >  NC_000913/4446245‑4446708

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