Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 167,292 C→A 11.4% intergenic (‑386/+34) pntAA ← / ← gcdR pyridine nucleotide transhydrogenase subunit alpha/gcdH transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947167,2920CA11.4% 101.3 / 2.9 35intergenic (‑386/+34)pntAA/gcdRpyridine nucleotide transhydrogenase subunit alpha/gcdH transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (14/17);  new base A (2/2);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

AGATTAGGAGAAAAAATCAAAAAAAAAGCGGCCAATTCGCTTACGACTAGTTGGATGTCGCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGACTCCACCTTGCGCTCTGGAATCATCAGA  >  NC_002947/167143‑167440
                                                                                                                                                     |                                                                                                                                                    
aGATTAGGAGAAAAAAGCAAAAAAAAAGCGGCCAATTCGCTTACGACTAGTTGGATGTCGCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATatacatac                                                                                                                                                      <  2:1480903/150‑1 (MQ=255)
       gagaAAAAATCAAAAAAAAAGCGGCCAATTCGCTTACGACTAGTTGGATGTCGCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTa                                                                                                                                               <  1:1342304/150‑1 (MQ=255)
                tCAAAAAAAAAGCGGCCAATTCGCTTACGACTAGTTGGATGTCGCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTaa                                                                                                                                      >  1:1136631/1‑150 (MQ=255)
                  aaaaaaaaaGCGGCCAATTCGCTTACGACTAGTTGGATTTTTCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTaaaa                                                                                                                                    <  2:1620797/150‑1 (MQ=255)
                   aaaaaaaaGCGGACAATTCGCTTACGACTAGTTGGATGTCGCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAaat                                                                                                                                   >  2:745696/1‑150 (MQ=255)
                                cAATTCGCTTACGAATAGTTGGATGTCGCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAgg                                                                                                                      <  1:745696/150‑1 (MQ=255)
                                       cTTACGACTAGTTGGATGTCGCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAAtt                                                                                                               >  1:1396808/1‑150 (MQ=255)
                                          aCGACTAGTTGGATGTCGCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGc                                                                                                            <  1:1409164/150‑1 (MQ=255)
                                          aCGACTAGTTGGATGTCGCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGc                                                                                                            >  2:71133/1‑150 (MQ=255)
                                                     gATGTCGCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTcc                                                                                                 >  2:6004/1‑150 (MQ=255)
                                                         tCGCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGa                                                                                             <  1:1596517/150‑1 (MQ=255)
                                                                    cccGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGtt                                                                                                                                         >  1:1574361/1‑95 (MQ=255)
                                                                    cccGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGtt                                                                                                                                         <  2:1574361/95‑1 (MQ=255)
                                                                     ccGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAaata                                                                                                                                  <  1:133163/101‑1 (MQ=255)
                                                                     ccGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAaata                                                                                                                                  >  2:133163/1‑101 (MQ=255)
                                                                      cGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAg                                                                                <  2:903640/150‑1 (MQ=255)
                                                                                  tACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTgg                                                                    >  1:731749/1‑150 (MQ=255)
                                                                                   aCTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGc                                                                   <  2:1396808/150‑1 (MQ=255)
                                                                                       aGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTAc                                                               <  2:247424/150‑1 (MQ=255)
                                                                                       aGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATAAATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTAc                                                               >  1:1513487/1‑150 (MQ=255)
                                                                                          aCGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTgg                                                                    <  2:1246183/142‑1 (MQ=255)
                                                                                          aCGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTgg                                                                    >  1:1246183/1‑142 (MQ=255)
                                                                                          aCGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTAccag                                                            <  1:71133/150‑1 (MQ=255)
                                                                                                       gtTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATAAATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAg                                               <  2:1513487/150‑1 (MQ=255)
                                                                                                                 ttGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGct                                     >  2:700735/1‑149 (MQ=255)
                                                                                                                   gCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGAc                                   <  1:700735/150‑1 (MQ=255)
                                                                                                                      ttACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGcc                                >  2:316114/1‑150 (MQ=255)
                                                                                                                        aCAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGa                              <  1:6004/150‑1 (MQ=255)
                                                                                                                            tttGCGAATATGCGACATATACATAAATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGAc                             <  2:1088374/147‑1 (MQ=255)
                                                                                                                            tttGCGAATATGCGACATATACATAAATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGAc                             >  1:1088374/1‑147 (MQ=255)
                                                                                                                               gCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGACTCCAcc                       >  2:537749/1‑150 (MQ=255)
                                                                                                                                 gAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGACTCCa                         >  2:512746/1‑146 (MQ=255)
                                                                                                                                 gAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGACTCCa                         <  1:512746/146‑1 (MQ=255)
                                                                                                                                 gAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGACTCCACCtt                     <  2:731749/150‑1 (MQ=255)
                                                                                                                                   atatGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGACTCCACCTTgc                   <  1:316114/150‑1 (MQ=255)
                                                                                                                                    tatGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGACTCCACCTTgcg                  >  2:59592/1‑150 (MQ=255)
                                                                                                                                                 catacatCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGACTCCACCTTGCGCTCTGGAatcatc     >  1:1575484/1‑150 (MQ=255)
                                                                                                                                                  atacatCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGACTCCACCTTGCGCTCTGGAAtcatca    >  1:630643/1‑150 (MQ=255)
                                                                                                                                                  atacatCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGACTCCACCTTGCGCTCTGGAAtcatca    >  2:910327/1‑150 (MQ=255)
                                                                                                                                                    acatCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGACTCCACCTTGCGCTCTGGAATCATCAGa  >  2:411440/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGATTAGGAGAAAAAATCAAAAAAAAAGCGGCCAATTCGCTTACGACTAGTTGGATGTCGCAAAATAACCCGCAAAGCCGCGTACTGAGCACGGTTAAAGGGTGTTTTTGTGCTTGCGTTACAGTTTGCGAATATGCGACATATACATACATCTGTAGGTGTTTAAAATAATTGACTACAGGGTCAATTTGCCGCGTTTGTCCCTGAGCTATTTATGTAGGCCTGAGCCTGGCTTACCAGCCAGTCGCGAAAGGCGCGTAGCGACGCCGACTCCACCTTGCGCTCTGGAATCATCAGA  >  NC_002947/167143‑167440

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