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
                                                                                                                 |                                                                                                                       
aaaaaGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATgag                                                                                               <  1:674520/139‑1 (MQ=255)
 aaaaGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATgaga                                                                                              >  1:567225/1‑139 (MQ=255)
                   cAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTgg                                                                             <  1:599448/139‑1 (MQ=255)
                                   ccGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCTAGACATTGGTTGATGAATTGTTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCa                                                            <  1:438316/139‑1 (MQ=255)
                                       gCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGcc                                                        >  1:25003/1‑139 (MQ=255)
                                       gCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGcc                                                        <  1:411051/139‑1 (MQ=255)
                                                     gTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCa                                                       <  1:113603/126‑1 (MQ=255)
                                                     gTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCa                                                       >  2:113603/1‑126 (MQ=255)
                                                             tGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAAcaca                                  <  2:567225/139‑1 (MQ=255)
                                                              ggAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCg                                        <  1:212975/132‑1 (MQ=255)
                                                              ggAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCg                                        >  2:212975/1‑132 (MQ=255)
                                                               gAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCa                                                            >  1:406324/1‑111 (MQ=255)
                                                               gAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCa                                                            <  2:406324/111‑1 (MQ=255)
                                                                    ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGt                                                                                                            >  2:110668/1‑58 (MQ=255)
                                                                    ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGt                                                                                                            <  1:110668/58‑1 (MQ=255)
                                                                    ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGc                                                                                                       <  2:979065/63‑1 (MQ=255)
                                                                    ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGc                                                                                                       >  1:979065/1‑63 (MQ=255)
                                                                    ctcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTc                           >  1:56955/1‑139 (MQ=255)
                                                                     tcATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCa                          >  1:337193/1‑139 (MQ=255)
                                                                                 tGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGCGAATTTTTCACTGGAGTgg                                                                             >  2:751026/1‑76 (MQ=255)
                                                                                 tGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTgg                                                                             <  1:751026/76‑1 (MQ=255)
                                                                                       atgaATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATccc                     <  2:254889/126‑1 (MQ=255)
                                                                                       atgaATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATccc                     >  1:254889/1‑126 (MQ=255)
                                                                                         aaaTTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATggc      >  1:458552/2‑139 (MQ=255)
                                                                                             tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGAt                                             <  1:596395/96‑1 (MQ=255)
                                                                                             tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGAt                                             >  2:596395/1‑96 (MQ=255)
                                                                                             tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  >  2:955927/1‑139 (MQ=255)
                                                                                             tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  2:337193/139‑1 (MQ=255)
                                                                                             tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  2:458552/139‑1 (MQ=255)
                                                                                             tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  2:56955/139‑1 (MQ=255)
                                                                                                       aGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCt                            <  1:1070819/103‑1 (MQ=255)
                                                                                                       aGCTGAGAAG‑AAAAAAACATTGTCTGGAGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCt                            >  2:1070819/1‑103 (MQ=255)
                                                                                                            agaag‑aAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  <  2:354221/124‑1 (MQ=255)
                                                                                                            agaag‑aAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg  >  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: ‑