Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	216,568	C→T	A146T (GCC→ACC)	yaeF ←	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	216,568	0	C	T	97.6%	289.8 / ‑2.1	84	A146T (GCC→ACC)	yaeF	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (41/41);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

GAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAG  >  NC_000913/216473‑216658                                                                                                |                                                                                           gaACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAAc                                                                                        >  2:449716/1‑100 (MQ=255)  aaCGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAAc                                                                                        <  1:1078348/99‑1 (MQ=255)   aCGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAAcgc                                                                                      >  2:197841/1‑100 (MQ=255)   aCGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAAcgc                                                                                      >  2:739482/1‑100 (MQ=255)    cGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAAcgcg                                                                                     >  1:444628/1‑100 (MQ=255)      gatctaGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGt                                                                                   <  2:1726681/95‑1 (MQ=255)      gATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATCAATTCGCCAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGCTAAACGCgg                                                                                    >  2:533698/1‑99 (MQ=255)      gATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGt                                                                                   <  2:1842316/100‑1 (MQ=255)      gATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGt                                                                                   <  1:1781079/100‑1 (MQ=255)      gATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGt                                                                                   >  2:1958018/1‑100 (MQ=255)       cttCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTg                                                                                  <  2:425286/99‑1 (MQ=255)       aTTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTg                                                                                  <  2:1772282/100‑1 (MQ=255)       aTTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTg                                                                                  <  2:2246898/100‑1 (MQ=255)         tCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGat                                                                                <  1:1349274/100‑1 (MQ=255)         tCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGag                                                                                >  1:425288/1‑99 (MQ=255)          cAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGTCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGata                                                                               <  1:691262/100‑1 (MQ=255)          cAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGata                                                                               >  1:1545631/1‑100 (MQ=255)           aGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGtaga                                                                                >  1:1726686/1‑95 (MQ=255)               aGCACATCTGGCGAGTTACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgt                                                                          >  1:108758/1‑100 (MQ=255)                gCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgtg                                                                         >  2:1359219/1‑100 (MQ=255)                gCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgtg                                                                         <  1:900389/100‑1 (MQ=255)                gCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgtg                                                                         <  1:1053321/100‑1 (MQ=255)                gCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgtg                                                                         <  1:1359222/100‑1 (MQ=255)                gCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCtgtg                                                                         <  1:1773870/100‑1 (MQ=255)                   caTCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCt                                                                      >  1:777351/1‑100 (MQ=255)                   caTCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCt                                                                      >  2:1506548/1‑100 (MQ=255)                     tCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttg                                                                    <  2:1545627/100‑1 (MQ=255)                      cTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgt                                                                   >  2:801302/1‑100 (MQ=255)                      cTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgt                                                                   >  2:1335999/1‑100 (MQ=255)                        ggCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCttgttg                                                                 >  1:591527/1‑100 (MQ=255)                            aGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCggg                                                             <  1:2268421/100‑1 (MQ=255)                            aGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCgg                                                              >  2:68023/1‑99 (MQ=255)                             gTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCgggg                                                            <  1:739483/100‑1 (MQ=255)                             gTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCgggg                                                            <  1:868170/100‑1 (MQ=255)                             gTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCgggg                                                            <  1:68025/100‑1 (MQ=255)                               cACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTa                                                          >  1:2223637/1‑100 (MQ=255)                                aCCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTaa                                                         >  2:629873/1‑100 (MQ=255)                                 ccATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAg                                                        >  2:1573326/1‑100 (MQ=255)                                 ccATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAg                                                        >  1:2230934/1‑100 (MQ=255)                                  cATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGa                                                       >  1:895575/1‑100 (MQ=255)                                    tAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATc                                                     <  2:1869140/100‑1 (MQ=255)                                       aGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCgg                                                  <  1:1958025/100‑1 (MQ=255)                                        gggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCggg                                                 <  1:584629/100‑1 (MQ=255)                                         ggAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGa                                                >  1:2270284/1‑100 (MQ=255)                                          gAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGAc                                               >  1:723026/1‑100 (MQ=255)                                          gAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGAc                                               <  1:197843/100‑1 (MQ=255)                                            aTAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTc                                             >  1:1772919/1‑100 (MQ=255)                                             tAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCg                                            <  1:235711/100‑1 (MQ=255)                                              aaaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGt                                           >  1:541193/1‑100 (MQ=255)                                              aaaTTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGt                                           <  2:438325/100‑1 (MQ=255)                                                 ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAg                                        <  2:2270354/100‑1 (MQ=255)                                                 ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAg                                        <  1:220973/100‑1 (MQ=255)                                                 ttCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAg                                        <  1:89952/100‑1 (MQ=255)                                                    gACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACg                                     >  1:835332/1‑100 (MQ=255)                                                      cAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTGAGATCCGGGACTCGTAAGACGaa                                   >  1:578994/1‑100 (MQ=255)                                                        aTGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGaaaa                                 >  1:100141/1‑100 (MQ=255)                                                          gCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGc                               >  1:19068/1‑100 (MQ=255)                                                           ccGCGATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCt                              <  2:74093/100‑1 (MQ=255)                                                            cgcgATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCtt                             >  2:385403/1‑100 (MQ=255)                                                            cgcgATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCtt                             <  1:1792775/100‑1 (MQ=255)                                                             gcgATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTa                            <  2:444626/100‑1 (MQ=255)                                                              cgATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTACGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTa                            <  2:1246833/99‑1 (MQ=255)                                                               gATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATc                          <  2:1772914/100‑1 (MQ=255)                                                               gATAGTTATAACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATc                          <  2:777350/100‑1 (MQ=255)                                                                   ttccgatctgCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTa                      <  2:754534/91‑1 (MQ=255)                                                                       taACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCt                  <  1:813141/100‑1 (MQ=255)                                                                       taACCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCt                  <  1:736025/100‑1 (MQ=255)                                                                        aaCCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCtt                 >  2:1494411/1‑100 (MQ=255)                                                                         aCCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTc                <  1:801303/100‑1 (MQ=255)                                                                         aCCGCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTc                <  1:1336002/100‑1 (MQ=255)                                                                            gCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTaagatcggaa             >  1:754535/1‑91 (MQ=255)                                                                            gCTATCTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATc             <  2:2014536/100‑1 (MQ=255)                                                                                 cTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttt        >  1:877272/1‑100 (MQ=255)                                                                                 cTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttt        <  2:412242/100‑1 (MQ=255)                                                                                 cTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttt        >  2:1752602/1‑100 (MQ=255)                                                                                 cTTTGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttt        >  2:1261988/1‑100 (MQ=255)                                                                                  tttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCtttt       >  1:1778062/1‑100 (MQ=255)                                                                                   ttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttttt      >  2:1180032/1‑100 (MQ=255)                                                                                   ttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttttt      <  2:723025/100‑1 (MQ=255)                                                                                   ttGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCttttt      >  1:990594/1‑100 (MQ=255)                                                                                    tGATTTTATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCtttttt     <  2:108756/100‑1 (MQ=255)                                                                                       ttttATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAg  <  1:449718/100‑1 (MQ=255)                                                                                       ttttATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAg  >  1:1921007/1‑100 (MQ=255)                                                                                       ttttATTGGTAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAg  >  2:946386/1‑100 (MQ=255)                                                                                                |                                                                                           GAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCGCGATAGTTATAACCGCTATCTTTGATTTTATTGGCAAACGCGGTGATATCTGTGGCTTGTTGCGGGGTAAGATCCGGGACTCGTAAGACGAAAAGCTTATCACTATGCTTCATCGCTTTTTTAAG  >  NC_000913/216473‑216658

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