Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,487,373 A→G 10.5% D127G (GAC→GGC)  PP_2181 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,487,3730AG10.5% 98.8 / 3.3 38D127G (GAC→GGC) PP_2181Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (15/19);  new base G (2/2);  total (17/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 = 6.06e-01

GACTGGGTGACCCGCCCGCACGAGCGTCGCACGCTGCACCTGGCCGGTGGCATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCAGTTGCAACCCCACGCCCGTTTCCGTTACGCCAACCTGACCGACAAG  >  NC_002947/2487225‑2487518
                                                                                                                                                    |                                                                                                                                                 
gACTGGGTGACCCGCCCGCACGAGCGTCGCACGCTGCACCTGGCCGGTGGCATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGgac                                                                                                                                                  <  1:1033258/150‑1 (MQ=255)
     ggTGACCCGCCCGCACGAGCGTCGCACGCTGCACCTGGCCGGTGGCATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAg                                                                                                                                             >  1:427922/1‑150 (MQ=255)
               ccgcACGAGCGTCGCACGCTGCACCTGGCCGGTGGCATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCGTTTGCATCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAg                                                                                                                                   <  2:1276919/150‑1 (MQ=255)
                   aCGAGCGTCGCACGCTGCACCTGGCCGGTGGCATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGtt                                                                                                                               <  2:109872/150‑1 (MQ=255)
                       gCGTCGCACGCTGCACCTGGCCGGTGGCATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGGCGACAGCTCGGAGCAGGGTTGCtt                                                                                                                           >  2:747974/1‑150 (MQ=255)
                        cGTCGCACGCAGAACCTGGCCGGTGGCATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTc                                                                                                                          <  1:583586/150‑1 (MQ=255)
                             cACGCTGCACCTGGCCGGTGGCATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCt                                                                                                                     <  1:194919/150‑1 (MQ=255)
                                 cTGCACCTGGCCGGTGGCATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGGCGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGaa                                                                                                                 <  1:747974/150‑1 (MQ=255)
                                  tGCACCTGGCCGGTGGCATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAg                                                                                                                <  1:1062269/150‑1 (MQ=255)
                                       cTGGCCGGTGGCATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGaa                                                                                                           >  1:1144931/1‑150 (MQ=255)
                                                  cATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGtgtg                                                                                                <  2:1028125/150‑1 (MQ=255)
                                                             tgCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATg                                                                                     <  2:427922/150‑1 (MQ=255)
                                                               cTGGGATCGCCCTCAATAACGGGCGCGTTTTCCACGCTGGACAGCCATTTGCACCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTCTGGTTGGGCGATGCt                                                                                   <  2:303483/150‑1 (MQ=255)
                                                                        cccTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTc                                                                          >  1:443025/1‑150 (MQ=255)
                                                                         tcTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGaa                                                                                                                 >  1:44737/2‑110 (MQ=255)
                                                                         ccTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGaa                                                                                                                 <  2:44737/110‑1 (MQ=255)
                                                                              aTAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGc                                                                                         <  2:918858/129‑1 (MQ=255)
                                                                              aTAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGc                                                                                         >  1:918858/1‑129 (MQ=255)
                                                                              aTAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGc                                                                      >  1:395527/1‑148 (MQ=255)
                                                                              aTAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGc                                                                      <  2:395527/148‑1 (MQ=255)
                                                                                     gcgcGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGGCGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCa                                                             >  1:559689/1‑150 (MQ=255)
                                                                                       gcgTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAAt                                                           <  2:628290/150‑1 (MQ=255)
                                                                                         gTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGa                                                         >  2:766251/1‑150 (MQ=255)
                                                                                                  gCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCa                                                >  2:253204/1‑150 (MQ=255)
                                                                                                         aGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGgtgt                                                                                                 <  1:166628/94‑1 (MQ=255)
                                                                                                         aGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGgtgt                                                                                                 >  2:166628/1‑94 (MQ=255)
                                                                                                         aGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCAGTTGCaa                                         <  1:519349/150‑1 (MQ=255)
                                                                                                          gCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGGCGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCAGTTGCAAc                                        <  2:559689/150‑1 (MQ=255)
                                                                                                           ccATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGCGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCAGTTGCAAcc                                       >  2:363303/1‑150 (MQ=255)
                                                                                                            cATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCAGTTGCAAccc                                      >  2:764197/1‑150 (MQ=255)
                                                                                                                 gCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCAGTTGCAACCCCACGc                                 <  2:894222/150‑1 (MQ=255)
                                                                                                                      cGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCAGTTGCAACCCCACGCCCGtt                            <  2:443025/150‑1 (MQ=255)
                                                                                                                            ccAGCAGTAGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCAGTTGCAACCCCACGCCCGTTTCCGtt                      <  2:295517/150‑1 (MQ=255)
                                                                                                                                       gaggagTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACTCTTCACGCCAATGACAGTTTTCACTTGCAACCCCCCGCCCCTTTCCGTTCCGCCAACCTg           >  2:215767/1‑150 (MQ=255)
                                                                                                                                             tACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCAGTTGCAACCCCACGCCCGTTTCCGTTACGCCAACCTGACCGAc     >  1:879185/1‑150 (MQ=255)
                                                                                                                                               ccTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCAc                                                                 <  1:759354/88‑1 (MQ=255)
                                                                                                                                               ccTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCAc                                                                 >  2:759354/1‑88 (MQ=255)
                                                                                                                                                cTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCAGTTGCAACCCCACGCCCGTTTCCGTTACGCCAACCTGACCGACAAg  >  1:45047/1‑150 (MQ=255)
                                                                                                                                                cTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCAGTTGCAACCCCACGCCCGTTTCCGTTACGCCAACCTGACCGACAAg  <  2:1144931/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GACTGGGTGACCCGCCCGCACGAGCGTCGCACGCTGCACCTGGCCGGTGGCATCACCGATGTGCTGGCATCGCCCTCCATAACGGGCGCGTTTTCCATGCTCGACAGCCATTTGCAGCCGGGAGCCAGCAGTGGCGAGGAGTACCTGGACGACAGCTCGGAGCAGGGTTGCTTCGTGCTCGAAGGCGAACTGACGGTGTGGTTGGGCGATGCTGAGCCGGTCACGCTGCACGCCAATGACAGTTTTCAGTTGCAACCCCACGCCCGTTTCCGTTACGCCAACCTGACCGACAAG  >  NC_002947/2487225‑2487518

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: ‑