Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,625,208 A→C 15.5% 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,2080AC15.5% 89.9 / 5.7 43intergenic (+31/+26)PP_1425/nadBhypothetical protein/L‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base A (22/14);  new base C (4/3);  total (26/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.98e-01

TTGGCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGC  >  NC_002947/1625060‑1625345
                                                                                                                                                    |                                                                                                                                         
atgccTGGGTGAAAGCGACTGTATACCGCAGGATGCGCAAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGa                                                                                                                                           <  4:243076/145‑1 (MQ=255)
       ggggAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATTGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCcg                                                                                                                                      <  5:331199/147‑1 (MQ=255)
          gAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACGGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAg                                                                                                                                 >  8:194917/1‑149 (MQ=255)
          gAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCTCGCCCGCGAAg                                                                                                                                 <  6:140810/149‑1 (MQ=255)
                cAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCacg                                                                                                                           >  6:233816/1‑149 (MQ=255)
                   tgtATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGGTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGGGAGGCGGGCATAGCTCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGACGGACacgacg                                                                                                                        >  4:234872/1‑149 (MQ=255)
                          cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGCGCCCGCGAAGGCCACGCCGCCGatat                                                                                                                 >  6:155908/1‑149 (MQ=255)
                          cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGcc                                                                                                                              >  1:197569/1‑136 (MQ=255)
                          cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGcc                                                                                                                              <  2:197569/136‑1 (MQ=255)
                          cTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGatat                                                                                                                 >  8:93216/1‑149 (MQ=255)
                           tCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGatat                                                                                                                 >  3:74169/1‑148 (MQ=255)
                            cAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCcgc                                                                                                                                     <  3:115430/127‑1 (MQ=255)
                            cAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCcgc                                                                                                                                     >  4:115430/1‑127 (MQ=255)
                                        ggggCTGACGAGCATCGGGGGTTGCGGGTGGGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGTGCCGGTCATAGATCGAGACCGGGTTTGCTTCTTCGGGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTaa                                                                                                    <  4:74169/148‑1 (MQ=255)
                                                 gAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGATTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTaa                                                                                                    >  1:56220/1‑139 (MQ=255)
                                                   gCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCt                                                                                        >  8:253248/1‑149 (MQ=255)
                                                        ggCGCTTCCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGATGCGGTGAGCAGGGCATAGATCGAGACCGGGTTAGCTTCTTAGCGGGCGATCCCGCGAAGGCCACGACGCCGATATCAGCTCGGGTTAAGGAGGttgtccttgtcc                                                                                   >  4:26131/1‑149 (MQ=255)
                                                          cgcTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTc                                                                                 >  7:36082/1‑149 (MQ=255)
                                                               gCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCa                                                                            <  5:233816/149‑1 (MQ=255)
                                                                cGGGTGCGGTTCATGTGCCGCGGGGGGCGGGGGGCGGGGGTGCGCCGGGAATAGAGCGAGACGGGGTTAGCTGCTTCGCGGGCGCCCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCa                                                                            <  7:194917/148‑1 (MQ=255)
                                                                 gggTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAg                                                                           >  5:247667/1‑148 (MQ=255)
                                                                    tGCGGTTGAAGTTCAGCGGGCGGCGGTGGGCTGCGGGGAGCGGGGCATAGATCGACACCGGGTTAGCTTCTTGGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATc                                                                       <  3:224987/149‑1 (MQ=255)
                                                                     gCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATcc                                                                      >  2:55789/1‑149 (MQ=255)
                                                                       ggTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGAt                                                                                                >  5:229495/1‑121 (MQ=255)
                                                                       ggTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGAt                                                                                                <  6:229495/121‑1 (MQ=255)
                                                                           gAAGTTCAGCCGGCGGAGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTAAGCTACTTCCCGGGCGAGCCCGCGAAGGACACGACGCCGATATCAGATCGGGTAAAGGATGGTGTCCTTGGCCTCGTCCAGCATCACTGGa                                                                 <  2:246412/148‑1 (MQ=255)
                                                                               ttCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATc                                                            >  1:7169/1‑149 (MQ=255)
                                                                                  aGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCa                                                          <  7:93216/148‑1 (MQ=255)
                                                                                        cggcgCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTa                                                   <  7:253248/149‑1 (MQ=255)
                                                                                                       gTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCt                                    >  7:16437/1‑149 (MQ=255)
                                                                                                        tGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGGGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCtt                                   <  7:298527/149‑1 (MQ=255)
                                                                                                            cggggattAGATCGAGACCGGGTTTGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCt                               <  2:7169/142‑1 (MQ=255)
                                                                                                               ctatctgATCGATACCGGGATAGCTTATTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTg                             <  2:216985/142‑1 (MQ=255)
                                                                                                                      aTCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTg                             >  1:216985/1‑141 (MQ=255)
                                                                                                                         gagaCCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGc                  >  2:231206/1‑149 (MQ=255)
                                                                                                                           gaCCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTg                >  3:251132/1‑149 (MQ=255)
                                                                                                                             ccGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGa               >  1:13940/1‑148 (MQ=255)
                                                                                                                                       ttcttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCATTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCATCGGCTTTCCTTTCGCTGCATGGCTGACAGGAtcatca    >  6:189749/1‑149 (MQ=255)
                                                                                                                                         cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAg   >  8:210975/1‑148 (MQ=255)
                                                                                                                                         cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGc  >  4:310393/1‑149 (MQ=255)
                                                                                                                                         cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAAACCTCCGCTTTCCTTGCGCTGCATGGCTGaaa             >  4:454/1‑136 (MQ=255)
                                                                                                                                         attggcgGGCTCGCCGGGGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACa             <  3:454/134‑1 (MQ=255)
                                                                                                                                                 gCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTcc                                <  3:244096/111‑1 (MQ=255)
                                                                                                                                                 gCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTcc                                >  4:244096/1‑111 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
TTGGCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGCGCTTGCGGGTGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGC  >  NC_002947/1625060‑1625345

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