Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,642,084 G→A A254V (GCG→GTG)  rodZ ← MreB assembly cytoskeletal protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,642,0840GA97.9% 161.6 / ‑4.3 47A254V (GCG→GTG) rodZMreB assembly cytoskeletal protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (23/23);  total (23/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.67e-01

GGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGAAGGCGAAACAACCGCATTCTGTTGCG  >  NC_000913/2641943‑2642231
                                                                                                                                             |                                                                                                                                                   
ggCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGAt                                                                                                                                            >  1:479699/1‑151 (MQ=255)
          ttttCAGTTTGCACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCTGCACTTTCGGATCAGCCACCg                                                                                                                                   <  1:653041/150‑1 (MQ=255)
                       cGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCa                                                                                                                     <  2:410841/151‑1 (MQ=255)
                              tGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGAt                                                                                                              >  2:524537/1‑151 (MQ=255)
                                  cGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCgg                                                                                                           >  2:466994/1‑150 (MQ=255)
                                  cGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGt                                                                                                          >  2:135021/1‑151 (MQ=255)
                                   ggTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGtt                                                                                                         >  1:182134/1‑151 (MQ=255)
                                    gTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTg                                                                                                        <  1:283246/151‑1 (MQ=255)
                                    gTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTg                                                                                                        <  1:4187/151‑1 (MQ=255)
                                    gTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTg                                                                                                        >  2:315603/1‑151 (MQ=255)
                                     ttAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTgg                                                                                                       >  1:405397/1‑151 (MQ=255)
                                        aGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGt                                                                                                          >  2:739452/1‑145 (MQ=255)
                                                  aCCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCgcagc                                                                                           >  2:398819/1‑150 (MQ=255)
                                                   ccGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCAc                                                                                         <  1:510954/151‑1 (MQ=255)
                                                    cGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCAcc                                                                                        >  2:142218/1‑151 (MQ=255)
                                                      tCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCAt                                                                                      >  1:481938/1‑151 (MQ=255)
                                                         ttACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATc                                                                                     <  2:104032/149‑1 (MQ=255)
                                                                       gCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGcc                                                                     <  1:619396/151‑1 (MQ=255)
                                                                           aaCAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGt                                                                   >  1:599039/1‑149 (MQ=255)
                                                                              aaTTTTTTACCTGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGTTCACGCCAGCCTGATCTGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTg                                                              <  1:343818/151‑1 (MQ=255)
                                                                                ttttttACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGcc                                                            >  2:578225/1‑151 (MQ=255)
                                                                                ttttttACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGAGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTg                                                              >  1:701220/1‑149 (MQ=255)
                                                                                 tttttACCGGTAGCATCAGTGACCTCCAGCCAGCAATCCGCCGTAAAGTTCATCACCAGCACATTCGGCTCAGCCACCGGCGTGGTCACGCTAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTGGCTGCCGTTGGTg                                                              >  1:614544/1‑148 (MQ=255)
                                                                                  ttttACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCgg                                                          >  1:604571/1‑151 (MQ=255)
                                                                                      aCCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTc                                                      <  2:498628/151‑1 (MQ=255)
                                                                                       ccGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTc                                                      >  1:823151/1‑150 (MQ=255)
                                                                                        cGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTcgc                                                    <  1:787626/151‑1 (MQ=255)
                                                                                        cGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTcgc                                                    <  1:391519/151‑1 (MQ=255)
                                                                                            aGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCgcggcg                                                >  2:771632/1‑151 (MQ=255)
                                                                                            aGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCgcggc                                                 <  1:139746/150‑1 (MQ=255)
                                                                                                      aCCTCCAGCCATCAAGCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCTCCGGCGTGGTGACGCCAGCCTGATCGGTTGGCAAAGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGTGTCGCGGCGGTa                                             <  1:216392/144‑1 (MQ=255)
                                                                                                       ccTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACAt                                      <  2:512753/150‑1 (MQ=255)
                                                                                                          ccagccagCAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGc                                  <  2:205523/151‑1 (MQ=255)
                                                                                                           cagccagcAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGcc                                 >  2:621851/1‑151 (MQ=255)
                                                                                                                agcAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGa                            <  1:578187/151‑1 (MQ=255)
                                                                                                                 gcAATCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGaa                           <  2:820666/151‑1 (MQ=255)
                                                                                                                    aTCGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGAAGGc                        <  2:435367/151‑1 (MQ=255)
                                                                                                                      cGGCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGAAGGCGa                      <  1:608472/151‑1 (MQ=255)
                                                                                                                        gCAGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGAAGGCGaaa                    <  1:744535/151‑1 (MQ=255)
                                                                                                                          aGTAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGAAGGCGAaaca                  <  1:470208/151‑1 (MQ=255)
                                                                                                                            tAAAGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGAAGGCGAaacaac                >  2:118008/1‑151 (MQ=255)
                                                                                                                              aaGTTCATCACCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGAAGGCGAAACAAcc               >  1:518301/1‑150 (MQ=255)
                                                                                                                                gTTCATCACCAGCACATTCGGATCAGCCACCGGCGGGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGAAGGCGAAACAACCGCa            <  1:212519/151‑1 (MQ=255)
                                                                                                                                  tcatcaCCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGAAGGCGAAACAACCGCAtt          >  2:266652/1‑151 (MQ=255)
                                                                                                                                    atcaCCAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGAAGGCGAAACAACCGCATTCt        <  1:184192/151‑1 (MQ=255)
                                                                                                                                         cAGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTTCCTGCGAAGGCGAAACAACCGCATTCTGTTGc   >  1:227622/1‑151 (MQ=255)
                                                                                                                                          aGCACATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGAAGGCGAAACAACCGCATTCTGTTGCg  <  2:462126/151‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GGCGCACCAATTTTCAGTTTGTACGGTGCCTGGCCGGTTAAGTTCAAATTACCGTCTTTACGCTGCATACCGCTAAACAATTTTTTACCGGTAGCATCAGTGACCTCCAGCCAGCAATCGGCAGTAAAGTTCATCACCAGCGCATTCGGATCAGCCACCGGCGTGGTCACGCCAGCCTGATCGGTTGGCAACGGCGCAGCACCATCTGGCGTTGTTGCTGCCGTTGGTGCCGGGGTCGCGGCGGTATCAACATTTGCCTGCGAAGGCGAAACAACCGCATTCTGTTGCG  >  NC_000913/2641943‑2642231

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