Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,343,779	A→C	S131A (TCC→GCC)	yfaL ←	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,343,779	0	A	C	81.8%	107.5 / 8.1	44	S131A (TCC→GCC)	yfaL	adhesin
Reads supporting (aligned to +/- strand):  ref base A (0/8);  new base C (27/9);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 1.37e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCT  >  NC_000913/2343634‑2343921                                                                                                                                                  |                                                                                                                                               ttaATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGcaaaa                                                                                                                                            >  2:440840/1‑150 (MQ=255)     tGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATaa                                                                                                                                        >  1:825301/1‑150 (MQ=255)           aaaTTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGt                                                                                                                                 <  2:64319/151‑1 (MQ=255)           aaaTTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGt                                                                                                                                 >  1:693992/1‑151 (MQ=255)            aaTTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGt                                                                                                                                 >  1:592734/1‑150 (MQ=255)              ttGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTa                                                                                                                                >  2:605927/1‑149 (MQ=255)                gcgcCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGa                                                                                                                             <  1:804868/150‑1 (MQ=255)                  gcCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTc                                                                                                                          <  1:397739/151‑1 (MQ=255)                    cACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTc                                                                                                                          >  1:471301/1‑149 (MQ=255)                     aCCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGc                                                                                                                       >  1:415913/1‑151 (MQ=255)                         ttGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGt                                                                                                                     <  1:537400/149‑1 (MQ=255)                          tGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTg                                                                                                                   >  2:133121/1‑150 (MQ=255)                           gCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAg                                                                                                                 <  1:794470/151‑1 (MQ=255)                               tcatTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGtttt                                                                                                             >  2:107254/1‑151 (MQ=255)                                catTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGtttt                                                                                                             <  1:594069/150‑1 (MQ=255)                                  ttAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCtt                                                                                                          >  1:361214/1‑151 (MQ=255)                                       gATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGc                                                                                                      >  1:471189/1‑150 (MQ=255)                                       gATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGc                                                                                                      >  2:341024/1‑150 (MQ=255)                                          gTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCaaa                                                                                                   <  1:697384/150‑1 (MQ=255)                                            tATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAata                                                                                                 <  2:211810/150‑1 (MQ=255)                                             aTTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAatat                                                                                                <  1:70817/150‑1 (MQ=255)                                              ttGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATgg                                                                                              >  1:302585/1‑151 (MQ=255)                                                  gAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGcccc                                                                                          >  1:486483/1‑151 (MQ=255)                                                    aaCATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCccgc                                                                                        >  2:916577/1‑151 (MQ=255)                                                      cATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCccgccg                                                                                      >  1:47536/1‑151 (MQ=255)                                                       aTGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGt                                                                                     >  1:115190/1‑151 (MQ=255)                                                                        tAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGGCGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGt                                                                     >  2:591465/1‑150 (MQ=255)                                                                         aaaCCGCCGCCACCATTATCTTAAGAACAAAAAAAATAGATTCCCCCCCCAACGCCCCCTGCGCCTTTCCCGAAAAAAAAAACATCATAAAGTTCCGGGTTTGTTTTTTCTTTAGAAAATATGCCCCCGCCGTTATTATATTCTCCTGt                                                                     <  2:471963/149‑1 (MQ=255)                                                                          aaTCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGTAAc                                                                  >  2:813109/1‑151 (MQ=255)                                                                              gACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGgtg                                                              >  1:352911/1‑151 (MQ=255)                                                                                   ccccATTATCTTTAAAACAAGAAAAAAAGTTTCCCCCCCCAGACCCGCCTGCGACTTTCCCGAAAAAAAAAACATCAGAAATTTTGGGTTTTTTTTTTTTTTTAAAAAATAGGCCCCCGCCGTTATTATATTCTCCTGTAACg                                                                 <  1:480771/141‑1 (MQ=255)                                                                                    aCCAGTATCATTAGTACAAGAAAAATAGATCCCCCCCCCACACCCGCCTGCGCCGTTACCGAAAAAAAAAACACCAGTAAGTTTCAGTTTTGTTTTTTTTTTAGAAAATAGGGCCCCGCCGTTATTATATTCTCCTGTAACGgtgt                                                             <  1:475306/146‑1 (MQ=255)                                                                                     cgattATCGTAAGAACAAAAAAAATAGTTGGCCCCCCCATCGCCCCTTGCGACTTAACCAAAAAAAAAAACATCAGTAAGATTCTGCGTTTTTTTTTTTTGAGAAAATAGGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGtt                                                            <  1:818976/142‑1 (MQ=255)                                                                                         tATCTTTAATACAAAAAAAATGAATTGCCCCCCAATACCCGCCTGGGACGTAGCCGAAAAAAAAACCATCATTAAGTTTCGATTTTGTTTTTTATAGAGAAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTg                                                        <  1:41458/146‑1 (MQ=255)                                                                                                  tACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCAt                                          >  2:301003/1‑151 (MQ=255)                                                                                                        aaaaaTATATCGCCCCACAATCCCCCCTTGCTACTTTCCCGAAAAAAAAAACATAATTAATTTTCGGTTTTTTTTTTTCTTTAAAAAATGGGCCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTcc                                       <  1:227052/145‑1 (MQ=255)                                                                                                            aTAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTc                                 >  1:700787/1‑150 (MQ=255)                                                                                                                       caccaTAGCCGCCTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTa                      >  2:739607/1‑150 (MQ=255)                                                                                                                         ccaTAGCCGCCGGCGACTTTCCCGAAAAAAAAAACATAAAAAAGTTTCAGTTTTTTTTTTTTTTTAGAAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAAt                    <  2:7814/150‑1 (MQ=255)                                                                                                                                 gccTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAAt           >  2:510402/1‑151 (MQ=255)                                                                                                                                  ccTGCGACGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAAtt          >  1:452153/1‑151 (MQ=255)                                                                                                                                       gACGTTACCGAAAAAAATAACATCAGTAAGATTCAGTGTTTTTTTTTTTTAAGAAAATAGGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTcc        <  2:311641/148‑1 (MQ=255)                                                                                                                                        aCGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTc    >  2:554395/1‑151 (MQ=255)                                                                                                                                         cGTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTcc   >  2:885428/1‑151 (MQ=255)                                                                                                                                          gTTACCGGCAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCt  >  1:55604/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                               TTAATGGTATAAATTGCGCCACCTTTGCCATCATTAGCGATGTTATTGCGAAACATGGCGTTAGTGACACGTAAATCGACGGCACCAGTATCGTTAGTACCAGAAGAATAGATTGCGCCACCATAGCCGCCTGCGACGTTACCGGAAAAAATAACATCAGTAAGATTCAGCGTTGAGTTTTCTTTAGCAAATATGGCCCCGCCGTTATTATATTCTCCTGTAACGGTGTTATTGGCAAACAGAGTCATTCCAGTTTCATTTTCTGGTAATAGTGAAAATTCCGCTCCT  >  NC_000913/2343634‑2343921

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