Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	953,650	T→C	9.6%	T273A (ACA→GCA)	focA ←	formate channel FocA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	953,650	0	T	C	9.6%	89.0 / 3.7	42	T273A (ACA→GCA)	focA	formate channel FocA
Reads supporting (aligned to +/- strand):  ref base T (22/16);  new base C (2/2);  total (24/18)
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.10e-01

GTAAAACCTTCCCAGGCTGTGGCTAACTTTTCATTAAGCTCGGACATGTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGAT  >  NC_000913/953508‑953786                                                                                                                                               |                                                                                                                                         gTAAAACCTTCCCAGGCATTGCCTAACTCTTCATTAAGCTCGGACATGTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACGTCGACAACCAGTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGt                                                                                                                                          <  2:254805/143‑1 (MQ=255)        cTTCCCAGGCTGTGGCTAACTTTTCATTAAGCTCGGACATGTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGAACCAATATGTctacgca                                                                                                                                   >  1:97520/1‑143 (MQ=255)              aGGCTGTGGCTAACTTTTCATTAAGCTCGGACATGTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCcaaccaac                                                                                                                             >  1:315688/1‑143 (MQ=255)               ggCTGTGGCTAACTTTTCATTAAGCTCGGACATGTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACcaaca                                                                                                                            >  1:211057/1‑143 (MQ=255)                      gCTAACTTTTCATTAAGCTCGGACATGTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGGTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAaccacc                                                                                                                     >  1:168463/1‑143 (MQ=255)                            ttttCATTAAGCTCGGACATGTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGAt                                                                                                               >  2:498540/1‑143 (MQ=255)                            ttttCATTAAGCTCGGACATGTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAAAAAACCACAACAGac                                                                                                               >  1:563524/1‑142 (MQ=255)                               tCATTAAGCTCGGACATGTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAAt                                                                                                            <  1:335796/143‑1 (MQ=255)                                   tAAGCTCGGACATGTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTg                                                                                                        <  1:27250/143‑1 (MQ=255)                                    aaGCTCGGACATGTCACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTACTGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGc                                                                                                       <  2:377306/143‑1 (MQ=255)                                         cGGACATGTAACACCTACCATATAAAGAGGAATTTCTATTTACTGCTGACTTCGACAACCATCAGTGGTCGTCCTTTTCACGCAGCTAAATGACCCAGTATGTCAACCCAACCAACCACACACCACCACTAAGGGTGCCCATc                                                                                                  >  2:412728/1‑143 (MQ=255)                                            aCATGTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGCTAATGTTGCCGATCGTa                                                                                               >  1:358243/1‑143 (MQ=255)                                                gTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAACTGACCCAGTATGCCGACCCAACCAACGCACCACCACCGATAATGTTGCCGATCGTAACCg                                                                                           >  1:172933/1‑143 (MQ=255)                                                gTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGCCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCg                                                                                           >  1:58748/1‑143 (MQ=255)                                                 tAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCgg                                                                                          >  1:233214/1‑143 (MQ=255)                                                    cacCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAAt                                                                                       <  1:548594/143‑1 (MQ=255)                                                           cttcttAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTCCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATCATGTTGCCTATCGTAGCCGGAATGAGGTTa                                                                                >  1:550370/1‑143 (MQ=255)                                                            ttcttcAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGTTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTAt                                                                               <  2:225193/143‑1 (MQ=255)                                                                         tttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCAAGTATGTCAACCCAACCAACAAACCACCACCGATAATATTGCCGATCGTAACCGGAAGCAGGTTATTAGTGAGGAAAtt                                                                  >  1:413359/1‑143 (MQ=255)                                                                           tttatttaCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATtca                                                                <  1:498540/143‑1 (MQ=255)                                                                            ttatttaCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATtcat                                                               >  1:443863/1‑143 (MQ=255)                                                                                ttaCTGCGTACTTCGACAACCATTAATGGTGGTCGATTTCACGCAGGTAAATGACCCAATATGTCAAAGTATCTGACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACg                                                           >  1:426024/1‑143 (MQ=255)                                                                                  aCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCCAGGCGGTTTATGTACCAGAAAATCTAGTCGGCAGAGAAACCACCACCGATAATGTTGCCGATCGTAACCCGAATCAGGTTATCAGTGATGAAATTCATCACGGt                                                         >  1:497457/1‑143 (MQ=255)                                                                                   cTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACGCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAACCAGGCTACAAGAGACGAAATCCAACACGGTc                                                        >  2:654096/1‑143 (MQ=255)                                                                                     tcGTACTGCGACGACCATTAATGGTGGTCGGTTTCACGCAGATAAATGACCTAGTATGTCAACCCAACCAACAAAGCACCACCGAGAATGTTGCCGAACGTAACCGGAAGCAGGTTATCAGTGATGAAATTCATCACGGTCAg                                                      <  2:566933/142‑1 (MQ=255)                                                                                         agttCGACAACACGTAATGGTGGTCGTCTTCACGAAGGTAAATGACCCAGTATGCCGAACCCCCATAAGAACCACCACCGATAATGTTGCCGATCGTAGACGGTAGCAGGGTAGCAGTGATGAAATTCATCACGGTCAGGTga                                                  <  2:172933/141‑1 (MQ=255)                                                                                         aCTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGCCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTga                                                  <  2:58748/143‑1 (MQ=255)                                                                                                 aaCCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGGAAAAtg                                          >  1:37372/1‑142 (MQ=255)                                                                                                   ccATTAGTGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAAtttt                                        <  2:634637/143‑1 (MQ=255)                                                                                                      ttgaTGGTGGTCCTTTGCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCg                                     <  2:315688/140‑1 (MQ=255)                                                                                                        aaTGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCCCCGCACCCAACAAACCACCACCGATAATGTGGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGt                                   <  1:642101/143‑1 (MQ=255)                                                                                                         attgtggtCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAAGATTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTg                                  <  2:366982/140‑1 (MQ=255)                                                                                                                tCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGGAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCGGGTGAGAAGACTTTTCCGGGGCAGAAcc                           >  1:192118/1‑143 (MQ=255)                                                                                                                   ttttCACGCAGGTAAATGACCCAGTATGTCAACCCAACAAACAAACCACCACCGATAATATTGCCGATCGTAACCGGAATCAGGTTAGCAGTGGTGACATTCAGCACGGTCAGGTGCGAAAAATTTTCCGATGCAGAACCGAc                        >  1:42184/1‑143 (MQ=255)                                                                                                                         cacccgTAAATGACCCAGTATGTCAACCCATCCAACAAACCACCACCGATAATGATGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGt                  <  1:268619/138‑1 (MQ=255)                                                                                                                         cGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGt                  >  1:135213/1‑143 (MQ=255)                                                                                                                             ggTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTGATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCaa              <  1:253202/143‑1 (MQ=255)                                                                                                                               tAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCaaaa            >  1:72045/1‑143 (MQ=255)                                                                                                                                 attGACCCAGTATGTCAACTCAACCAAAAAACCACCACCGATAATGTTGCCGAGCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAAtt          <  1:186173/141‑1 (MQ=255)                                                                                                                                   tgccccAGGATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGCTGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGATCCGACTGCGGTCCAAAATTcc        <  2:630651/140‑1 (MQ=255)                                                                                                                                    gACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCg       <  2:234916/143‑1 (MQ=255)                                                                                                                                       ccAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGCTACCAGAGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTGCGGTGCAGAACCGCCTGCCGTCCAAAAAACCgggg    >  2:484272/1‑143 (MQ=255)                                                                                                                                         aGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCACCATTCCGTGGAt  >  1:45762/1‑143 (MQ=255)                                                                                                                                               |                                                                                                                                         GTAAAACCTTCCCAGGCTGTGGCTAACTTTTCATTAAGCTCGGACATGTAACACCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACCGGAATCAGGTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGAT  >  NC_000913/953508‑953786

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