Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,625,207 G→T 12.2% 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,2070GT12.2% 104.3 / 5.4 46intergenic (+30/+27)PP_1425/nadBhypothetical protein/L‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base T (2/4);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

GCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATC‑GGCGCTTGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGC  >  NC_002947/1625063‑1625350
                                                                                                                                                  |                                                                                                                                               
gCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATCGGGCGCTTGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGc                                                                                                                                              >  1:93662/1‑149 (MQ=255)
      ggAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATC‑GGCGCTTGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCTCGCCCGCGaa                                                                                                                                       <  4:147126/149‑1 (MQ=255)
      ggAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATC‑GGCGCTTGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCcgcg                                                                                                                                         >  7:168949/1‑147 (MQ=255)
        aaaGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATC‑GGCGCTTGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCTCGCCCGCGAAga                                                                                                                                     <  4:8148/149‑2 (MQ=255)
              aGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATC‑GGCGCTTGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCacga                                                                                                                               >  3:131562/1‑149 (MQ=255)
                           ggATGCGCTAGGGGCTGACGAGCATC‑GGCGCTTGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAg                                                                                                                                      <  7:10462/129‑1 (MQ=255)
                           ggATGCGCTAGGGGCTGACGAGCATC‑GGCGCTTGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAg                                                                                                                                      >  8:10462/1‑129 (MQ=255)
                           ggATGCGCTAGGGGCTGACGAGCATC‑GGCGCTTGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCa                                                                                                                    >  3:126962/1‑147 (MQ=255)
                            gATGCGCTAGGGGCTGACGAGCATC‑GGCGCTTGCGGGTTTCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCaga                                                                                                                  <  1:114067/149‑1 (MQ=255)
                             aTGCGCTCGGGGCTGACGAGCATC‑GGCGGTTGCGGG‑TGCGGTTTAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGAATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGc                                                                                                                <  8:168949/149‑1 (MQ=255)
                                        gCTGACGAGCATC‑GGCGCTTGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGAt                                                                                                     >  8:238929/1‑149 (MQ=255)
                                                   tC‑GGCGCTTGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTgg                                                                                          <  1:259016/149‑1 (MQ=255)
                                                   tC‑GGCGCTTGCGGG‑TGCGGTTGAAGTTAAGCCGGCGGCGCTGGGCTGCGTTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGACCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTgg                                                                                          >  3:273815/1‑149 (MQ=255)
                                                      ggCGCTTGCGGG‑GGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCTAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGcc                                                                                        <  1:262530/149‑1 (MQ=255)
                                                            tGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCGGGGAATAGATCGAGACCGGTTTTGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTc                                                                                   <  2:1923/148‑1 (MQ=255)
                                                             gCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGAt                                                                                                     >  8:190024/1‑129 (MQ=255)
                                                              gggg‑TGGTGTTGATGTTCAGCCGGCGGGGCTGGGGTGGGGTGTGCGGGTATTATTTCGAGACGGGGTTTGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCa                                                                                 <  4:131562/147‑1 (MQ=255)
                                                               ggg‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTg                                                                                           >  5:173095/1‑137 (MQ=255)
                                                               ggg‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTg                                                                                           <  6:173095/137‑1 (MQ=255)
                                                                   tGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATc                                                                            >  5:312373/1‑149 (MQ=255)
                                                                      ggTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCt                                                                         >  8:40160/1‑149 (MQ=255)
                                                                       gTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCAGCGAAGGCCACGACGCCGATATCAGAGCGCGTTAAGGAAGGTGTCCTTGGCCTCGTCCAGCATCCCtc                                                                        >  4:112143/1‑148 (MQ=255)
                                                                           aaGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGAt                                                                     >  6:145304/1‑148 (MQ=255)
                                                                            aGTTCAGCCGGCGGCGCTGGGCTGTGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCCAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATaa                                                                   <  5:116502/149‑1 (MQ=255)
                                                                                   cgggcggcGCTGGGCGGCGGTTAGCGGGGAATTGATCTTTACGTGGTTAGCTTCTTCGCGGTCGAGCCCGCGAAGGCCACGACGCCTATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTTGATAATCCa                                                               <  4:273815/144‑1 (MQ=255)
                                                                                     gggggcgCTGGGTTGGGGTTTGCGGGGAATTGATCTAGACCGGGTTTGCTTTTTGGCGGGCTAGCCGGCGAAGGCCACGACGCCTATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAgtgtg                                                          <  2:151549/146‑1 (MQ=255)
                                                                                         gcgcTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAAt                                                      <  7:40160/149‑1 (MQ=255)
                                                                                          cgcTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATg                                                     <  7:238929/149‑1 (MQ=255)
                                                                                                     ggTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTgg                                                                                          >  6:151447/1‑101 (MQ=255)
                                                                                                     ggTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTgg                                                                                          <  5:151447/101‑1 (MQ=255)
                                                                                                         aGCCGGGAATAGATCGCCCCCGGGGTATCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTc                                      <  4:266027/149‑1 (MQ=255)
                                                                                                           ccGGGCATTGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCt                                    <  6:312373/149‑1 (MQ=255)
                                                                                                              ggCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTgc                                 <  2:93662/149‑1 (MQ=255)
                                                                                                                   aGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTg                             >  4:229399/1‑148 (MQ=255)
                                                                                                                   aGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTg                             >  8:36960/1‑148 (MQ=255)
                                                                                                                    gATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGc                            <  5:145304/148‑1 (MQ=255)
                                                                                                                     aTCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATc                                                                 >  6:82353/1‑110 (MQ=255)
                                                                                                                      tCGAGTCCGGGTTAGCTTCTTCGGGGGCCTGCCGGCGAAGGCCACGACGCCGATATCAGTGCGGGTTAAGGTTGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATg                         <  3:112143/149‑1 (MQ=255)
                                                                                                                         tgagggggTTAGGTTGTTCGCGGGGGAGCCCGGGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCt                      <  3:276702/144‑1 (MQ=255)
                                                                                                                                gTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGa               <  7:36960/149‑1 (MQ=255)
                                                                                                                                      ttcttcGGGGGAGAGCCCGGGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGAtcatca         <  6:295049/149‑1 (MQ=255)
                                                                                                                                       tcttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAg        >  8:37841/1‑149 (MQ=255)
                                                                                                                                        cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGc       >  3:9956/1‑149 (MQ=255)
                                                                                                                                        cttcGCGGGCTCGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGc       >  1:166717/1‑149 (MQ=255)
                                                                                                                                         ttcGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCt      >  7:312505/1‑149 (MQ=255)
                                                                                                                                              gggCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGc  >  3:202457/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCTGGGGGAAAGCCAGTGTATACCTCAGGATGCGCTAGGGGCTGACGAGCATC‑GGCGCTTGCGGG‑TGCGGTTGAAGTTCAGCCGGCGGCGCTGGGCTGCGGTGAGCCGGGCATAGATCGAGACCGGGTTAGCTTCTTCGCGGGCGAGCCCGCGAAGGCCACGACGCCGATATCAGAGCGGGTTAAGGATGGTGTCCTTGGCCTCGTCCAGCATCCCTGGATAATCCAGTGTGTAATGCAACCCTCGGCTTTCCTTGCGCTGCATGGCTGACAGGATCATCAGCTCGGC  >  NC_002947/1625063‑1625350

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