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.97.5% 155.8 / ‑2.0 40coding (169/3633 nt)yehIDUF4132 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (21/18);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

AATATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGGCGGCGC  >  NC_000913/2200328‑2200567
                                                                                                                       |                                                                                                                        
aaTATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGaaaaaaaca                                                                                                                  >  2:160518/1‑126 (MQ=255)
aaTATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGaaaaaaaca                                                                                                                  <  1:160518/128‑3 (MQ=255)
           ggCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGaaaaaaac                                                                                                                   >  2:195544/1‑115 (MQ=255)
           ggCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGaaaaaaac                                                                                                                   <  1:195544/116‑2 (MQ=255)
            gCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATgag                                                                                                <  1:1100344/133‑1 (MQ=255)
            gCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATgag                                                                                                >  2:1100344/1‑133 (MQ=255)
               aaaCGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTc                                                                                        <  2:86715/139‑1 (MQ=255)
               aaaCGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCa                                                                                       >  1:757036/1‑139 (MQ=255)
                            tCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGt                                                                          <  2:560525/139‑1 (MQ=255)
                                      tGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCAc                                                                >  1:458706/1‑139 (MQ=255)
                                       ggCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACg                                                               <  2:458706/139‑1 (MQ=255)
                                                            ttATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGc                                          >  1:17490/1‑139 (MQ=255)
                                                               tCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGaa                                       <  2:445595/139‑1 (MQ=255)
                                                                     gAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGaa                                                                                              >  2:327385/1‑78 (MQ=255)
                                                                     gAAGTCTCATCCAGACAGTTGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGaa                                                                                              <  1:327385/78‑1 (MQ=255)
                                                                          ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGaa                                       <  1:345274/128‑1 (MQ=255)
                                                                          ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGaa                                       >  2:345274/1‑128 (MQ=255)
                                                                           tcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGTCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCa                           >  1:353695/1‑139 (MQ=255)
                                                                           tcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCg                                         <  1:753421/125‑1 (MQ=255)
                                                                           tcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCg                                         >  2:753421/1‑125 (MQ=255)
                                                                             aTCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAAt                         >  2:719219/1‑139 (MQ=255)
                                                                             aTCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAAt                         >  2:372406/1‑139 (MQ=255)
                                                                              tCCAGACATTGGGTGATGAATTTCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATc                        <  1:1152055/139‑1 (MQ=255)
                                                                                 aGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATcccg                     <  1:372406/139‑1 (MQ=255)
                                                                                   aCATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCt                                                                                                           <  2:289106/51‑1 (MQ=255)
                                                                                   aCATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCt                                                                                                           >  1:289106/1‑51 (MQ=255)
                                                                                   aCATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATcccggc                   >  2:485575/1‑139 (MQ=255)
                                                                                       tGGGTGATGAATTGCTACAAAAAGCTGAGAAGaaaaaaac                                                                                                                   >  1:642449/1‑39 (MQ=255)
                                                                                       tGGGTGATGAATTGCTACAAAAAGCTGAGAAGaaaaaaac                                                                                                                   <  2:642449/40‑2 (MQ=255)
                                                                                       tGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATcccggcccgg               >  2:1060133/1‑139 (MQ=255)
                                                                                         ggTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGAc             >  1:973388/1‑139 (MQ=255)
                                                                                         ggTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGAc             >  1:1117512/1‑139 (MQ=255)
                                                                                            gatgaATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAAt          >  1:1113321/1‑139 (MQ=255)
                                                                                              tgaATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTgg                                                                              <  2:1182855/69‑1 (MQ=255)
                                                                                              tgaATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTgg                                                                              >  1:1182855/1‑69 (MQ=255)
                                                                                              tgaATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATgg        <  2:973388/139‑1 (MQ=255)
                                                                                                   ttcTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg   <  2:140758/137‑1 (MQ=255)
                                                                                                   tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAAtt                                                                                            <  2:180433/50‑1 (MQ=255)
                                                                                                   tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAAtt                                                                                            >  1:180433/1‑50 (MQ=255)
                                                                                                   tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGTCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg   <  2:353695/139‑1 (MQ=255)
                                                                                                   tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg   <  2:1117512/139‑1 (MQ=255)
                                                                                                   tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg   >  1:880761/1‑139 (MQ=255)
                                                                                                   tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg   <  1:485575/139‑1 (MQ=255)
                                                                                                    tgctaaaaaaGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCTCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGGCGgcgc  >  1:709986/6‑139 (MQ=255)
                                                                                                       aCAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGa                                               >  2:1115678/1‑91 (MQ=255)
                                                                                                       aCAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGa                                               <  1:1115678/91‑1 (MQ=255)
                                                                                                                       |                                                                                                                        
AATATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGGCGGCGC  >  NC_000913/2200328‑2200567

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: