Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,478,893 T→C 18.4% intergenic (+103/+342) 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,8930TC18.4% 66.7 / 7.4 32intergenic (+103/+342)PP_2172/PP_2173hypothetical protein/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (16/10);  new base C (3/3);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-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: