Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,625,249 G→T 17.8% 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,2490GT17.8% 61.2 / 4.7 27intergenic (+30/+27)PP_1425/nadBhypothetical protein/L‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base G (11/11);  new base T (2/3);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGCCA  >  NC_002947/1625105‑1625394
                                                                                                                                                |                                                                                                                                                  
gCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGccc                                                                                                                                               >  2:588041/1‑150 (MQ=255)
         aaGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGcc                                                                                                                                      <  2:9066/150‑1 (MQ=255)
                    taCCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGat                                                                                                                           >  1:21314/1‑150 (MQ=255)
                    taCCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGat                                                                                                                           <  2:731944/150‑1 (MQ=255)
                       cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCTCGCCCGCGAAgaa                                                                                                                                       <  2:694048/135‑3 (MQ=255)
                       cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCTCGCCCGCGAAgaa                                                                                                                                       >  1:694048/1‑133 (MQ=255)
                                                    cGGCGCTGGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCAAAGATGGAGCCCGGGTTAGCGTCTTCGCGGCCGAGCCCGCGAAGGCCACGACGCTGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGcc                                                                                          <  1:676003/150‑1 (MQ=255)
                                                           tGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCa                                                                                   <  2:21314/150‑1 (MQ=255)
                                                               ggTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCAt                                                                               >  2:263967/1‑150 (MQ=255)
                                                                       tGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATg                                                                                                       <  1:117911/119‑1 (MQ=255)
                                                                       tGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATg                                                                                                       >  2:117911/1‑119 (MQ=255)
                                                                          aGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAt                                                                    >  2:374633/1‑150 (MQ=255)
                                                                            ttCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATcc                                                                  >  2:244077/1‑150 (MQ=255)
                                                                                   ggcggcGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgtgtgt                                                           >  1:744529/1‑150 (MQ=255)
                                                                                      ggcgCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCTCGCCCGCGAAgaa                                                                                                                                       <  2:333782/72‑3 (MQ=37)
                                                                                      ggcgCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCTCGCCCGCGAAgaa                                                                                                                                       >  1:333782/1‑70 (MQ=37)
                                                                                       gcgcTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCACCATCCCTGGATAATCCGGTGTGTaagg                                                       <  2:421980/150‑3 (MQ=255)
                                                                                          cTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTg                                                                                             <  1:410234/109‑1 (MQ=255)
                                                                                          cTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTg                                                                                             >  2:410234/1‑109 (MQ=255)
                                                                                                tGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCg                                              <  1:263967/150‑1 (MQ=255)
                                                                                                             gCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTgcgc                                 <  1:244077/150‑1 (MQ=255)
                                                                                                                      gagaCCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCt                        <  2:744529/150‑1 (MQ=255)
                                                                                                                         aCCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGAc                     >  1:307734/1‑150 (MQ=255)
                                                                                                                             ggTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGa                 <  1:588041/150‑1 (MQ=255)
                                                                                                                                 gtctccttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTTCTACGGGATCAAACCGTGGTCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGatcat             <  2:401676/145‑1 (MQ=255)
                                                                                                                                      cttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCt        >  1:356768/1‑150 (MQ=255)
                                                                                                                                            gggCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGCCa  >  1:484189/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGCCA  >  NC_002947/1625105‑1625394

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