Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,398,903 C→G 15.1% intergenic (+123/‑119) aroQ‑III → / → PP_3004 type II 3‑dehydroquinate dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,398,9030CG15.1% 80.8 / 6.4 33intergenic (+123/‑119)aroQ‑III/PP_3004type II 3‑dehydroquinate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/9);  new base G (2/3);  total (21/12)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GTGCGGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTCAA  >  NC_002947/3398758‑3399048
                                                                                                                                                  |                                                                                                                                                 
gtgcgGCGATGAGCATGCTCTGAGCAGGGCGCAGCGGTCGGATCTCCCAGCGCGTCTACCCGCGAAGCAGGTCGCGCGGGGCATGCGGCGACCGCGCCTGAAGTCACCGTAGGTGTCGTGGATGTCGTGCTGCCTACCCATCCAGCatg                                                                                                                                                >  4:173468/1‑149 (MQ=255)
     gcgATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCATGctggatgg                                                                                                                                           <  3:63351/149‑1 (MQ=255)
           aGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGCTACTGGCTTTCCATCCATGctggatgggtag                                                                                                                                       <  1:150311/148‑5 (MQ=38)
                       gCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACtt                                                                                                                          <  5:286867/148‑1 (MQ=255)
                       gCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACtt                                                                                                                          >  6:286867/1‑148 (MQ=255)
                                     tCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAg                                                                                                            >  6:166900/1‑148 (MQ=255)
                                               gAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGAGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGAATGCATAGTACGGTCCGcc                                                                                                 >  2:81943/1‑149 (MQ=255)
                                                            cgcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGaaa                                                                                                                                        >  6:70066/1‑97 (MQ=255)
                                                            cgcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGaaa                                                                                                                                        <  5:70066/97‑1 (MQ=255)
                                                             gcgAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCAcc                                                                                   >  5:169095/1‑149 (MQ=255)
                                                                         gcgcGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACACGCCACCCGGTTCGTTGcc                                                                       <  3:173468/149‑1 (MQ=255)
                                                                         gcgcGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGAATGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGc                                                                        <  1:81943/148‑1 (MQ=255)
                                                                                   tGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCa                                                             <  6:169095/149‑1 (MQ=255)
                                                                                                          ccGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTg                                      >  3:180086/1‑149 (MQ=255)
                                                                                                            gTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCACGGATGGAAAGCCAGTAGCGACTTCCGTGAc                                                                                                                   <  4:268739/70‑1 (MQ=255)
                                                                                                                 tgtgGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCAACCCGGTTCGATGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGACGTc                               >  4:213138/1‑149 (MQ=255)
                                                                                                                        gATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGa                        >  5:103540/1‑149 (MQ=255)
                                                                                                                          tGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCGCTCATGACATGGGAGGTCGAATGGACg                      >  2:238158/1‑149 (MQ=255)
                                                                                                                           gTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTgg                                     >  2:58514/1‑133 (MQ=255)
                                                                                                                            tCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTc                    >  5:97569/1‑149 (MQ=255)
                                                                                                                                tgctgcCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGc                  >  6:98427/1‑148 (MQ=255)
                                                                                                                                     ccTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGt                                                                               >  6:70241/1‑82 (MQ=255)
                                                                                                                                     ccTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACaa             <  4:180086/148‑1 (MQ=255)
                                                                                                                                     ccTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACaa             <  6:97569/148‑1 (MQ=255)
                                                                                                                                     ccTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTAAAACCCAACGGACCGGGGCTCGCATGACCTGGGAGGTCGAATGGACGACGCTACAAc            >  5:293301/1‑149 (MQ=255)
                                                                                                                                       tttccatccatGCTGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTa                >  2:85406/4‑143 (MQ=255)
                                                                                                                                       tttccatccatGCTGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTa                <  1:85406/140‑1 (MQ=255)
                                                                                                                                       tttccatccatGCTGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGa          >  2:126837/4‑149 (MQ=255)
                                                                                                                                       tACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGAc                       <  2:270959/136‑1 (MQ=255)
                                                                                                                                       tACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGAc                       >  1:270959/1‑136 (MQ=255)
                                                                                                                                        aCCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGAt         >  4:299044/1‑149 (MQ=255)
                                                                                                                                             tCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGt     >  5:74660/1‑148 (MQ=255)
                                                                                                                                                aGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTcaa  >  1:294274/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTGCGGCGATGAGCATGGTCTGAGCAGGGCGCAGGGGTCGGATCTCCGAGCGCGTCTACCCGCGAAGCAGGTGGCGCGGGGCATGCGGCGACGGCGCCTGAAGTGACCGTAGGTGTGGTGGATGTCGTGCTGCCTACCCATCCAGCATGGATGGAAAGCCAGTAGCGACTTCCGTGACTGCATAGTACGGTCCGCCTTCAACAAGCCACCCGGTTCGTTGCCTCCAACCCCACGGACCGGGGCTCTCATGACCTGGGAGGTCGAATGGACGTCGCTACAACGATCAGTCAA  >  NC_002947/3398758‑3399048

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: