Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,625,208 A→C 13.3% 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,2080AC13.3% 98.9 / 2.5 42intergenic (+31/+26)PP_1425/nadBhypothetical protein/L‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base A (25/11);  new base C (2/4);  total (27/15)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCT  >  NC_002947/1625060‑1625346
                                                                                                                                                     |                                                                                                                                          
ttGGCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGa                                                                                                                                            <  6:313890/149‑1 (MQ=255)
          gAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAg                                                                                                                                  <  4:323616/149‑1 (MQ=255)
           aaaGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCTCGCCCGCGAAga                                                                                                                                 <  3:106148/149‑2 (MQ=255)
           aaaGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCTCGCCCGCGAAg                                                                                                                                  >  4:106148/1‑148 (MQ=255)
             aGCCCGTGTATACCTCAGGATGCGCTCGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGcc                                                                                                                               <  2:135457/149‑1 (MQ=255)
               ccAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGTCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCac                                                                                                                             >  2:309176/1‑149 (MQ=255)
                     tataCCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGcc                                                                                                                       >  2:338570/1‑149 (MQ=255)
                        aCCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGa                                                                                                                     >  5:240181/1‑148 (MQ=255)
                        aCCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGa                                                                                                                     >  6:148274/1‑148 (MQ=255)
                        aCCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGa                                                                                                                     >  3:298094/1‑148 (MQ=255)
                          cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGatat                                                                                                                  >  5:193333/1‑149 (MQ=255)
                           acAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTAGCGCGCGAGCCAGCGAAGGCCACGACGCCGatat                                                                                                                  >  7:14065/2‑148 (MQ=255)
                                   cgcTAGGGGCTGACGAGCAGCGGCGCGTGCGGGGGCGGTTGAAGTTCAGCGGGGGGGGCTGGGCTGCGGTGTGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCgg                                                                                                          <  4:298094/148‑1 (MQ=255)
                                      tAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACCCCGATATCAGAGCGGGtt                                                                                                       >  5:138305/1‑148 (MQ=255)
                                       aGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGtt                                                                                                       >  2:352833/1‑147 (MQ=255)
                                            cTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATg                                                                                                <  1:269905/149‑1 (MQ=255)
                                             tGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATgg                                                                                               >  7:1785/1‑149 (MQ=255)
                                                     aTCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTg                                                                                       <  6:193333/149‑1 (MQ=255)
                                                      tCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTgg                                                                                      >  4:367155/1‑149 (MQ=255)
                                                                    tGCGGTTGAAGTTCAGCCGGCGGCGCTGCGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGTTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATc                                                                        >  1:207725/1‑149 (MQ=255)
                                                                           gAAGTTCCGCCGGCGGCGCGGGGCGGCGGTGAGCCGGGCATATATCGAGACCGGG‑TTAGCTGCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGa                                                                  <  3:116471/148‑1 (MQ=255)
                                                                             aGTTCAGCCGGCGGCGCTGGGCTGCGGTGCGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATa                                                                >  4:25130/1‑148 (MQ=255)
                                                                               ttCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATc                                                             >  3:290516/1‑149 (MQ=255)
                                                                               ttCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATc                                                             >  4:107893/1‑149 (MQ=255)
                                                                                       gcggcgCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgtgtgt                                                     <  8:1785/149‑1 (MQ=255)
                                                                                          gcgcTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAAt                                                  <  4:290516/149‑1 (MQ=255)
                                                                                              tGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCaa                                              >  4:301002/1‑149 (MQ=255)
                                                                                                       gTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTACAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCCGCt                                     >  7:307832/1‑149 (MQ=255)
                                                                                                        tGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCTAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCtt                                    <  5:148274/149‑1 (MQ=255)
                                                                                                          cgcggggCATAGTTCGAGTCGGGG‑TTAGTTTCTTCGGGGGCGCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTc                                  <  6:138305/145‑1 (MQ=255)
                                                                                                           gCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGTGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGAGGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTc                                  >  3:35459/1‑148 (MQ=255)
                                                                                                           cccGGGCATAGTTCGAGACCGGG‑TTAGCTTCTTCGCGCGTGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTc                                  <  4:35459/147‑1 (MQ=255)
                                                                                                                gCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTgcg                            >  5:304245/1‑149 (MQ=255)
                                                                                                                     gATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGc                        <  1:338570/148‑1 (MQ=255)
                                                                                                                          cgccgggg‑TTCGCTTCTTCGCGGCCGCCCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATgg                    <  4:358230/142‑1 (MQ=255)
                                                                                                                          agaCCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATgg                    >  3:358230/1‑147 (MQ=255)
                                                                                                                             ccGGGTTTAGCTTCTCCGCGGGCCCCCGCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGa                <  1:352833/149‑1 (MQ=255)
                                                                                                                                 g‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGa           >  8:325249/1‑149 (MQ=255)
                                                                                                                                     atcatcttaGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCCTGTTAAGGATGGTGTCCTTGGCCTCGGCCACAATCCATGGATTATCCAGAGTGGAATGAAACCCTCGGCTTTCCTTGCGCGGCATGGCTGACATGatca        >  6:171313/5‑149 (MQ=255)
                                                                                                                                      gCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGatcat       >  1:138861/1‑149 (MQ=255)
                                                                                                                                        ttcttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGatcat       >  6:254064/1‑147 (MQ=255)
                                                                                                                                           ttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGAAAGGATCATCAGCt  >  1:334627/1‑149 (MQ=255)
                                                                                                                                           ttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCTTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCt  >  8:216867/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
TTGGCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGG‑TTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCT  >  NC_002947/1625060‑1625346

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