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.96.9% 121.9 / ‑1.8 32coding (169/3633 nt)yehIDUF4132 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (15/16);  total (15/17)
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.51e-01

AAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGGCGGCG  >  NC_000913/2200334‑2200566
                                                                                                                 |                                                                                                                       
aaaaaGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATgag                                                                                               <  1:674520/139‑1 (MQ=255)
 aaaaGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATgaga                                                                                              >  1:567225/1‑139 (MQ=255)
                   cAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTgg                                                                             <  1:599448/139‑1 (MQ=255)
                                   ccGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCTAGACATTGGTTGATGAATTGTTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCa                                                            <  1:438316/139‑1 (MQ=255)
                                       gCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGcc                                                        >  1:25003/1‑139 (MQ=255)
                                       gCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGcc                                                        <  1:411051/139‑1 (MQ=255)
                                                     gTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCa                                                       <  1:113603/126‑1 (MQ=255)
                                                     gTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCa                                                       >  2:113603/1‑126 (MQ=255)
                                                             tGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAAcaca                                  <  2:567225/139‑1 (MQ=255)
                                                              ggAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCg                                        <  1:212975/132‑1 (MQ=255)
                                                              ggAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCg                                        >  2:212975/1‑132 (MQ=255)
                                                               gAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCa                                                            >  1:406324/1‑111 (MQ=255)
                                                               gAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCa                                                            <  2:406324/111‑1 (MQ=255)
                                                                    ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGt                                                                                                            >  2:110668/1‑58 (MQ=255)
                                                                    ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGt                                                                                                            <  1:110668/58‑1 (MQ=255)
                                                                    ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGc                                                                                                       <  2:979065/63‑1 (MQ=255)
                                                                    ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGc                                                                                                       >  1:979065/1‑63 (MQ=255)
                                                                    ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTc                           >  1:56955/1‑139 (MQ=255)
                                                                     tcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCa                          >  1:337193/1‑139 (MQ=255)
                                                                                 tGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGCGAATTTTTCACTGGAGTgg                                                                             >  2:751026/1‑76 (MQ=255)
                                                                                 tGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTgg                                                                             <  1:751026/76‑1 (MQ=255)
                                                                                       atgaATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATccc                     <  2:254889/126‑1 (MQ=255)
                                                                                       atgaATTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATccc                     >  1:254889/1‑126 (MQ=255)
                                                                                         aaaTTGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATggc      >  1:458552/2‑139 (MQ=255)
                                                                                             tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGAt                                             <  1:596395/96‑1 (MQ=255)
                                                                                             tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGAt                                             >  2:596395/1‑96 (MQ=255)
                                                                                             tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  >  2:955927/1‑139 (MQ=255)
                                                                                             tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  2:337193/139‑1 (MQ=255)
                                                                                             tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  2:458552/139‑1 (MQ=255)
                                                                                             tGCTACAAAAAGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  2:56955/139‑1 (MQ=255)
                                                                                                       aGCTGAGAAGAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCt                            <  1:1070819/103‑1 (MQ=255)
                                                                                                       aGCTGAGAAGAAAAAAACATTGTCTGGAGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCt                            >  2:1070819/1‑103 (MQ=255)
                                                                                                            agaagaAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  2:354221/124‑1 (MQ=255)
                                                                                                            agaagaAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  >  1:354221/1‑124 (MQ=255)
                                                                                                                 |                                                                                                                       
AAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGGCGGCG  >  NC_000913/2200334‑2200566

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