Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,254,980	A→G	F325L (TTT→CTT)	ycgV ←	putative adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,254,980	0	A	G	96.3%	178.1 / ‑1.8	54	F325L (TTT→CTT)	ycgV	putative adhesin
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (25/27);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

GTATTTTGCATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAAATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGTTGTACCACCACGAACT  >  NC_000913/1254883‑1255077                                                                                                  |                                                                                                  gTATTTTGCATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGAt                                                                                                 <  2:14091/100‑1 (MQ=255)  tATTTTGCATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGAtt                                                                                                <  2:776287/100‑1 (MQ=255)   aTTTTGCATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTg                                                                                               <  2:1300128/100‑1 (MQ=255)    ttttGAATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTg                                                                                               <  2:785148/99‑1 (MQ=255)     tttGCATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGAt                                                                                             <  1:241367/100‑1 (MQ=255)     tttGCATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGAt                                                                                             <  1:1535624/100‑1 (MQ=255)      ttGCATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAgattgatt                                                                                            <  1:1223069/100‑1 (MQ=255)      ttGCATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAgattgatt                                                                                            <  2:718914/100‑1 (MQ=255)       tGCATGTTGACAACAGCCGTTGCTGGCTGGGCCGAGTGACCATAAGAGCCGCCGGAAAAGAGGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGAGTGAt                                                                                             <  2:1250118/98‑1 (MQ=255)        gCATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTAt                                                                                          >  2:516195/1‑100 (MQ=255)         cATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGACAAGATGCTGTTTCGGTGCGCTGCCGTGCCCGTAAGATTGATTATc                                                                                         >  2:165361/1‑100 (MQ=255)         cATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATc                                                                                         <  1:899328/100‑1 (MQ=255)              tGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTc                                                                                    <  1:1514941/100‑1 (MQ=255)              tGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCAGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATccagta                                                                                    >  1:1481265/1‑95 (MQ=255)                   acaGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAAc                                                                               >  1:1114044/1‑100 (MQ=255)                    caGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAAcc                                                                              <  2:1114037/100‑1 (MQ=255)                    caGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAAcc                                                                              >  2:649780/1‑100 (MQ=255)                    caGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAAcc                                                                              >  2:649830/1‑100 (MQ=255)                      gCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCCTCTGAACCAc                                                                            >  2:643710/1‑100 (MQ=255)                       ccGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACt                                                                           <  1:411208/100‑1 (MQ=255)                        cGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTg                                                                          >  2:1139121/1‑100 (MQ=255)                                  tGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAAATTGATTATCGCGTCTGAACCACTGGTGACGAGGc                                                                >  2:61625/1‑100 (MQ=255)                                   gggCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGcc                                                               >  2:743148/1‑100 (MQ=255)                                     gCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCAc                                                             <  1:1286418/100‑1 (MQ=255)                                        gAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGc                                                          <  1:715933/100‑1 (MQ=255)                                            cACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTg                                                      >  2:700138/1‑100 (MQ=255)                                             aCCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTgc                                                     >  1:543113/1‑100 (MQ=255)                                               cATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTgc                                                     >  2:383088/1‑98 (MQ=255)                                                 tAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCgg                                                 <  2:563292/100‑1 (MQ=255)                                                  aaGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGga                                                <  2:97544/100‑1 (MQ=255)                                                   aGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGgaa                                               >  1:367909/1‑100 (MQ=255)                                                      aCCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGgaagaa                                            <  1:886865/100‑1 (MQ=255)                                                      aCCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGgaagaa                                            >  1:707528/1‑100 (MQ=255)                                                          ccgGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAatat                                        >  1:191488/1‑100 (MQ=255)                                                                 aGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATc                                 >  1:1087962/1‑100 (MQ=255)                                                                 aGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATc                                 <  2:707525/100‑1 (MQ=255)                                                                  gATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCt                                <  1:593910/100‑1 (MQ=255)                                                                     gCTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCa                             >  2:429089/1‑100 (MQ=255)                                                                      cTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCAc                            <  1:601478/100‑1 (MQ=255)                                                                      cTGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCAc                            <  1:602303/100‑1 (MQ=255)                                                                       tGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCAcc                           <  1:1400642/100‑1 (MQ=255)                                                                       tGTTTCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCAcc                           <  1:1313777/100‑1 (MQ=255)                                                                         tttCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGa                         >  1:281966/1‑100 (MQ=255)                                                                           tCGTTGCGCTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATa                       >  2:1567532/1‑100 (MQ=255)                                                                           tCGTTGCGCTGCCGGGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGaga                       <  2:184799/100‑3 (MQ=255)                                                                                gcgcTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGGGGAAGAAATATGGCTATCTGCACCGATAATGGt                  >  2:951181/1‑100 (MQ=255)                                                                                gcgcTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGt                  <  2:1425130/100‑1 (MQ=255)                                                                                gcgcTGCCGTGCCAGTAAAATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGt                  >  2:127163/1‑100 (MQ=255)                                                                                 cgcTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGtt                 >  2:1121958/1‑100 (MQ=255)                                                                                 cgcTGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGtt                 >  1:781739/1‑100 (MQ=255)                                                                                    tGCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGTTGTa              <  2:1310389/100‑1 (MQ=255)                                                                                     gCCGTGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGTTGTac             >  1:1466259/1‑100 (MQ=255)                                                                                         tGCCAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGTTGTaccacc         >  2:218130/1‑100 (MQ=255)                                                                                           ccAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGTTGTACcaccac       <  1:1531544/100‑1 (MQ=255)                                                                                           ccAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGTTGTACcaccac       >  1:737453/1‑100 (MQ=255)                                                                                            cAGTAAGATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGTTGTACCACCACg      <  1:1177916/100‑1 (MQ=255)                                                                                                aaaaTTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGTTGTACCACCACGAACt  <  2:1068848/100‑1 (MQ=255)                                                                                                  |                                                                                                  GTATTTTGCATGTTGACAACAGCCGTTGCCGTCTGGGCCGAGGCACCATAAGAACCGCCGGAAAAGATGCTGTTTCGTTGCGCTGCCGTGCCAGTAAAATTGATTATCGCGTCTGAACCACTGGTGACGAGGCCACCGCCCTGCGCGGAAGAAATATGGCTATCTGCACCGATAGTGGTTGTACCACCACGAACT  >  NC_000913/1254883‑1255077

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