Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,868,920 G→C 16.6% intergenic (+416/‑515) PP_2520 → / → PP_2522 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,868,9200GC16.6% 56.3 / 6.6 26intergenic (+416/‑515)PP_2520/PP_2522hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (8/12);  new base C (2/2);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

AGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTCGC  >  NC_002947/2868781‑2869059
                                                                                                                                           |                                                                                                                                           
aGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGt                                                                                                                                    >  2:338300/1‑149 (MQ=255)
aGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGt                                                                                                                                    <  1:338300/149‑1 (MQ=255)
   cTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATc                                                                                                                                <  2:45922/150‑1 (MQ=255)
              cattggccattCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATatcaa                                                                                                                     >  2:410490/1‑150 (MQ=255)
                    ccattCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTc                                                                                                               <  1:322777/150‑1 (MQ=255)
                                  ccATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGccc                                                                                                 >  1:428363/1‑150 (MQ=255)
                                                       aTGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGtt                                                                                          >  2:59820/1‑136 (MQ=255)
                                                       aTGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGtt                                                                                          <  1:59820/136‑1 (MQ=255)
                                                            ccTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTc                                                                       <  2:428363/150‑1 (MQ=255)
                                                                 ggggCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACg                                                                  <  2:383532/150‑1 (MQ=255)
                                                                  gggCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGAtt                                                                            <  1:34778/139‑1 (MQ=255)
                                                                  gggCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGAtt                                                                            >  2:34778/1‑139 (MQ=255)
                                                                     cAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGCCCGACGGCtt                                                              >  2:479886/1‑150 (MQ=255)
                                                                                 tCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGt                                                                        >  1:267015/1‑128 (MQ=255)
                                                                                 tCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGt                                                                        <  2:267015/128‑1 (MQ=255)
                                                                                     tGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTcaca                                              <  2:454063/150‑1 (MQ=255)
                                                                                             cTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGTGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAg                                                          >  1:202273/1‑130 (MQ=255)
                                                                                             cTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGTGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAg                                                          <  2:202273/130‑1 (MQ=255)
                                                                                                gAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCACCAATGCCCCGGCATTg                                                                                                                                      <  1:560728/51‑1 (MQ=25)
                                                                                                gAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCACCAATGCCCCGGCATTg                                                                                                                                      >  2:560728/1‑51 (MQ=25)
                                                                                                   ttCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATAcc                                <  1:410490/150‑1 (MQ=255)
                                                                                                            gggcggCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCt                       <  1:216991/150‑1 (MQ=255)
                                                                                                                 gcCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGCCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTc                  <  1:479886/150‑1 (MQ=255)
                                                                                                                              atcaccaATGCCCCGGCATTGATGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGttt     >  2:80213/6‑150 (MQ=255)
                                                                                                                               tcaccaATGCCCCGGCATTGATGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTc    <  1:80213/150‑1 (MQ=255)
                                                                                                                                 atcaATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTcgc  >  1:144101/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                           
AGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTCGC  >  NC_002947/2868781‑2869059

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: