Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,868,919 G→C 29.5% intergenic (+415/‑516) 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,9190GC29.5% 28.1 / 11.6 20intergenic (+415/‑516)PP_2520/PP_2522hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base C (3/3);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGGTGTTCAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTC  >  NC_002947/2868771‑2869057
                                                                                                                                                    |                                                                                                                                          
aaaGGTGTTCAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCgg                                                                                                                                           <  1:27925/150‑1 (MQ=255)
                  tCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTACCCCGGTCATCAATGCCGGGGCATTGGTGATCTGCGa                                                                                                                         <  1:580255/150‑1 (MQ=255)
                      gccattggccattCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCACCAATGCCCCGGCATtgatgaac                                                                                                                              <  1:210974/141‑3 (MQ=255)
                      gccattggccattCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCACCAATGCCCCGGCATtgatgaac                                                                                                                              >  2:210974/1‑139 (MQ=255)
                                         aaGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCgg                                                                                                  >  2:385633/1‑150 (MQ=255)
                                         aaGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCgg                                                                                                  >  1:52281/1‑150 (MQ=255)
                                                                      ccTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTc                                                                     <  1:385633/150‑1 (MQ=255)
                                                                              gCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCt                                                             >  2:312193/1‑150 (MQ=255)
                                                                                   cGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAg                                                        <  1:563932/150‑1 (MQ=255)
                                                                                        aaCTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATGTCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAAc                                                   >  1:247395/1‑150 (MQ=255)
                                                                                        aaCTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAAc                                                   >  2:118020/1‑150 (MQ=255)
                                                                                           tCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATGTCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCt                                                <  2:247395/150‑1 (MQ=255)
                                                                                                   gCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATAGGCGACGTCAATCAGCAACGCTTCGCGGCCCCCCCGTTGTCGCTGCGCGTTTTTGTCCGTCGGCTTTCAGGCAACCCTCACAtctc                                        <  1:505885/150‑1 (MQ=255)
                                                                                                                                           atcaATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGt                                                                      <  1:550747/80‑1 (MQ=255)
                                                                                                                                           atcaATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGt                                                                      >  2:550747/1‑80 (MQ=255)
                                                                                                                                           accaATGCCCCGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGtt    >  1:546387/3‑146 (MQ=255)
                                                                                                                                           accaATGCCCCGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGtt    >  1:546384/3‑146 (MQ=255)
                                                                                                                                           accaATGCCCCGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGtt    <  2:546384/144‑1 (MQ=255)
                                                                                                                                           accaATGCCCCGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGtt    <  2:546387/144‑1 (MQ=255)
                                                                                                                                               aTGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTc  <  2:272157/144‑1 (MQ=255)
                                                                                                                                               aTGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTc  >  1:272157/1‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
AAAGGTGTTCAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTTC  >  NC_002947/2868771‑2869057

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