Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 217,037 A→T 6.4% intergenic (‑34/+20) yaeF ← / ← proS putative lipoprotein/prolyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913217,0370AT 6.4% 216.2 / 3.0 61intergenic (‑34/+20)yaeF/proSputative lipoprotein/prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base T (2/2);  total (31/30)
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.00e+00

TGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCAC  >  NC_000913/216911‑217179
                                                                                                                              |                                                                                                                                              
tGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCa                                                                                                                          <  4:1436547/149‑1 (MQ=255)
    gggTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCt                                                                                                                      >  1:617281/1‑149 (MQ=255)
          ggctggctGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATc                                                                                                                <  2:617281/149‑1 (MQ=255)
                  gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGttt                                                                                                           <  1:863428/146‑1 (MQ=255)
                  gATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGttt                                                                                                           >  2:863428/1‑146 (MQ=255)
                      tCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGa                                                                                                    >  2:9988/1‑149 (MQ=255)
                        aaCGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGata                                                                                                  >  2:588211/1‑149 (MQ=255)
                         aCGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGatat                                                                                                 >  3:1350189/1‑149 (MQ=255)
                           ggTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTc                                                                                               <  4:1350189/149‑1 (MQ=255)
                            gTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTgg                                                                                                                               >  1:1351687/1‑116 (MQ=255)
                            gTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTgg                                                                                                                               <  2:1351687/116‑1 (MQ=255)
                                   ccGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgac                                                                                            <  2:93810/144‑1 (MQ=255)
                                   ccGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgac                                                                                            >  1:93810/1‑144 (MQ=255)
                                   ccGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGt                                                                                       >  2:150613/1‑149 (MQ=255)
                                        aGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCa                                                                                  <  1:9988/149‑1 (MQ=255)
                                         gTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGaaa                                                                                                                                                <  1:246800/86‑1 (MQ=255)
                                         gTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGaaa                                                                                                                                                >  2:246800/1‑86 (MQ=255)
                                           aaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacga                                                                                          <  4:577864/138‑1 (MQ=255)
                                           aaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacga                                                                                          >  3:577864/1‑138 (MQ=255)
                                               gaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaa                                                                           <  1:150613/149‑1 (MQ=255)
                                                     tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAAc                                                                      >  2:1315681/1‑148 (MQ=255)
                                                     tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAAc                                                                      <  1:1315681/148‑1 (MQ=255)
                                                      ggggAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAAt                                                                                                                                               <  3:673752/74‑1 (MQ=255)
                                                      ggggAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAAt                                                                                                                                               >  4:673752/1‑74 (MQ=255)
                                                         gagaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGAccc                                                                                                                                      >  3:1416627/1‑80 (MQ=255)
                                                         gagaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGAccc                                                                                                                                      <  4:1416627/80‑1 (MQ=255)
                                                         gagaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCa                                                                                                      <  1:846932/112‑1 (MQ=255)
                                                         gagaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCa                                                                                                      >  2:846932/1‑112 (MQ=255)
                                                              aCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCt                                                                                                                      <  2:1354869/91‑1 (MQ=255)
                                                              aCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCt                                                                                                                      >  1:1354869/1‑91 (MQ=255)
                                                              aCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                                                              <  1:844496/99‑1 (MQ=255)
                                                              aCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                                                              >  2:844496/1‑99 (MQ=255)
                                                              aCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGc                                                            >  1:138598/1‑149 (MQ=255)
                                                                aGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATc                                                                                                                           <  2:185866/84‑1 (MQ=255)
                                                                aGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATc                                                                                                                           >  1:185866/1‑84 (MQ=255)
                                                                    ggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGc                                                                                                                              >  1:795084/1‑77 (MQ=255)
                                                                    ggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGc                                                                                                                              <  2:795084/77‑1 (MQ=255)
                                                                         tACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgacga                                                 <  2:724908/149‑1 (MQ=255)
                                                                               tttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGt                                                                                                             <  1:1187870/83‑1 (MQ=255)
                                                                               tttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGt                                                                                                             >  2:1187870/1‑83 (MQ=255)
                                                                               tttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaa                                                                           >  3:1463084/1‑117 (MQ=255)
                                                                               tttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaa                                                                           <  4:1463084/117‑1 (MQ=255)
                                                                                ttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTtat                                          <  2:138598/149‑1 (MQ=255)
                                                                                   gTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGt                                                                                                             >  4:519271/1‑79 (MQ=255)
                                                                                   gTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGt                                                                                                             <  3:519271/79‑1 (MQ=255)
                                                                                          cATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaat                                                                          <  1:203348/107‑1 (MQ=255)
                                                                                          cATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaat                                                                          >  2:203348/1‑107 (MQ=255)
                                                                                           aTTTCATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATAtcg                               >  4:741715/1‑149 (MQ=255)
                                                                                            ttttATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGttt                                                                                                           >  2:1168138/1‑72 (MQ=255)
                                                                                            ttttATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGttt                                                                                                           <  1:1168138/72‑1 (MQ=255)
                                                                                                  aCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgacga                                                 >  2:1382146/1‑124 (MQ=255)
                                                                                                  aCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgacga                                                 <  1:1382146/124‑1 (MQ=255)
                                                                                                     aTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGa                     <  1:588211/149‑1 (MQ=255)
                                                                                                            tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTtata                                         >  2:1021151/1‑122 (MQ=255)
                                                                                                            tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTtata                                         <  1:1021151/122‑1 (MQ=255)
                                                                                                            tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGt                        >  3:401960/1‑139 (MQ=255)
                                                                                                            tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGt                        <  4:401960/139‑1 (MQ=255)
                                                                                                                     ggtccTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGa                                    <  1:493744/115‑1 (MQ=255)
                                                                                                                     ggtccTGAATTTCAGGGCCCTTTCTGGCATAAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGa                                    >  2:493744/4‑118 (MQ=255)
                                                                                                                       tccTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCtta                                                                            <  3:1116873/75‑1 (MQ=39)
                                                                                                                       tccTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCtta                                                                            >  4:1116873/2‑76 (MQ=39)
                                                                                                                        ccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGcac  >  3:1207014/1‑149 (MQ=255)
                                                                                                                        ccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGcac  >  1:687902/1‑149 (MQ=255)
                                                                                                                              |                                                                                                                                              
TGATGGGTCTGGCTGGCTGATATCAACGGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATATTCGATATCGTCGTTGTCGAGGTTACGGTCGCCCAGCAC  >  NC_000913/216911‑217179

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