Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	715,849	A→G	L146L (TTA→TTG)	ybfP →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	715,849	0	A	G	97.1%	224.5 / ‑3.4	68	L146L (TTA→TTG)	ybfP	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (29/37);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACAT  >  NC_000913/715752‑715943                                                                                                  |                                                                                               gcgaGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGTATTATTGAGTTgcg                                                                                              <  1:792392/97‑1 (MQ=255) gCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcg                                                                                              >  2:1743179/1‑100 (MQ=255)   cAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcgac                                                                                            >  1:216008/1‑100 (MQ=255)   cAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTgcgac                                                                                            >  2:45057/1‑100 (MQ=255)      ggCAAGCCCGGCACTCGACCGATAGTCAATGAGGCTTATGTCAGCGAAAGCGGATTTTTTGTCTCGATTATGCATGTCGGGATTATTGAGTTGCGACGc                                                                                          >  2:1414457/1‑99 (MQ=255)       gCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTg                                                                                        >  2:1521395/1‑100 (MQ=255)       gCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTg                                                                                        <  2:1403310/100‑1 (MQ=255)       gCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTg                                                                                        >  1:1527850/1‑100 (MQ=255)          aGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCa                                                                                     >  1:869303/1‑100 (MQ=255)          aGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCa                                                                                     >  1:1719833/1‑100 (MQ=255)           gCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAg                                                                                    <  2:1342619/100‑1 (MQ=255)             tcGGCACTCGACCGATAGTTAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTc                                                                                  <  2:180376/99‑1 (MQ=255)              cGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCa                                                                                 >  2:793933/1‑100 (MQ=255)               ggCACTCGACTGATAGTCAATGAGGATTATTTCAGCGAAAGGGGATTTTTTGGCACGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTc                                                                                  <  2:680961/98‑1 (MQ=255)                 cACTCGACCGATAGTCAATGAGGATTATGTCAGAGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAg                                                                                    <  2:429822/94‑1 (MQ=255)                      gACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTg                                                                         <  1:492879/100‑1 (MQ=255)                       aCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGc                                                                        >  2:133575/1‑100 (MQ=255)                        ccGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGcc                                                                       <  1:1279900/100‑1 (MQ=255)                        ccGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGcc                                                                       <  1:220936/100‑1 (MQ=255)                         cGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTTTTTTTTTTTTTCATGTCGGGATTATTGAGTTGTGACGCTGTGAGTTAGAAAATTGcga                                                                      >  1:1498442/1‑98 (MQ=255)                         cGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCa                                                                      >  1:1488622/1‑100 (MQ=255)                                   tGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTg                                                            >  2:8495/1‑100 (MQ=255)                                        aTTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAg                                                        <  1:1508844/99‑1 (MQ=255)                                         ttATGTCATCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCg                                                      <  2:288882/100‑1 (MQ=255)                                         ttATGTCAGCGAAAGCGGATTTTTTGTCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGATAATTGCCAGCAAGAATTGAAAGCg                                                      <  1:690314/100‑1 (MQ=255)                                         ttATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCg                                                      <  1:45057/100‑1 (MQ=255)                                           aTGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAt                                                    >  2:87185/1‑100 (MQ=255)                                           aTGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAt                                                    >  2:570568/1‑100 (MQ=255)                                           aTGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAt                                                    <  2:209125/100‑1 (MQ=255)                                           aTGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAt                                                    <  1:1537714/100‑1 (MQ=255)                                             gTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATaa                                                  >  1:936644/1‑100 (MQ=255)                                             gTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATaa                                                  >  1:1137922/1‑100 (MQ=255)                                             gTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGACAGCAAGAATTGAAAGCGATaa                                                  >  2:602514/1‑100 (MQ=255)                                                 gCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGTATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATac                                              <  2:835251/100‑1 (MQ=255)                                                 gCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATac                                              <  2:1805713/100‑1 (MQ=255)                                                 gCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGATGCCAGTCAGAAAATTGACAGCAAGAATTGAAAGCGATAAATac                                              <  1:898394/100‑1 (MQ=255)                                                  cGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATaca                                             <  1:467400/100‑1 (MQ=255)                                                  cGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATaca                                             <  2:630167/100‑1 (MQ=255)                                                   gAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATacac                                            <  1:1422023/100‑1 (MQ=255)                                                      aGCGGATTTTTTTTCTCTATGATGCATGTCAGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACAtt                                         >  1:1010454/1‑100 (MQ=255)                                                           aTTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAAt                                    >  2:1052600/1‑100 (MQ=255)                                                           aTTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAAt                                    <  2:563093/100‑1 (MQ=255)                                                           aTTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAAt                                    >  2:1426179/1‑100 (MQ=255)                                                           aTTTTTTGGCTCGATGATGAATGTCGGGATTATTGAGTTGCGACGCTGCCAGTAAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACAAATTAAAAt                                    >  1:675690/1‑100 (MQ=255)                                                             tttttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAg                                  <  1:734712/100‑1 (MQ=255)                                                               tttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTg                                <  2:1212952/100‑1 (MQ=255)                                                               tttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTg                                <  1:1740200/100‑1 (MQ=255)                                                               tttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAAGTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTg                                <  2:436891/100‑1 (MQ=255)                                                               tttGGCTCGATGATGCATATCGGGATTATTTAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGt                                 <  1:1798348/99‑1 (MQ=255)                                                                ttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGa                               <  2:905582/100‑1 (MQ=255)                                                                ttGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGa                               >  1:702073/1‑100 (MQ=255)                                                                 tGGCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAg                              <  2:1291099/100‑1 (MQ=255)                                                                 atgccggatgatGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAg                              <  1:1803923/94‑1 (MQ=255)                                                                   gCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGt                             >  2:1760890/1‑99 (MQ=255)                                                                   gCTCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTc                            >  2:1585275/1‑100 (MQ=255)                                                                     tCGATGATGCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCg                          >  1:425320/1‑100 (MQ=255)                                                                        atgatgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGa                       <  2:156444/100‑1 (MQ=255)                                                                         tgatgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGaa                      >  2:1310892/1‑100 (MQ=255)                                                                          gatgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAAt                     <  2:1498453/100‑1 (MQ=255)                                                                          gatgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAAt                     <  2:1527861/100‑1 (MQ=255)                                                                           atgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATg                    <  1:558779/100‑1 (MQ=255)                                                                           atgCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGACAGCAAGAATTGAAAGCAATAAATATACATTAAAATAGTGAGTCCGGGAATg                    <  1:441046/100‑1 (MQ=255)                                                                            tgtaTGTCGGGATTATTGAGTTGCGCTGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATTCACATTAAAAGAGTGAGTCCGGGAATg                    <  1:35938/96‑1 (MQ=255)                                                                             gCATGTCGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGaa                  >  1:460749/1‑100 (MQ=255)                                                                                   cGGGATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCg             >  1:223374/1‑99 (MQ=255)                                                                                      gATTATTGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGAcc         <  2:1488633/100‑1 (MQ=255)                                                                                         tattGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGa      <  1:1521384/100‑1 (MQ=255)                                                                                          attGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGaa     >  2:1333661/1‑100 (MQ=255)                                                                                           ttGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAAc    >  2:1353269/1‑100 (MQ=255)                                                                                            tGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACa   >  2:1098626/1‑100 (MQ=255)                                                                                            tGAGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACa   >  1:129850/1‑100 (MQ=255)                                                                                              aGTTGCGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACAt  <  2:629894/99‑1 (MQ=255)                                                                                                  |                                                                                               GCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAGTGAGTCCGGGAATGAAACCCGGACCTGAACAT  >  NC_000913/715752‑715943

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