Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,200,454 (A)8→7 coding (176/3633 nt) yehI → DUF4132 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,200,4470A.93.5% 107.8 / 2.9 31coding (169/3633 nt)yehIDUF4132 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (14/15);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

ACTGGAACTGAAATATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGGCGGCG  >  NC_000913/2200317‑2200566
                                                                                                                                   |                                                                                                                       
aCTGGAACTGAAATATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAAGCTGAGAAGaaaaaaac                                                                                                                  >  1:1141300/1‑138 (MQ=255)
      aCTGAAATATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCt                                                                                                          <  2:1464431/139‑1 (MQ=255)
           aaTATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGt                                                                                                     >  1:519578/1‑139 (MQ=255)
                               ccgcTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTgg                                                                                  <  2:1141300/139‑1 (MQ=255)
                                  cTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTg                                                                              >  1:726121/1‑139 (MQ=255)
                                   tCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTgg                                                                             <  2:519578/139‑1 (MQ=255)
                                        cATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTc                                                                        <  2:650434/139‑1 (MQ=255)
                                             tcgCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCa                                                                   <  2:726121/139‑1 (MQ=255)
                                                 tGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCAc                                                               >  2:1246801/1‑139 (MQ=255)
                                                        gCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGcc                                                        <  2:846111/139‑1 (MQ=255)
                                                                               ggAAGTCTCATCCAGACATTGGGTGATGAATTACTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACaa                                 >  2:911275/1‑139 (MQ=255)
                                                                                gAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTc                                                                                       >  1:169577/1‑84 (MQ=255)
                                                                                gAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTc                                                                                       <  2:169577/84‑1 (MQ=255)
                                                                                    tctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCt                            <  1:1234937/139‑1 (MQ=255)
                                                                                        aTCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGACAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCaa                         >  2:485458/1‑139 (MQ=255)
                                                                                            aGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATcccg                    <  1:154731/139‑1 (MQ=255)
                                                                                              aCATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCt                                                                                                          <  2:728255/51‑1 (MQ=255)
                                                                                              aCATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCt                                                                                                          >  1:728255/1‑51 (MQ=255)
                                                                                                  tGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCg                                                                      >  2:939106/1‑83 (MQ=255)
                                                                                                  tGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCg                                                                      <  1:939106/83‑1 (MQ=255)
                                                                                                  tGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATcccggcccgg              >  1:626028/1‑139 (MQ=255)
                                                                                                    ggTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGAc            >  1:795429/1‑139 (MQ=255)
                                                                                                       gatgaATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAAt         >  1:629029/1‑139 (MQ=255)
                                                                                                        atgaATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCAGGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATg        >  2:1194674/1‑139 (MQ=255)
                                                                                                        atgaATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAgg                                                          <  2:241604/89‑1 (MQ=255)
                                                                                                        atgaATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAgg                                                          >  1:241604/1‑89 (MQ=255)
                                                                                                             ttGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATggcggc   >  1:475023/1‑139 (MQ=255)
                                                                                                              tGCTACAACAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCTCGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  1:911275/139‑1 (MQ=255)
                                                                                                              tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  2:649311/139‑1 (MQ=255)
                                                                                                               gCTACAAAAAGCTGACAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  1:485458/139‑1 (MQ=255)
                                                                                                                         aGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATcccggcc                 >  2:1284515/1‑114 (MQ=255)
                                                                                                                         aGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATcccggcc                 <  1:1284515/114‑1 (MQ=255)
                                                                                                                                   |                                                                                                                       
ACTGGAACTGAAATATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGGCGGCG  >  NC_000913/2200317‑2200566

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