Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,478,892 T→G 16.1% intergenic (+102/+343) PP_2172 → / ← PP_2173 hypothetical protein/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,478,8920TG16.1% 62.8 / 8.2 32intergenic (+102/+343)PP_2172/PP_2173hypothetical protein/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (16/10);  new base G (2/3);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

CTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGACGGCACGA  >  NC_002947/2478749‑2479034
                                                                                                                                               |                                                                                                                                              
cTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCa                                                                                                                                          >  1:260766/1‑150 (MQ=255)
        cGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGtgttt                                                                                                                                  <  2:143524/150‑1 (MQ=255)
           tGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGCGCTCGCAGGgtgtggggg                                                                                                                               <  2:466636/150‑6 (MQ=38)
                      cttcAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACca                                                                                                                    <  2:209928/150‑1 (MQ=255)
                       ttcAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACcat                                                                                                                   <  2:532202/150‑1 (MQ=255)
                       ttcAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCGGCTTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCCGGGTGCTTTGTTTTGGAGACcat                                                                                                                   >  1:480775/1‑150 (MQ=255)
                        tcAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACcatc                                                                                                                  <  2:260766/150‑1 (MQ=255)
                             aaaTGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCAc                                                                                                             >  2:482529/1‑150 (MQ=255)
                                     gatgatTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGCGCTCGCAGGgt                                                                                                                                      <  2:833205/117‑1 (MQ=38)
                                     gatgatTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAAAGCTCAGCGCTCGCAGGgt                                                                                                                                      >  1:833205/1‑117 (MQ=38)
                                      atgatTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTc                                                                                                    >  1:243068/1‑150 (MQ=255)
                                              aCCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTgttttgtttt                                                                                                                            >  1:11380/1‑118 (MQ=255)
                                              aCCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTgttttgtttt                                                                                                                            <  2:11380/118‑1 (MQ=255)
                                                   aTCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATTATCACGAGCGAGTCAGCACGTGGCggg                                                                                       >  2:833689/1‑150 (MQ=255)
                                                         gCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGgatcgat                                                                                 <  2:537836/150‑1 (MQ=255)
                                                                  gCTGTAGCGTGCACGAGAGCTCCGCGTCAGCCCGACCAGGGTAACGTGTTACTCCCCACACCATGCGAGCGCTGAGCCCTCTAGCTGACATTTTCTTGGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGCACCTATGTCGGTAAc                                                                        >  1:619558/1‑150 (MQ=255)
                                                                              aCGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGga                                                                                                                         >  1:232344/1‑89 (MQ=255)
                                                                              aCGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGga                                                                                                                         <  2:232344/89‑1 (MQ=255)
                                                                                     cTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATg                                                     >  2:846496/1‑150 (MQ=255)
                                                                                      tCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATgg                                                    >  2:658065/1‑150 (MQ=255)
                                                                                            tAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAgg                                              <  2:243068/150‑1 (MQ=255)
                                                                                                ccGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACg                                          >  1:425998/1‑150 (MQ=255)
                                                                                                        aTTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATg                                                     <  2:118312/131‑1 (MQ=255)
                                                                                                        aTTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATg                                                     >  1:118312/1‑131 (MQ=255)
                                                                                                              cGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCGGGCATGGCCAAGGTACGTTGGTCCAGGCgtt                            <  2:480775/150‑3 (MQ=255)
                                                                                                                       ccccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGCCCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCaa                   >  1:745517/1‑150 (MQ=255)
                                                                                                                        cccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGCGTCAGCACGTGCCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTGCAAGCGGTGGTTCGCaag                  >  1:318836/1‑150 (MQ=255)
                                                                                                                        cccACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCaag                  >  2:111937/1‑150 (MQ=255)
                                                                                                                        aaaacacCCTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCaag                  >  2:784630/4‑150 (MQ=38)
                                                                                                                               ccTGCGAAAGCTCAGCGCTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGc                     <  2:297931/140‑1 (MQ=255)
                                                                                                                                        gcTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGACGGCACGa  >  2:661655/1‑150 (MQ=255)
                                                                                                                                        gcTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGACGGCACGa  >  1:209807/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CTCAGCAACGCTGTCGGGTCTTCTTCAAAAAATGATCGATGATTTGACCGCATCCATGCAGTACCGGCTGTAGCGTGCACGAGAGCTCCGCGTAGGCCGTGCTAATTGAACGGGTTACCCCCCACACCCTGCGAGCGCTGAGCTTTCGCAGGGTGTTTTGTTTTGGAGACCATCATCACGAGCGAGTCAGCACGTGGCGGGATCGATCTCGGTAACCACAGTTACCTCCGGCATGGCCAAGGTACGTTGGTCCAAGCGGTGTTTCGCAAGAAGGTAGACGGCACGA  >  NC_002947/2478749‑2479034

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: