Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,844,682 A→T 10.2% intergenic (‑42/+38) PP_2495 ← / ← PP_2496 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,844,6820AT10.2% 129.1 / 6.3 49intergenic (‑42/+38)PP_2495/PP_2496hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/23);  new base T (2/3);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

GCGGATCTGCAGCCAGTCCTCGCGGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAACCAGGCGGAGTCGAGGCGGTATTCGGTGCTGCGTGCCCGC  >  NC_002947/2844539‑2844822
                                                                                                                                               |                                                                                                                                            
gcGGATCTGCAGCCAGTCCTCGCGGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCgg                                                                                                                                         >  7:118224/1‑149 (MQ=255)
         cAGCCAGTCCTCGCGGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGc                                                                                                                                  <  5:30635/147‑1 (MQ=255)
         cAGCCAGTCCTCGCGGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCCGCTTCGTTGCAGGTCggg                                                                                                                                   <  5:59672/146‑1 (MQ=255)
          aGCCAGTCCTCGCGGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCCGCTTCGTTGCAGGTCggg                                                                                                                                   >  4:155868/1‑145 (MQ=255)
          aGCCAGTCCTCGCGGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCCGCTTCGTTGCAGGTCggg                                                                                                                                   <  3:155868/145‑1 (MQ=255)
           gCCAGTCCTCGCGGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAGTAGAAACGCCCGACCTGCAACGAAGCTGGTCGGGCGTTa                                                                                                                              >  5:129675/1‑149 (MQ=255)
             cAGTCCTCGCGGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCCGCTTCGTTGCAGGTCGGGCGt                                                                                                                                <  4:86366/145‑1 (MQ=255)
             cAGTCCTCGCGGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCCGCTTCGTTGCAGGTCGGGCGt                                                                                                                                >  3:86366/1‑145 (MQ=255)
                tCCTCGCGGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATgggg                                                                                                                         >  1:287384/1‑149 (MQ=255)
                    cgcgGGGTAGCTCCCTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTgg                                                                                                                      >  5:268503/1‑148 (MQ=255)
                     gcgGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTTACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTgg                                                                                                                      <  1:88266/147‑1 (MQ=255)
                     gcgGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCTGAAATAGAAACGCCCGACCTCCAACGAAGCGGGACGGGCGTAATGGGGTGGGa                                                                                                                    >  6:208584/1‑149 (MQ=255)
                        gggTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTg                                                                                                                 <  4:137413/149‑1 (MQ=255)
                        gggTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCAAAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTg                                                                                                                 <  6:288385/149‑1 (MQ=255)
                         ggTAGCTCCCTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTgg                                                                                                                <  6:268503/149‑1 (MQ=255)
                              cTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTaaaa                                                                                                           <  8:118224/149‑1 (MQ=255)
                                 cTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAg                                                                                                          <  5:157321/147‑1 (MQ=255)
                                    ggTGCTAACCAGGGGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCa                                                                                                      <  4:290722/148‑1 (MQ=255)
                                      tGCTATCCTGGGGCAGAGCGGGACGGATTTTGCTGGTGCTGCTATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGc                                                                                                   <  2:287384/149‑1 (MQ=255)
                                         tAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGAc                                                                                                >  6:233682/1‑149 (MQ=255)
                                          aaCCCGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGAcc                                                                                               <  4:98613/149‑1 (MQ=255)
                                           aCCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGAc                                                                                                <  7:121714/147‑1 (MQ=255)
                                            ccAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCaa                                                                                             >  3:124003/1‑149 (MQ=255)
                                                      agCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCa                                                                                   <  5:233682/149‑1 (MQ=255)
                                                                      cTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGtt                                                                   >  5:100494/1‑149 (MQ=255)
                                                                                 aTTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAGTAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTgg                                                        <  6:129675/149‑1 (MQ=255)
                                                                                 aTTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTgg                                                        >  5:86437/1‑149 (MQ=255)
                                                                                 aTTAAATGATGTGGCTTTCATAGTGTTCACCTTGCAGAATTAGAAACGGCCGAACTGCAACGAAGCGGGTCGGGCGTTATGGGGTGTGATTGGTAAAAGGTCAGGCTACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTgg                                                        <  7:128838/149‑1 (MQ=255)
                                                                                    cAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGc                                                      >  8:255318/1‑148 (MQ=255)
                                                                                    cAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGc                                                      <  7:255318/148‑1 (MQ=255)
                                                                                     aaTGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCAc                                                    >  1:69163/1‑149 (MQ=255)
                                                                                     aaTGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGAACTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGAGACCAATTGCACGGCAGTATCCAATGCGCGTTGCTTCATTTGGGCAc                                                    >  8:244344/1‑149 (MQ=255)
                                                                                         aaGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTg                                                >  6:47511/1‑149 (MQ=255)
                                                                                                tttCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGTGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAAc                                         >  4:177714/1‑149 (MQ=255)
                                                                                                ttgcATTGTGTTCACCTTGCAGAAATCGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAAc                                         <  5:208584/146‑1 (MQ=255)
                                                                                                  tCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAAc                                         <  3:177714/147‑1 (MQ=255)
                                                                                                    aTAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGc                                              <  7:158955/140‑1 (MQ=255)
                                                                                                    aTAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGc                                              >  2:67590/1‑140 (MQ=255)
                                                                                                    aTAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGc                                              >  8:158955/1‑140 (MQ=255)
                                                                                                    aTAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGc                                              <  1:67590/140‑1 (MQ=255)
                                                                                                     tAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAACCAGGc                                    <  4:124003/149‑1 (MQ=255)
                                                                                                      aGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAACCAGGCg                                   >  7:248077/1‑149 (MQ=255)
                                                                                                         tttCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAACCAGGCGGa                                 <  7:244344/147‑1 (MQ=255)
                                                                                                         gtTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAACCAGGCGGAg                                >  4:23064/1‑149 (MQ=255)
                                                                                                                         tAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAACCAGGCGGAGTCGAGGCGGTATTCgg                <  6:86437/149‑1 (MQ=255)
                                                                                                                                      ccTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAACCAGGCGGAGTCGAGGCGGTATTCGGTGCTGCGTGCccg   >  3:159767/1‑149 (MQ=255)
                                                                                                                                       cTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAACCAGGCGGAGTCGAGGCGGTATTCGGTGCTGCGTGCccgc  >  2:221897/1‑149 (MQ=255)
                                                                                                                                              gAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAACCAGGCGGAGTCGAGGCGGTATTCgg                <  6:109364/128‑1 (MQ=255)
                                                                                                                                              gAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAACCAGGCGGAGTCGAGGCGGTATTCgg                >  5:109364/1‑128 (MQ=255)
                                                                                                                                               |                                                                                                                                            
GCGGATCTGCAGCCAGTCCTCGCGGGGTAGCTCCTTGGTGCTAACCAGGCGCAGAGCGGGACGGATTTTGCTGGTGCTGCGATTCAATGAAGTGGCTTTCATAGTGTTCACCTTGCAGAAATAGAAACGCCCGACCTGCAACGAAGCGGGTCGGGCGTTATGGGGTGGGATTGGTAAAAGGTCAGGCGACCAATTGCAGGGCAGTATCCAAGGCGCGTTGCTTGATTTGGGCACCCTGGCCGAACCAGGCGGAGTCGAGGCGGTATTCGGTGCTGCGTGCCCGC  >  NC_002947/2844539‑2844822

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