Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,403,590 G→A 21.8% intergenic (‑56/+25) PP_1227 ← / ← PP_1228 membrane protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,403,5900GA21.8% 59.4 / 9.6 32intergenic (‑56/+25)PP_1227/PP_1228membrane protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base G (14/11);  new base A (3/4);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCC  >  NC_002947/1403450‑1403734
                                                                                                                                            |                                                                                                                                                 
tGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCcgcg                                                                                                                                           <  2:108544/149‑1 (MQ=255)
   aCAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGaa                                                                                                                                         <  8:338585/148‑1 (MQ=255)
            cggcgATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGa                                                                                                                               <  2:321785/149‑1 (MQ=255)
            cggcgATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGGGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAAGAGGCCCTGa                                                                                                                               <  3:8559/149‑1 (MQ=255)
                              aaaGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACt                                                                                                             >  2:272749/1‑149 (MQ=255)
                                                    gTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGTAACTGCTCCATCAGCGCctgctgctg                                                                                      >  1:61718/1‑149 (MQ=255)
                                                       gaagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAaga                                                                                                                                       >  6:136751/1‑98 (MQ=255)
                                                       gaagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCACGCCCGCGAaga                                                                                                                                       <  5:136751/98‑1 (MQ=255)
                                                        aagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGTAACTGCTCCATCAGCGCCTGCTGCTGGTTg                                                                                  <  2:61718/149‑1 (MQ=255)
                                                              gTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGc                                                                            >  6:350360/1‑149 (MQ=255)
                                                                     gggTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCa                                                                     <  5:350360/149‑1 (MQ=255)
                                                                       gTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGt                                                                                                                                                  >  7:180499/1‑71 (MQ=255)
                                                                       gTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGt                                                                                                                                                  <  8:180499/71‑1 (MQ=255)
                                                                       gTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAgg                                                                   >  4:248385/1‑149 (MQ=255)
                                                                                ggCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCt                                                           >  7:181042/1‑148 (MQ=255)
                                                                                        ttGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACTCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGCCCAGCCGGTTCAGGGCTTgg                                                             >  7:233096/1‑138 (MQ=255)
                                                                                        ttGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGCCCAGCCGGTTCAGGGCTTgg                                                             <  8:233096/138‑1 (MQ=255)
                                                                                              tCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGc                                                      >  5:171993/1‑139 (MQ=255)
                                                                                              tCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGc                                                      <  6:171993/139‑1 (MQ=255)
                                                                                               cAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCg                                                                          >  4:133334/1‑118 (MQ=255)
                                                                                               cAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCg                                                                          <  3:133334/118‑1 (MQ=255)
                                                                                               cAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTc                                            >  4:249075/1‑148 (MQ=255)
                                                                                                aaGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCAGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCAt                                          >  3:52404/1‑149 (MQ=255)
                                                                                                   aaTAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTCGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCg                                        <  1:272749/149‑1 (MQ=255)
                                                                                                                  gtgcCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGcc                                                 <  4:199776/124‑1 (MQ=255)
                                                                                                                  gtgcCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGcc                                                 >  3:199776/1‑124 (MQ=255)
                                                                                                                           aGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCa               >  2:322705/1‑149 (MQ=255)
                                                                                                                                  ctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTc                                                                                                         <  4:172349/53‑1 (MQ=17)
                                                                                                                                  ctTCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTc                                                                                                         >  3:172349/1‑53 (MQ=17)
                                                                                                                                   ttCGCGGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg        >  1:87653/1‑148 (MQ=255)
                                                                                                                                      gcgGGCGTGCCCGCGAAGAGGCCCTGACAGCCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATg    >  1:30497/1‑149 (MQ=255)
                                                                                                                                        gggCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGcc  <  3:248385/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
TGAACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGCCTGTCAGGGCCTCTTCGCGGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCC  >  NC_002947/1403450‑1403734

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: