Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,625,207 G→T 18.0% intergenic (+30/+27) PP_1425 → / ← nadB hypothetical protein/L‑aspartate oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,625,2070GT18.0% 50.2 / 5.9 27intergenic (+30/+27)PP_1425/nadBhypothetical protein/L‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base G (14/8);  new base T (2/3);  total (16/11)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

CTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCG  >  NC_002947/1625086‑1625348
                                                                                                                         |                                                                                                                                             
cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCa                                                                                                                                >  5:65360/1‑137 (MQ=255)
cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCa                                                                                                                                <  6:65360/137‑1 (MQ=255)
     gATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCagag                                                                                                               >  8:181318/1‑149 (MQ=255)
           cTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGt                                                                                                          >  3:51956/1‑148 (MQ=255)
            tAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGtt                                                                                                         >  7:76523/1‑148 (MQ=255)
              ggggCGGACGTGCTGCGGCGCTTGGGGGGGCGTTTGAATTTCAGCCGGCGGCGCTGGGCGGCGGTGAGCCGGGCATCGATCGAGCCCGGGTTAGCTTCTTCGGGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTaa                                                                                                       <  1:234250/148‑1 (MQ=255)
                  cTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATg                                                                                                  <  7:181318/149‑1 (MQ=255)
                       gAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAg                                                                                                                                    <  3:178041/110‑1 (MQ=255)
                       gAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAg                                                                                                                                    >  4:178041/1‑110 (MQ=255)
                                              gTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTcc                                                                                >  1:224265/1‑139 (MQ=255)
                                              gTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTcc                                                                                <  2:224265/139‑1 (MQ=255)
                                                          ccGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgtg                                                          >  3:93212/1‑149 (MQ=255)
                                                           cggcggCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGCCCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgtgt                                                         <  8:76523/149‑1 (MQ=255)
                                                                 cgcTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATg                                                   >  6:227931/1‑149 (MQ=255)
                                                                           cGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCgg                                         <  4:93212/149‑1 (MQ=255)
                                                                                 gCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAg                                                            >  6:243582/1‑124 (MQ=255)
                                                                                              cGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATgg                      >  5:88818/1‑149 (MQ=255)
                                                                                                 gaCCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCt                    >  7:251538/1‑148 (MQ=255)
                                                                                                  aCCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGa                  >  7:12535/1‑149 (MQ=255)
                                                                                                          agattattcGGGGGCTAGCCGGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGatc           <  2:128234/142‑1 (MQ=255)
                                                                                                           gCTTCTGCCCGGGCGCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGatca          <  8:251538/148‑1 (MQ=255)
                                                                                                            cttcttCGCGGTCGAGCCCGAGAAGGCCACGACGCAGATATCAGAGCGGCTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGatcat         >  2:13109/1‑148 (MQ=255)
                                                                                                               cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATc                                                               <  7:204587/91‑1 (MQ=255)
                                                                                                               cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATc                                                               >  8:204587/1‑91 (MQ=255)
                                                                                                                ttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCt    >  6:38877/1‑149 (MQ=255)
                                                                                                                  cgcgGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCg  <  5:38877/149‑1 (MQ=255)
                                                                                                                  cgcgGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCg  >  6:159570/1‑149 (MQ=255)
                                                                                                                         |                                                                                                                                             
CTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCG  >  NC_002947/1625086‑1625348

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: