Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,076,654 C→T A170V (GCG→GTG)  yihY → putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,076,6540CT98.1% 182.3 / ‑5.6 53A170V (GCG→GTG) yihYputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (28/24);  total (29/24)
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.82e-01

GATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGC  >  NC_000913/4076577‑4076731
                                                                             |                                                                             
gatgatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgtgag                                                                            >  1:822725/1‑81 (MQ=255)
 atgatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgtgagc                                                                           <  2:822651/81‑1 (MQ=255)
  tgatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgtgagcg                                                                          <  1:822653/81‑1 (MQ=255)
  tgatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgtgagcg                                                                          <  1:822652/81‑1 (MQ=255)
   gatTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgtgagcga                                                                         >  2:822738/1‑81 (MQ=255)
    atTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGgtgagcga                                                                         <  1:822655/80‑1 (MQ=255)
    atTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGAt                                                                        <  2:822654/81‑1 (MQ=255)
    atTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGAt                                                                        <  2:822656/81‑1 (MQ=255)
     ttttAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGAtc                                                                       <  2:822657/81‑1 (MQ=255)
       ttAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGAtct                                                                      <  2:822658/80‑1 (MQ=255)
         aaCGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCaa                                                                   >  1:822749/1‑81 (MQ=255)
         aaCGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCa                                                                    <  2:822659/80‑1 (MQ=255)
          aCGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCaa                                                                   >  1:822735/1‑80 (MQ=255)
          aCGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAAt                                                                  <  1:822660/81‑1 (MQ=255)
            gCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATAc                                                                >  1:822741/1‑81 (MQ=255)
             cTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACt                                                               <  1:822661/81‑1 (MQ=255)
              tGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTg                                                              >  2:822758/1‑81 (MQ=255)
               ggggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGt                                                             >  1:822732/1‑81 (MQ=255)
                gggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTc                                                            <  2:822662/81‑1 (MQ=255)
                 ggCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCa                                                           <  2:822663/81‑1 (MQ=255)
                   ccGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATc                                                         <  2:822664/81‑1 (MQ=255)
                    cGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCg                                                        <  1:822665/81‑1 (MQ=255)
                     gctgctGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGa                                                       >  1:822747/1‑81 (MQ=255)
                        gctgGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACaa                                                    >  2:822751/1‑81 (MQ=255)
                         ctgGCAGGGGCTAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAAc                                                   <  1:822669/81‑1 (MQ=255)
                         ctgGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACaa                                                    <  1:822666/80‑1 (MQ=255)
                         ctgGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACa                                                     <  1:822668/79‑1 (MQ=255)
                         ctgGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACa                                                     <  1:822667/79‑1 (MQ=255)
                          tgGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACg                                                  >  1:822743/1‑81 (MQ=255)
                           ggCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAAc                                                   >  1:822777/1‑79 (MQ=255)
                           ggCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGt                                                 >  1:822713/1‑81 (MQ=255)
                             cAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGtgc                                               <  1:822670/81‑1 (MQ=255)
                              aGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGtgct                                              >  1:822757/1‑81 (MQ=255)
                              aGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGtgct                                              <  2:822671/81‑1 (MQ=255)
                                gggCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGtgctgc                                            <  2:822672/81‑1 (MQ=255)
                                gggCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGCCATCGACAACGtgctgc                                            >  2:822746/1‑81 (MQ=255)
                                   ccAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGTTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTa                                         >  1:822745/1‑81 (MQ=255)
                                   ccAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTa                                         >  2:822792/1‑81 (MQ=255)
                                      gTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAttt                                      >  2:822773/1‑81 (MQ=255)
                                        cTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAttttt                                    >  2:822770/1‑81 (MQ=255)
                                        cTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTAttttt                                    >  1:822748/1‑81 (MQ=255)
                                             gATCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGtt                               >  1:822767/1‑81 (MQ=255)
                                              aTCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTg                              <  1:822673/81‑1 (MQ=255)
                                               tCAGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGc                             >  2:822762/1‑81 (MQ=255)
                                                 aGTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtg                           >  1:822740/1‑81 (MQ=255)
                                                  gTTCCTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgt                          >  2:822744/1‑81 (MQ=255)
                                                   ttCCTATTTGCTCTCCCTGCGCTGGGTGAACGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCtgtt                         <  1:822674/81‑1 (MQ=255)
                                                      cTATTTGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTc                      >  2:822768/1‑81 (MQ=255)
                                                           tGCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGa                 >  2:822737/1‑81 (MQ=255)
                                                            gCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGa                 <  2:822675/80‑1 (MQ=255)
                                                            gCTCTCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGAt                >  2:822731/1‑81 (MQ=255)
                                                             ctctCCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGAtc               <  2:822676/81‑1 (MQ=255)
                                                                tcCCTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTcc            >  1:822756/1‑81 (MQ=255)
                                                                   cTGCGCTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTc         <  2:822677/81‑1 (MQ=255)
                                                                     gcgcTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTCTTCCGTTGCTGTTGTCGTGGATCTCCTTCTg       >  2:822763/1‑81 (MQ=255)
                                                                      cgcTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCt        <  2:822678/79‑1 (MQ=255)
                                                                      cgcTGGGTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTgg      >  2:822753/1‑81 (MQ=255)
                                                                          gggTGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGc  >  2:822766/1‑81 (MQ=255)
                                                                          gggCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGc  >  2:822742/1‑81 (MQ=255)
                                                                             |                                                                             
GATGATTTTAACGCTGGGGCCGCTGCTGGCAGGGGCCAGTCTGGCGATCAGTTCCTATTTGCTCTCCCTGCGCTGGGCGAGCGATCTCAATACTGTCATCGACAACGTGCTGCGTATTTTTCCGTTGCTGTTGTCGTGGATCTCCTTCTGGTTGC  >  NC_000913/4076577‑4076731

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