Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,625,208 A→C 25.4% intergenic (+31/+26) 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,2080AC25.4% 53.5 / 24.0 43intergenic (+31/+26)PP_1425/nadBhypothetical protein/L‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base A (20/12);  new base C (5/6);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-01

GGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGCCACTT  >  NC_002947/1625068‑1625355
                                                                                                                                            |                                                                                                                                                    
gggAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCTCGCCCGCGa                                                                                                                                              <  3:280802/149‑1 (MQ=255)
 ggAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCTCGCCCGCGaa                                                                                                                                             <  1:246725/149‑1 (MQ=255)
  gAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCTCGCCCGCGAAg                                                                                                                                            <  8:224365/149‑1 (MQ=255)
                  cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGatat                                                                                                                            >  5:1350/1‑149 (MQ=255)
                  cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGata                                                                                                                             >  7:232677/1‑148 (MQ=255)
                  cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCAGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCcgc                                                                                                                                                <  6:299077/129‑1 (MQ=255)
                  cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCAGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCcgc                                                                                                                                                >  5:299077/1‑129 (MQ=255)
                           cgcTAGGGGCTGACGAGCTTCGGCGCTTGCGGGTGGGGTTGAAGTTCTGCCGGCGGCGCTGGGCTGGGGTGAGCCGGGAATAGATCGAGACCGGTTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCgg                                                                                                                    <  6:1350/148‑1 (MQ=255)
                                 gggCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAg                                                                                                              >  4:51743/1‑148 (MQ=255)
                                       aCGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGgtg                                                                                                       >  4:350938/1‑149 (MQ=255)
                                               cGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCg                                                                                                                     >  5:337376/1‑127 (MQ=255)
                                                           gTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCAt                                                                                   <  3:51743/149‑1 (MQ=255)
                                                               ggTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTc                                                                                         >  7:191025/1‑139 (MQ=255)
                                                                gTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTg                                                                              >  5:194333/1‑149 (MQ=255)
                                                                      gTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATTTCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAt                                                                        >  1:306802/1‑149 (MQ=255)
                                                                       ttCAGCCGGCGGCGCTGGGCTGCGGTCAGCCGGGCATAGATCGAGACCGGGATAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCGGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑c                                                                      >  2:332528/1‑149 (MQ=255)
                                                                            ccGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGCCCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAgtg                                                                 <  8:349731/149‑1 (MQ=255)
                                                                                gggcgCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTa                                                             <  4:252227/148‑1 (MQ=255)
                                                                                 ggcgCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTaa                                                            >  1:19909/1‑149 (MQ=255)
                                                                                            gCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCg                                                 >  7:263156/1‑149 (MQ=255)
                                                                                             cgggtcgCCGTGAATAGATCGATACCGGTTTTGCTTCTTCGGGGGCGACCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAATGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCgg                                                <  2:19909/143‑1 (MQ=255)
                                                                                                       gtcATAGATCGCGACCGGTTTAGCTTCTTCGCGGGCGCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTgc                                      <  7:270447/147‑1 (MQ=255)
                                                                                                          aTAGTTCGATCCCGGGTTAGCTTCTTCGCGGGCGAGCCCGGGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCt                                   <  6:140580/149‑1 (MQ=255)
                                                                                                           tAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTg                                  >  2:67035/1‑149 (MQ=255)
                                                                                                            aGATCGAGACCGGGTTAGCTTCTTCGCGGGCGCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTg                                  >  8:142823/1‑148 (MQ=255)
                                                                                                             gATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGATCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGc                                 <  1:67035/149‑1 (MQ=255)
                                                                                                               tCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCg                                                 >  2:260860/1‑130 (MQ=255)
                                                                                                               tCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCg                                                 <  1:260860/130‑1 (MQ=255)
                                                                                                                      cGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGa                         >  6:219649/1‑147 (MQ=255)
                                                                                                                               ttcttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGatcat                >  5:373253/1‑147 (MQ=255)
                                                                                                                                tcttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAg             >  1:369631/1‑149 (MQ=255)
                                                                                                                                 cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATgg                             <  6:219325/132‑1 (MQ=255)
                                                                                                                                 cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGatcatc               <  3:156667/146‑1 (MQ=255)
                                                                                                                                 cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGatcatc               >  4:156667/1‑146 (MQ=255)
                                                                                                                                 cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGc            >  6:164624/1‑149 (MQ=255)
                                                                                                                                 cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTGGCGCTGCATgg                             >  5:219325/1‑132 (MQ=255)
                                                                                                                                  ttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCt           >  8:248/1‑149 (MQ=255)
                                                                                                                                  ttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCt           >  1:166344/1‑149 (MQ=255)
                                                                                                                                   tcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTc          >  3:250004/1‑149 (MQ=255)
                                                                                                                                     gcgGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCg         <  8:263156/148‑1 (MQ=255)
                                                                                                                                      cgGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGc       <  6:194333/149‑1 (MQ=255)
                                                                                                                                         ccGACCCGGCCAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGCCAc    <  2:166344/148‑1 (MQ=255)
                                                                                                                                           gAGCCAGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTAGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGCCACtt  >  1:64495/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAAT‑CCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGCCACTT  >  NC_002947/1625068‑1625355

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: ‑