Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,070,268 G→C 43.2% intergenic (+65/+79) aroP‑I → / ← PP_0928 aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,070,2680GC43.2% 23.4 / 60.7 44intergenic (+65/+79)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base G (11/14);  new base C (11/8);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

ATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATA  >  NC_002947/1070120‑1070415
                                                                                                                                                    |                                                                                                                                                   
aTCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTAGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTTAGATCGAGCGCCGCgcggg                                                                                                                                                    >  2:1167940/1‑150 (MQ=255)
     ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGgcggcg                                                                                                                                               >  2:574757/1‑150 (MQ=255)
        aaTGCTGCTGCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGTTCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                            <  2:1114077/150‑1 (MQ=255)
          tGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                          <  1:574757/150‑1 (MQ=255)
           gCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                         <  1:2017174/150‑1 (MQ=255)
            cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                                        <  1:1557731/150‑1 (MQ=255)
             tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                       >  1:1000395/1‑150 (MQ=255)
               aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTTAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                     <  1:592158/150‑1 (MQ=255)
               aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTAAAGATGAAGGAGGCTCGCCAGGCCAAGGTCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCCCGGGCGGCGCTCGATCTca                                                                                                                                     <  2:2200860/150‑1 (MQ=39)
                                tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTa                                                                                                                    <  2:1000395/150‑1 (MQ=255)
                                    ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                            >  2:1190797/1‑122 (MQ=255)
                                    ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                            <  1:1190797/122‑1 (MQ=255)
                                              cTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                     >  2:1230786/1‑119 (MQ=255)
                                              cTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                     <  1:1230786/119‑1 (MQ=255)
                                                                cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGCATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                                                                    <  1:933594/150‑1 (MQ=255)
                                                                        gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCTTTAGGCTTTCATCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                          >  2:588779/1‑88 (MQ=14)
                                                                        gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCTTTAGGCTTTCATCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                          <  1:588779/88‑1 (MQ=14)
                                                                                                          tGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                                                                  >  2:1193329/1‑78 (MQ=255)
                                                                                                          tGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                                                                  <  1:1193329/78‑1 (MQ=255)
                                                                                                                          tcgcCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                    <  1:1205476/123‑1 (MQ=255)
                                                                                                                          tcgcCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                    >  2:1205476/2‑124 (MQ=255)
                                                                                                                            gcCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtg                                                                 <  2:1538341/109‑1 (MQ=255)
                                                                                                                            gcCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtg                                                                 >  1:1538341/1‑109 (MQ=255)
                                                                                                                                 tgAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCaa                             <  2:1834422/140‑1 (MQ=255)
                                                                                                                                 tgAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCaa                             >  1:1834422/1‑140 (MQ=255)
                                                                                                                                  gagaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCCGAAAACCTAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGa                  >  2:158504/1‑150 (MQ=25)
                                                                                                                                   agaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCCGAAAACCTAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAg                 <  1:158504/150‑1 (MQ=25)
                                                                                                                                     aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa                                 >  1:2408032/1‑132 (MQ=255)
                                                                                                                                     aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa                                 <  2:2408032/132‑1 (MQ=255)
                                                                                                                                      tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCAc              >  1:117406/1‑150 (MQ=255)
                                                                                                                                      tCGAGCGCCGCGAGGGCGGCGCCCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCAc              >  2:1895957/1‑150 (MQ=255)
                                                                                                                                       caaGCGCCGCTCGGCGGGCGCTCGATCTCACAGGGGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCAcc             <  1:1167940/148‑1 (MQ=255)
                                                                                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTc                                                                                     <  1:2815074/77‑1 (MQ=25)
                                                                                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTc                                                                                     >  2:2815074/1‑77 (MQ=25)
                                                                                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATTCCCCTCAATGCTGCCTCCGCCAAGCACGCACCGTGGTGTACAGCAACACCGTAGCGAGCACCg            >  2:313332/1‑150 (MQ=255)
                                                                                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa                                   <  2:1524419/127‑1 (MQ=255)
                                                                                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa                                   >  1:1524419/1‑127 (MQ=255)
                                                                                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTa                      <  1:552000/140‑1 (MQ=255)
                                                                                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTa                      >  2:552000/1‑140 (MQ=255)
                                                                                                                                        gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCg            >  1:617694/1‑150 (MQ=255)
                                                                                                                                           cgccgcCCGCGCGGCGCTCGATCTCACAGGCGCAGAAAACCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCa         >  2:923760/1‑150 (MQ=38)
                                                                                                                                             ccgcGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgca                                                 <  2:1179809/108‑1 (MQ=255)
                                                                                                                                             ccgcGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgca                                                 >  1:1179809/1‑108 (MQ=255)
                                                                                                                                             ccgcGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGt       >  2:2934705/1‑150 (MQ=255)
                                                                                                                                              cgcgcgGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTa      >  2:1079709/1‑150 (MQ=255)
                                                                                                                                                  cgggcgtCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCATCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATa  <  1:1870610/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATA  >  NC_002947/1070120‑1070415

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