Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,018,201 G→A 100% intergenic (+144/‑36) mvaB → / → PP_mr44 hydroxymethylglutaryl‑CoA lyase/CrcZ ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,018,2010GA100.0% 52.3 / NA 18intergenic (+144/‑36)mvaB/PP_mr44hydroxymethylglutaryl‑CoA lyase/CrcZ ncRNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/6);  total (12/6)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGATTTTT  >  NC_002947/4018073‑4018335
                                                                                                                                |                                                                                                                                      
caTTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCaa                                                                                                                                 <  2:46087/136‑1 (MQ=255)
       gtgGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTgg                                                                                                                          <  2:261685/136‑1 (MQ=255)
          ggTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcac                                                                                                                       <  1:102925/136‑1 (MQ=255)
           gTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacg                                                                                                                      >  2:127864/1‑136 (MQ=255)
                  gtTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTg                                                                                                               >  2:263808/1‑136 (MQ=255)
                  gtTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTg                                                                                                               >  1:376186/1‑136 (MQ=255)
                        tCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTatat                                                                                                         >  1:380159/1‑136 (MQ=255)
                              acaTTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTTTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTgg                                                                                                   >  1:63536/1‑136 (MQ=255)
                              acaTTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTgg                                                                                                   >  1:108345/1‑136 (MQ=255)
                                             aaaTCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACaaaaaa                                                                                    >  2:499658/1‑136 (MQ=255)
                                                      aaCAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATa                                                                           >  1:406574/1‑136 (MQ=255)
                                                                 aaaTCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATaaaaa                                                                >  2:364745/1‑136 (MQ=255)
                                                                 aaaTCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATaaaaa                                                                <  1:127864/136‑1 (MQ=255)
                                                                        cAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAgg                                                         <  2:376186/136‑1 (MQ=255)
                                                                                        gtggtgTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATa                                         >  2:271865/1‑136 (MQ=255)
                                                                                          ggagtCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAAc                                       >  2:317977/4‑136 (MQ=255)
                                                                                                      aaaTCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACg                           <  2:48681/136‑1 (MQ=255)
                                                                                                               cTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCg                  >  1:29749/1‑136 (MQ=255)
                                                                                                                               ggTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGAttttt  >  1:209815/1‑136 (MQ=255)
                                                                                                                                |                                                                                                                                      
CATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGATTTTT  >  NC_002947/4018073‑4018335

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