Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,424,543 Δ1 bp coding (249/276 nt) yjfY ← YhcN family protein, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,424,5430C.98.3% 231.7 / ‑2.5 60coding (249/276 nt)yjfYYhcN family protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (27/32);  minor base A (0/1);  total (27/33)
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

ATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGA  >  NC_000913/4424420‑4424660
                                                                                                                           |                                                                                                                     
aTAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTcc                                                                                                       <  2:1072243/139‑1 (MQ=255)
          tAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCg                                                                                             <  1:950969/139‑1 (MQ=255)
                      cGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTAc                                                                                 >  2:1852874/1‑139 (MQ=255)
                       gTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                <  2:905498/139‑1 (MQ=255)
                        tttttCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGGGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGt                                                                               >  1:55023/1‑139 (MQ=255)
                         ttttCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGTTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTg                                                                              <  1:1104086/139‑1 (MQ=255)
                                        ccccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGcc                                                               <  1:1019496/139‑1 (MQ=255)
                                        ccccTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGcc                                                               <  1:1852874/139‑1 (MQ=255)
                                           cTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTg                                                            <  1:360532/139‑1 (MQ=255)
                                             gCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgcg                                                          <  2:55023/139‑1 (MQ=255)
                                                        aaagcaaaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTAATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttct                                               <  1:337576/139‑1 (MQ=255)
                                                              aaCTATGGCAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCAttt                                         <  2:976119/139‑1 (MQ=255)
                                                                      cAATGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGc                                    <  2:1322646/136‑1 (MQ=255)
                                                                      cAATGTTACGTTGGTTATACCTGGCGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGCGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGc                                    >  1:1322646/1‑136 (MQ=255)
                                                                        aTGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGc                                    >  2:539821/1‑134 (MQ=255)
                                                                        aTGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGc                                    <  1:539821/134‑1 (MQ=255)
                                                                        aTGTTTCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc                               <  1:427053/139‑1 (MQ=255)
                                                                           tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCAt                                                                                          >  2:1591123/1‑77 (MQ=255)
                                                                           tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCAt                                                                                          <  1:1591123/77‑1 (MQ=255)
                                                                           tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                <  1:289870/87‑1 (MQ=255)
                                                                           tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                >  2:289870/1‑87 (MQ=255)
                                                                           tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAAt                                      <  2:1460274/129‑1 (MQ=255)
                                                                           tttCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAAt                                      >  1:1460274/1‑129 (MQ=255)
                                                                             tCGTTGGTTATACCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCt                          <  1:661315/139‑1 (MQ=255)
                                                                                     tataCCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCt                                                     >  2:1813365/1‑104 (MQ=255)
                                                                                     tataCCTGGTGTTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCt                                                     <  1:1813365/104‑1 (MQ=255)
                                                                                             gtgttAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgc                                                           >  2:632870/5‑90 (MQ=255)
                                                                                             gtgttAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGc                                    >  2:44461/5‑113 (MQ=255)
                                                                                             gtgttAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGc                                    <  1:44461/109‑1 (MQ=255)
                                                                                             gtgttAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCg           >  1:742810/5‑138 (MQ=255)
                                                                                             gtgttAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGa          >  2:262106/5‑139 (MQ=255)
                                                                                             gtgttAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGa          >  2:1597062/5‑139 (MQ=255)
                                                                                             gtgtTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCt                                                              <  1:865065/87‑1 (MQ=255)
                                                                                             gtgtTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCt                                                              >  2:865065/1‑87 (MQ=255)
                                                                                             gtgtTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgc                                                           <  1:632870/90‑1 (MQ=255)
                                                                                             gtgtTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCg           <  2:742810/138‑1 (MQ=255)
                                                                                             gtgtTAGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGa          >  2:1350195/1‑139 (MQ=255)
                                                                                                  aGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTcc                                                                                               >  2:441762/1‑49 (MQ=39)
                                                                                                  aGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTcc                                                                                               <  1:441762/49‑1 (MQ=39)
                                                                                                  aGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                <  2:638134/64‑1 (MQ=255)
                                                                                                  aGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACg                                                                                >  1:638134/1‑64 (MQ=255)
                                                                                                  aGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTgcg                                                                 <  1:65533/79‑1 (MQ=255)
                                                                                                  aGCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTgcg                                                                 >  2:65533/1‑79 (MQ=255)
                                                                                                   gCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc                               >  2:779814/1‑112 (MQ=255)
                                                                                                   gCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGAc                               <  1:779814/112‑1 (MQ=255)
                                                                                                   gCGATAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAt    >  2:1214036/1‑139 (MQ=255)
                                                                                                       tAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCg                                                                                                    <  2:1611914/39‑1 (MQ=37)
                                                                                                       tAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCg                                                                                                    >  1:1611914/1‑39 (MQ=37)
                                                                                                       tAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgc                                                           <  2:1672848/80‑1 (MQ=255)
                                                                                                       tAAATATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgc                                                           >  1:1672848/1‑80 (MQ=255)
                                                                                                        aaaTATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGcc                                                               <  2:469994/75‑1 (MQ=255)
                                                                                                        aaaTATCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGcc                                                               >  1:469994/1‑75 (MQ=255)
                                                                                                             tCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTgc                                                           >  2:1042649/1‑74 (MQ=255)
                                                                                                             tCCGCGCTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGCTAGTACGTTGCGCCTTgc                                                           <  1:1042649/74‑1 (MQ=255)
                                                                                                                 cgcTGGCGTGATATTGCCGTTACTCCCCGGTTccc                                                                                              >  1:898248/1‑35 (MQ=37)
                                                                                                                 cgcTGGCGTGATATTGCCGTTACTCCCCGGTTccc                                                                                              <  2:898248/35‑1 (MQ=37)
                                                                                                                 cgcTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtga  >  2:62565/1‑127 (MQ=255)
                                                                                                                 cgcTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTAtga  <  1:62565/127‑1 (MQ=255)
                                                                                                                   cTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttctt                                              <  2:34215/81‑1 (MQ=255)
                                                                                                                   cTGGCGTGATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTttctt                                              >  1:34215/1‑81 (MQ=255)
                                                                                                                           |                                                                                                                     
ATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAATGCCTGACGTGCTTCGCTTTCAGTGGCGAAATTATGA  >  NC_000913/4424420‑4424660

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: