Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,057,851 G→A 12.1% intergenic (+67/+61) PP_1832 → / ← PP_1833 oxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,057,8510GA12.1% 81.4 / 3.1 33intergenic (+67/+61)PP_1832/PP_1833oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (14/15);  new base A (2/2);  total (16/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.39e-01

GTATGGCAGCGAGCAGGGCATGCAGGCGCGGTTTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGT  >  NC_002947/2057703‑2057994
                                                                                                                                                    |                                                                                                                                               
gTATGGCAGCGAGCAGGGCATGCAGGCGCGGATTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAAACCGCGACGAAAATGACTTAGCGATTCCGAAGGgcggcg                                                                                                                                                 >  3:97183/1‑149 (MQ=255)
       aGCGAGCTGGGCATGCAGGCGCGGTTTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCgt                                                                                                                                          <  5:102146/149‑1 (MQ=255)
          gAGCAGGGCATGCAGGCGCGGTTTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCACAACCgtcgg                                                                                                                                       >  3:288707/1‑149 (MQ=255)
            gCAGGGCATGCAGGCGCGGTTTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCgtcggtc                                                                                                                                     <  2:42817/149‑1 (MQ=255)
               gggCATGCAGGCGCGGTTTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTa                                                                                                                                  <  1:145187/149‑1 (MQ=255)
               gggCATGCAGGCGCGGTTTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTa                                                                                                                                  <  7:307155/149‑1 (MQ=255)
               gggCATGCAGGCGCGGTTTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGGGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTa                                                                                                                                  <  8:260658/149‑1 (MQ=255)
                   aTGCAGGCGCGGTTTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAg                                                                                                                              >  3:119899/1‑149 (MQ=255)
                      cAGGCGCGGTTTACTGGCGATGGTAGTGCGCGGGCATTTGTTGGTATCGACGAGCTGTAATCGCATTGCCGGCGGCTTGGCGGGGGAAGACGCTCCTACATTGACTGCGCGATGGCTAAGGGCGGCGCAACAGTCGGTCGTAGGAGCg                                                                                                                            <  4:200572/148‑1 (MQ=255)
                      cAGGCGCGGTTTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCgg                                                                                                                           <  6:147034/149‑1 (MQ=255)
                           gcgGTTTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCa                                                                                                                      <  4:119899/149‑1 (MQ=255)
                              gTTTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCACAACCGTCGGTCGTAGGAGCGGGTTCAccc                                                                                                                   <  4:288707/149‑1 (MQ=255)
                                    gggCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAAACCCCTCCGACATTTACTGCGCGCTGCCGAAGGGCGGGGCAACCGGCGGGGGTAGGAGCGGGGTGACCcgcg                                                                                                                >  3:303361/1‑146 (MQ=255)
                                             gACTGCGCGGGCATTTCTTGGCATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGtt                                                                                                                        >  3:84181/1‑129 (MQ=255)
                                             gACTGCGCGGGCATTTCTTGGCATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGtt                                                                                                                        <  4:84181/129‑1 (MQ=255)
                                                            tCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCGCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACAccc                                                                                      <  7:268909/148‑1 (MQ=255)
                                                               tGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCCCTCCTTCATTGACTGCGCGTTGCCTAAGGGCGGCGCAACCGTCGGTCGGAGGAGCGGGTTCACCCCCGAAGAGGCCAGTTGAGGCAATACAccccggc                                                                                  >  5:345493/1‑146 (MQ=255)
                                                                      gacgaGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTa                                                                            <  4:303361/148‑1 (MQ=255)
                                                                       acgaGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCAGCGaa                                                                                                              <  1:117566/113‑1 (MQ=255)
                                                                       acgaGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCCGCGCAAACGTCGGTCGTAGGAGCGGGTTCACCAGAGaa                                                                                                              >  2:117566/1‑113 (MQ=39)
                                                                                                      cgGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTg                                              >  8:124990/1‑146 (MQ=255)
                                                                                                      cgGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCACAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCaaa                                                             >  2:220945/1‑131 (MQ=255)
                                                                                                      cgGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCACAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCaaa                                                             <  1:220945/131‑1 (MQ=255)
                                                                                                         gTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCAt                                        >  4:11628/1‑149 (MQ=255)
                                                                                                             aCCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATc                                       >  4:18169/1‑146 (MQ=255)
                                                                                                                 gCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGa                                >  8:49222/1‑149 (MQ=255)
                                                                                                                  cTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGAc                               <  3:18169/149‑1 (MQ=255)
                                                                                                                         aTTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAAt                        >  3:288336/1‑149 (MQ=255)
                                                                                                                             aCTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCa                     <  4:288336/148‑1 (MQ=255)
                                                                                                                                 cgggATGCATAAGGGCGGGGCAACCTTCGGTCGTAGGAGGGGGTTAACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTg                  <  3:296535/144‑1 (MQ=255)
                                                                                                                                      tGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCtt           >  8:51749/1‑149 (MQ=255)
                                                                                                                                         cTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGAAATTCACTTGCCTCTTCgg        >  1:279048/1‑149 (MQ=255)
                                                                                                                                          tAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGc       >  7:314818/1‑149 (MQ=255)
                                                                                                                                              ggcggcGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAgg   >  8:111628/1‑149 (MQ=255)
                                                                                                                                               gcggcgCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGt  >  6:311987/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTATGGCAGCGAGCAGGGCATGCAGGCGCGGTTTACGGGCGATGTGACTGCGCGGGCATTTCTTGGTATCGACGAGCTGTAATCGCATTGCCGGCCGCTTCGCGGGTGAACCCGCTCCTACATTGACTGCGCGATGCCTAAGGGCGGCGCAACCGTCGGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGTTCAGGCAATACACCCCTCAGCGGTAGGTCGGCAAGGCAAACCGCTGCTGGCTCTGCAGCATCGAGATGACCGGCAATTCACTGGCCTCTTCGGCCAGGT  >  NC_002947/2057703‑2057994

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: