Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,030,147 A→C 9.1% intergenic (+104/‑590) pcaY → / → mgtA aromatic acid chemoreceptor/ATP‑dependent magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,030,1470AC 9.1% 118.0 / 3.6 44intergenic (+104/‑590)pcaY/mgtAaromatic acid chemoreceptor/ATP‑dependent magnesium transporter
Reads supporting (aligned to +/- strand):  ref base A (27/13);  new base C (2/2);  total (29/15)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GTGGGGGAGATGACGGCGTTGGTGGGGCGGTTCAAGGTTTGAGGCGCTGCGCCTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCTGGTCGGGCAAGCCGCGCCAGAGCCTGATGC  >  NC_002947/3030002‑3030296
                                                                                                                                                 |                                                                                                                                                     
gtggGGGAGATGACGGCGTTGGTGGGGCGGTTCAAGGTTTGAGGCGCTGCGCCTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAgcgc                                                                                                                                                   >  2:94904/1‑150 (MQ=255)
 tggGGGAGATGACGGCGTTGGTGGGGCGGTTCAAGGTTTGAGGCGCTGCGCCTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGcgcg                                                                                                                                                  >  1:954507/1‑150 (MQ=255)
        gaTGACGGCGTTGGTGGGGCGGTTCAAGGTTTGAGGCGCTGCGCCTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGCGCTCGGGGcggt                                                                                                                                           <  1:987759/150‑4 (MQ=255)
            aaagcgcTGGTGGGGCGGTTCAAGGTTTGAGGCGCTGCGCCTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATACAGTTCGACCGCCCCGAGCGCGGGGCCtttttt                                                                                                                                       >  2:1064550/4‑150 (MQ=255)
                   tggtggGGCGGTTCAAGGTTTGAGGCGCTGCGCCTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTAATTACTATGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCa                                                                                                                                <  2:1058778/150‑1 (MQ=255)
                    ggtggGGCGGTTCAAGGTTTGAGGCGCTGCGCCTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCAc                                                                                                                               >  1:1211403/1‑150 (MQ=255)
                         ggCGGTTCAAGGTTTGAGGCGCTGCGCCTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTgaag                                                                                                                          >  1:224158/1‑150 (MQ=255)
                            ggTTCAAGGTTTGAGGCGCTGCGCCTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAAt                                                                                                                       <  1:335441/150‑1 (MQ=255)
                                                  gcCTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCtcgtc                                                                                                 >  2:1389029/1‑150 (MQ=255)
                                                    cTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCTACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCTTCGc                                                                                               >  1:98153/1‑150 (MQ=255)
                                                    cTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGc                                                                                               >  2:901615/1‑150 (MQ=255)
                                                          ggAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTg                                                                                         >  1:328837/1‑150 (MQ=255)
                                                           gAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGa                                                                                        >  2:76502/1‑150 (MQ=255)
                                                                       aaGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAAt                                                                            <  2:1381690/150‑1 (MQ=255)
                                                                          gTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGACGAATggg                                                                         >  1:48804/1‑149 (MQ=255)
                                                                           ttAAACTTATTTAAAGTTCCCCTTGCGGGTTAACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGcc                                                                        >  2:874123/1‑150 (MQ=255)
                                                                                 ttATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCa                                                                  >  1:1277674/1‑150 (MQ=255)
                                                                                   aTTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAgg                                                                <  2:1277674/150‑1 (MQ=255)
                                                                                    tttAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTc                                                                                                                                      <  1:108990/79‑1 (MQ=255)
                                                                                    tttAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTc                                                                                                                                      >  2:108990/1‑79 (MQ=255)
                                                                                    tttAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGc                                                               >  1:205052/1‑150 (MQ=255)
                                                                                            tCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAAc                                                       >  2:669312/1‑150 (MQ=255)
                                                                                               ccTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCgg                                                    <  1:669312/150‑1 (MQ=255)
                                                                                                cTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGa                                                   >  2:827205/1‑150 (MQ=255)
                                                                                                        ttGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGt                                           <  1:1389029/150‑1 (MQ=255)
                                                                                                         tGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTc                                          <  1:901615/150‑1 (MQ=255)
                                                                                                             ttCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTc                                          >  1:1189131/1‑146 (MQ=255)
                                                                                                             ttCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTc                                          <  2:1189131/146‑1 (MQ=255)
                                                                                                             ttCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGa                                      >  1:560960/1‑150 (MQ=255)
                                                                                                                  tactacGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACAccc                                 >  1:751162/1‑150 (MQ=255)
                                                                                                                   actacGCTTGAATACAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCt                                <  2:328837/150‑1 (MQ=255)
                                                                                                                        gCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCTGGTCg                           >  1:567295/1‑150 (MQ=255)
                                                                                                                             aaTAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCTGGTCGGGCaa                      >  2:1345325/1‑150 (MQ=255)
                                                                                                                             aaTAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCTGGTCGGGCaa                      >  2:203470/1‑150 (MQ=255)
                                                                                                                               tAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCg                                                     <  2:1336850/117‑1 (MQ=255)
                                                                                                                               tAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCg                                                     >  1:1336850/1‑117 (MQ=255)
                                                                                                                               tAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCTGGTCGGGCAAGc                    >  1:1019478/1‑150 (MQ=255)
                                                                                                                                    ttCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCTGGTCGGGCAAGCCgcgc               >  2:78742/1‑150 (MQ=255)
                                                                                                                                    ttCGACCGCCCCGAGCGCGGGGCCTTTTTTACCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCTGGTCGGGCAAGCCgcgc               >  1:1302457/1‑150 (MQ=255)
                                                                                                                                    aaaaaaggCCCCGCGCTCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCTGGTCGGGCAAGCCgcgc               >  1:483292/8‑150 (MQ=255)
                                                                                                                                        aaggCCCCGCGCTCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAGCCGGAACACCGGTCCCGACACCCTGGTCg                           <  2:229074/131‑1 (MQ=255)
                                                                                                                                        aaggCCCCGCGCTCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAGCCGGAACACCGGTCCCGACACCCTGGTCg                           >  1:229074/4‑134 (MQ=255)
                                                                                                                                              ccGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCTGGTCGGGCAAGCCGCGCCAGAGCCTGa     <  2:224158/150‑1 (MQ=255)
                                                                                                                                               cGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCTGGTCGGGCAAGCCGCGCCAGAGCCTGAt    <  2:560960/150‑1 (MQ=255)
                                                                                                                                                 aGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCTGGTCGGGCAAGCCGCGCCAGAGCCTGATGc  <  1:203470/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GTGGGGGAGATGACGGCGTTGGTGGGGCGGTTCAAGGTTTGAGGCGCTGCGCCTGCGTGGAACATTAAGATAAGGTTAAACTTATTTAAAGTTCCCCTTGCGGGTTGACTTCGATACTACGCTTGAATAAAGTTCGACCGCCCCGAGCGCGGGGCCTTTTTTCCCTTCACTGAAGAATTCAATGGACACAGTATCTCGTCGCTGCCTGACCTCTGCCGAATGGCCTGGTCAGGCACTTGAACCGGAACACCGGTCCCGACACCCTGGTCGGGCAAGCCGCGCCAGAGCCTGATGC  >  NC_002947/3030002‑3030296

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