Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,997,194 Δ1 bp 8.7% intergenic (‑163/+37) PP_5252 ← / ← PP_5253 hypothetical protein/arylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,997,1940C. 8.7% 121.7 / 7.3 46intergenic (‑163/+37)PP_5252/PP_5253hypothetical protein/arylesterase
Reads supporting (aligned to +/- strand):  ref base C (17/24);  new base . (2/2);  total (20/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 = 1.23e-01

GGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAG  >  NC_002947/5997049‑5997338
                                                                                                                                                  |                                                                                                                                                
ggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGt                                                                                                                                               <  7:215779/149‑1 (MQ=255)
  ccGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAg                                                                                                                                             <  3:27173/149‑1 (MQ=255)
        gggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGT‑TGTAGGAGCgg                                                                                                                                       <  8:76859/149‑1 (MQ=255)
            gACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCAACAAAA‑ACGGCTGTAGGAGCGGGttt                                                                                                                                   >  4:58230/1‑149 (MQ=255)
               tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTAcc                                                                                                                                >  2:54152/1‑149 (MQ=255)
                gCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTAccc                                                                                                                               >  2:3175/1‑149 (MQ=255)
                 cGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTCAATCGATGTGCGGTGGGTTGAAGAAGGTTAGACCTGGACAAGCTTTAATCTGCCTCGGACGGCCTCTTCGCGGCTAAACCCCCTCCTACA‑ACGTCTGTAGTAGCCGGTTTaaacg                                                                                                                              >  5:295930/1‑145 (MQ=255)
                           cGGTTGATGCAAGGCCTGGGCTTGAATGGATGGGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTc                                                                                                                    <  3:111242/149‑1 (MQ=255)
                                  tGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAAt                                                                                                             <  1:3175/149‑1 (MQ=255)
                                          cTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTg                                                                                                     <  3:58230/149‑1 (MQ=255)
                                                   aaTGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGc                                                                                            >  1:193898/1‑149 (MQ=255)
                                                         tgtgCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGc                                                                                        >  4:81652/1‑147 (MQ=255)
                                                          ttgCGGTGGGTGTAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGTAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCaa                                                                                      <  8:197867/147‑1 (MQ=255)
                                                            gctgTGGTTTGTAGTAGGTTAGCCCTGGTCAATCTTTTTGGTGGCTGTGCCGGCGTGTTGGGGTGTAAACCCGCTCCTAGA‑ACGTCTGTTGGTGAGGGTTTACCCGCGAATGTTTCCGACAATTACCTCTGTATGTACGCCAGCAAGt                                                                                    <  6:295930/145‑1 (MQ=255)
                                                               gTGGGTTTAAGAAGTTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTc                                                                                <  1:274841/149‑1 (MQ=255)
                                                                tGGGTTGAAGAAGGTTAGCCCTGGGCAAGCTTTTTGCGGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCa                                                                               <  4:93649/149‑1 (MQ=255)
                                                                    ttGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCa                                                                           <  2:193898/149‑1 (MQ=255)
                                                                       aagaagGTTAGCCTTGGTCAAGCTTTTTGCTGCGTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCACACAATTACCGCTGAAGGAACGCCAGCaa                                                                                      <  5:38062/135‑1 (MQ=38)
                                                                          aagGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAgctgct                                                                     <  3:81652/149‑1 (MQ=255)
                                                                           agGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGctgctg                                                                    <  1:54152/149‑1 (MQ=255)
                                                                              atAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTAAGCTGCTgcg                                                                  >  1:109532/2‑148 (MQ=255)
                                                                                aGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTc                                                                                   <  6:243628/129‑1 (MQ=255)
                                                                                aGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTc                                                                                   >  5:243628/1‑129 (MQ=255)
                                                                                  cccTGTCCAAGGTTTTTGCGCCCTGGGCCGGCCTCTGCGCGGGAAACCCCGCCCCAACA‑CCGTCTTTCGTCGGGGTTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGa                                                             <  5:39336/149‑1 (MQ=255)
                                                                                       gTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCAc                                                        >  2:53208/1‑149 (MQ=255)
                                                                                       gTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCAc                                                        >  6:80910/1‑149 (MQ=255)
                                                                                             cTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCAc                                                        <  4:114574/143‑1 (MQ=255)
                                                                                             cTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCAc                                                        >  3:114574/1‑143 (MQ=255)
                                                                                              tttttGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCAt                                                               <  3:293849/135‑1 (MQ=255)
                                                                                              tttttGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACA‑ACGTGTGTAGGAGCGGGTTTACCCGCGAATGTTTACCACAATTACCGCGGAAGGAACGCCCGCAAGTCCTCATTAAGCTTCTGCGCAt                                                               >  4:293849/1‑135 (MQ=37)
                                                                                                       ccTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGTAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgcg                                          <  3:170120/147‑1 (MQ=255)
                                                                                                         tGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTTCCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCaa                                                   >  5:325320/1‑136 (MQ=255)
                                                                                                         tGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCaa                                                   <  6:325320/136‑1 (MQ=255)
                                                                                                               ggccTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGc                                <  2:318544/149‑1 (MQ=255)
                                                                                                                gactctTCGCGCGTATACCCGCTCCTACT‑ACGTCTGTAGGTGCGGGTTTACCCGCGAATGTTTCCCACCATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGGCACTGCTAAACCGTGCGGCGCCCCGGCa                               <  6:321759/147‑1 (MQ=255)
                                                                                                                          gggTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc                                           <  3:182389/127‑1 (MQ=255)
                                                                                                                          gggTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc                                           <  3:178128/127‑1 (MQ=255)
                                                                                                                          gggTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc                                           >  4:178128/1‑127 (MQ=255)
                                                                                                                          gggTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc                                           >  4:182389/1‑127 (MQ=255)
                                                                                                                                  ccGCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCGGCGCATGAGTCACTGCAAACCCGTGCGGCGCCACGGCCTACACCTTCAGCTCggag             >  5:182491/1‑147 (MQ=38)
                                                                                                                                    gCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc              >  7:205013/1‑146 (MQ=255)
                                                                                                                                    gCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc              <  8:205013/146‑1 (MQ=255)
                                                                                                                                           cA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATc    <  2:109532/149‑1 (MQ=255)
                                                                                                                                            a‑aCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAt     >  2:8146/1‑147 (MQ=255)
                                                                                                                                              aCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAg  >  2:304913/1‑149 (MQ=255)
                                                                                                                                              aCGTCTGTAGGAGCGGGTTTACCCGCGAATGTGTCCAACAATTACCGCTGAAGGAACGCCAGCAAGTCCACATTCAGCTGCTGCTCATGAGTCACTGCAAAACCGTGCTGCGCCCCGGCATACACATTCACCTCGGCGCCACGAATCac  >  3:39979/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAG  >  NC_002947/5997049‑5997338

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