Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,557,779 A→G 11.5% intergenic (‑99/+232) PP_4888 ← / ← purA methyl‑accepting chemotaxis transducer/adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,557,7790AG11.5% 91.4 / 3.5 37intergenic (‑99/+232)PP_4888/purAmethyl‑accepting chemotaxis transducer/adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base A (14/18);  new base G (2/2);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

AAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGT  >  NC_002947/5557630‑5557923
                                                                                                                                                     |                                                                                                                                                 
aaacaCAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCaa                                                                                                                                                   <  1:234833/150‑1 (MQ=255)
     caGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCt                                                                                                                                              <  2:362477/150‑1 (MQ=255)
           gCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAAt                                                                                                                                        >  2:442454/1‑150 (MQ=255)
              agagGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATggg                                                                                                                                     <  1:215630/150‑1 (MQ=255)
                  gcggcTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGatt                                                                                                                                 >  1:392440/1‑150 (MQ=255)
                     gcTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCgt                                                                                                                              <  2:76032/150‑1 (MQ=255)
                          gTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTGTTCTGGCGCCACTATTTTGATTGCATCTACAGATACTTCTTACAGCAAGGTTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCa                                                                                                                        <  2:26811/150‑1 (MQ=255)
                           tAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAg                                                                                                                       <  2:553117/150‑1 (MQ=255)
                            aaCAGTGAGAGGATGGCAGACACGTTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAgg                                                                                                                      >  2:594129/1‑150 (MQ=255)
                                  gagaGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCATGCCCAtt                                                                                                                <  2:392440/150‑1 (MQ=255)
                                                      gtAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTa                                                                                            <  1:42403/150‑1 (MQ=255)
                                                             tCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGc                                                                                     <  2:395546/150‑1 (MQ=255)
                                                                 gtTTTCTGGCGCCACTATTTTGATGGCATTTATAGAACGTTCTCACAGAAAGGGTTGTGCTGCAGGCAGCGTGGACTTGCCCCAAGGTTTAGCAATGGGGTTTCGCGTGCAGGCCCATTCGCGGGTGAACCCGTTCCTACAACGGCCttt                                                                                 <  2:418611/150‑3 (MQ=255)
                                                                         gcgcCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAAcc                                                                         >  1:246957/1‑150 (MQ=255)
                                                                           gcCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCt                                                                                             <  2:455023/128‑1 (MQ=255)
                                                                           gcCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCt                                                                                             >  1:455023/1‑128 (MQ=255)
                                                                              aCTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTcgc                                                                                                             >  2:441040/1‑109 (MQ=255)
                                                                              aCTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTcgc                                                                                                             <  1:441040/109‑1 (MQ=255)
                                                                               cTATTTTGATCGCAGGTAGAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGtt                                                                   >  2:305493/1‑150 (MQ=255)
                                                                                 tttttGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTtgt                                                                 <  2:553438/149‑1 (MQ=255)
                                                                                            aTCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATa                                                      >  2:290184/1‑150 (MQ=255)
                                                                                             tCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAg                                                     >  1:590120/1‑150 (MQ=255)
                                                                                               tACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGAGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGa                                                   >  1:417210/1‑150 (MQ=255)
                                                                                                  aGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCg                                                >  2:462873/1‑150 (MQ=255)
                                                                                                   gAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCgg                                               <  2:590120/150‑1 (MQ=255)
                                                                                                   gAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGAGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCgg                                               <  2:417210/150‑1 (MQ=255)
                                                                                                    aaaCTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGCCATGCGACATCGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCggg                                              <  1:442454/150‑1 (MQ=39)
                                                                                                        ttcttACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTa                                          >  1:118826/1‑150 (MQ=255)
                                                                                                           ttACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTAccc                                       >  1:420165/1‑150 (MQ=255)
                                                                                                             aCAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCcgc                                     <  1:462873/150‑1 (MQ=255)
                                                                                                               aGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTAACCCGCGa                                   >  1:548412/1‑150 (MQ=255)
                                                                                                                    gggTTGTGCCGCAGGCAGCGGGGAAATGCGCCAGGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACg                                                                                       <  1:400375/93‑1 (MQ=255)
                                                                                                                    gggTTGTGCCGCAGGCAGCGGGGAAATGCGCCAGGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACg                                                                                       >  2:400375/1‑93 (MQ=255)
                                                                                                                    gggTTGTGCCGCAGGCAGCGGGGAAATGCGCCAGGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACg                                                                                       >  2:400347/1‑93 (MQ=255)
                                                                                                                    gggTTGTGCCGCAGGCAGCGGGGAAATGCGCCAGGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACg                                                                                       <  1:400347/93‑1 (MQ=255)
                                                                                                                          tgCCGCAGGCAGCGGGGAAATGCGCCAACGTCACCCACAGGATCTCGTCCAGCAGTTCCATTCGCGGGTCACACCGGTCGTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAgcg                        <  1:594129/150‑1 (MQ=255)
                                                                                                                                              tGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCg    <  1:290184/150‑1 (MQ=255)
                                                                                                                                                cgcCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGt  <  1:1154/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGGTTGTGCCGCAGGCAGCGGGGAAATGCGCCAAGGTCTAGCAATGGGGTTTCGTGTGCAGGCCCATTCGCGGGTAAACCCGCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGT  >  NC_002947/5557630‑5557923

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: