Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,820,598 G→C 9.1% intergenic (+399/+113) PP_2473 → / ← PP_2474 hypothetical protein/glutathione S‑transferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,820,5980GC 9.1% 133.0 / 4.1 44intergenic (+399/+113)PP_2473/PP_2474hypothetical protein/glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base G (20/20);  new base C (2/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

ATTGAAGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTGCTCGCTTGTAAC  >  NC_002947/2820452‑2820745
                                                                                                                                                  |                                                                                                                                                   
aTTGAAGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGCCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGaa                                                                                                                                                   >  3:158798/1‑149 (MQ=255)
       ggATGCCCTGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATTATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGt                                                                                                                                            <  2:352789/149‑1 (MQ=255)
          tGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGTCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGtttt                                                                                                                                         <  4:158798/149‑1 (MQ=255)
            cccAGGGGTTAGTTCGTAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTTCGACCCGGGtttttg                                                                                                                                       <  7:190033/149‑3 (MQ=255)
                  ggTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGATCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAg                                                                                                                                 <  4:36977/149‑1 (MQ=255)
                           gAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGt                                                                                                                                >  1:271033/1‑141 (MQ=255)
                           gAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGt                                                                                                                                <  2:271033/141‑1 (MQ=255)
                              aGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGa                                                                                                                     <  7:238468/149‑1 (MQ=255)
                              aGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGa                                                                                                                     <  5:114649/149‑1 (MQ=255)
                               gTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAg                                                                                                                    <  3:235504/149‑1 (MQ=255)
                                  tctcCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTcc                                                                                                                 >  1:44555/1‑149 (MQ=255)
                                    tcctcGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGtt                                                                                                                       <  4:131665/141‑1 (MQ=255)
                                    tcctcGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGtt                                                                                                                       >  3:131665/1‑141 (MQ=255)
                                         gTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTggg                                                                                                          <  2:44555/149‑1 (MQ=255)
                                         gTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTggg                                                                                                          <  3:5018/149‑1 (MQ=255)
                                                    atatCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTc                                                                                                                           <  4:282474/121‑1 (MQ=255)
                                                    atatCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTc                                                                                                                           >  3:282474/1‑121 (MQ=255)
                                                       tCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTTCGACCCGGGtttttgctt                                                                                                                                    <  1:181888/109‑6 (MQ=38)
                                                       tCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTTCGACCCGGGtttttgctt                                                                                                                                    >  2:181888/1‑104 (MQ=38)
                                                         ccATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGt                                                                                                                                <  8:195421/111‑1 (MQ=255)
                                                         ccATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGt                                                                                                                                >  7:195421/1‑111 (MQ=255)
                                                         ccATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGATCaa                                                                                           <  3:237998/148‑1 (MQ=255)
                                                          cATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAAtg                                                                                         >  7:314969/1‑149 (MQ=255)
                                                                         gTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCgg                                                                          >  8:2738/1‑149 (MQ=255)
                                                                         gTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCgg                                                                          >  1:187261/1‑149 (MQ=255)
                                                                                aGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCtttt                                                                                                             >  2:248455/1‑107 (MQ=255)
                                                                                aGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCtttt                                                                                                             <  1:248455/107‑1 (MQ=255)
                                                                                 gAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAAcc                                                                                                                                                 >  6:255805/1‑70 (MQ=255)
                                                                                 gAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAAcc                                                                                                                                                 <  5:255805/70‑1 (MQ=255)
                                                                                      tGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCtt                                                              <  2:187261/148‑1 (MQ=255)
                                                                                           tATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGc                                                        >  1:208783/1‑149 (MQ=255)
                                                                                                    gctgATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCtt                                                              >  5:81564/1‑134 (MQ=255)
                                                                                                    gctgATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCtt                                                              <  6:81564/134‑1 (MQ=255)
                                                                                                      tgATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGa                                                           >  5:295893/1‑135 (MQ=255)
                                                                                                      tgATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGa                                                           <  6:295893/135‑1 (MQ=255)
                                                                                                          ataGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCttttttt                                                 >  2:181043/1‑141 (MQ=255)
                                                                                                          ataGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCttttttt                                                 <  1:181043/141‑1 (MQ=255)
                                                                                                                aaaGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGtt                                    >  5:319676/1‑148 (MQ=255)
                                                                                                                             aaaaaaGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCC‑TTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGATTTTTTTGCGAAGCCAGGTTTAATTTTCATCtata                     >  6:199474/1‑149 (MQ=255)
                                                                                                                              aaaaaGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCtata                     >  8:344174/1‑149 (MQ=255)
                                                                                                                               aaaaGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTatat                    <  6:319676/149‑1 (MQ=255)
                                                                                                                                      aaaaCCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATGTTTATCCAGTGc             >  7:190825/1‑149 (MQ=255)
                                                                                                                                          cccGGGTTCGACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTGCTCGc         >  7:193973/1‑149 (MQ=255)
                                                                                                                                           ccGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTGCTCGCt        >  5:123305/1‑149 (MQ=255)
                                                                                                                                                  gAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTGCTCGCTTGTAAc  <  8:314969/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATTGAAGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGGTCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAGGCTTTTTTTGCGAAGCCAGGTTTAATTGGCATCTATATCCAGTGCTCGCTTGTAAC  >  NC_002947/2820452‑2820745

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