Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,018,201 G→A 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,2010GA97.2% 98.4 / ‑5.1 36intergenic (+144/‑36)mvaB/PP_mr44hydroxymethylglutaryl‑CoA lyase/CrcZ ncRNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/19);  minor base C (1/0);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

AGGTCACACGGGTGTCACATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGATTTTTTTGGAGTCGA  >  NC_002947/4018056‑4018345
                                                                                                                                                 |                                                                                                                                                
aGGTCACACGGGTGTCACATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCAATAATTTTAATGGAttgt                                                                                                                                              <  2:150249/150‑3 (MQ=255)
     acacGGGTGTCACATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCCCTTCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCaaaa                                                                                                                                         >  1:678082/1‑150 (MQ=255)
         gggTGTCACATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAGCCGTGACTCATTAATTTTAAAGGATTTTCAAAAGTTg                                                                                                                                     <  1:304633/150‑1 (MQ=255)
           gtgtCACATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgc                                                                                                                                   <  1:343609/150‑1 (MQ=255)
             gtCACATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcac                                                                                                                                 <  2:350973/150‑1 (MQ=255)
              tCACATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacg                                                                                                                                >  1:432320/1‑150 (MQ=255)
                     ggTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGACTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTg                                                                                                                         <  2:973163/150‑1 (MQ=255)
                      gtgtGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGc                                                                                                                        >  2:146692/1‑150 (MQ=255)
                          gggTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCtata                                                                                                                    >  1:116021/1‑150 (MQ=255)
                                      aCCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTacaa                                                                                                        >  2:1230193/1‑150 (MQ=255)
                                         tCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTacaacaa                                                                                                     >  1:1225270/1‑150 (MQ=255)
                                             gCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacgca                                                                                                                              <  2:356608/121‑1 (MQ=255)
                                             gCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacgca                                                                                                                              >  1:356608/1‑121 (MQ=255)
                                              cacaTTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACaaaa                                                                                                <  2:432320/150‑1 (MQ=255)
                                                        aGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTCAAGGCTTTTAGGAGATCGGCAGACACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGAt                                                                                      >  1:1260108/1‑150 (MQ=255)
                                                                          aGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCaaaa                                                                                                                                         <  1:518498/81‑1 (MQ=255)
                                                                          aGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCaaaa                                                                                                                                         >  2:518498/1‑81 (MQ=255)
                                                                          aGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAaacaac                                                                      <  2:996539/148‑1 (MQ=255)
                                                                          aGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAaacaac                                                                      >  1:996539/1‑148 (MQ=255)
                                                                           gggAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAgg                                                                   >  2:261460/1‑150 (MQ=255)
                                                                                     tCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGctct                                                         <  2:1225270/150‑1 (MQ=255)
                                                                                        gCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGctc                                                          >  1:1170333/1‑146 (MQ=255)
                                                                                        gCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGctc                                                          <  2:1170333/146‑1 (MQ=255)
                                                                                                ttaGGCTGGGCGGTGGCGACCAAAAAACTTAAATCATCAATATTAAAGGATTTACAAAAGTTGGCACTCATCCTACTACATCTAATATAAAACAAAAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAg                                              <  2:271594/148‑1 (MQ=255)
                                                                                                  aGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGaa                                            <  2:711230/150‑1 (MQ=255)
                                                                                                                 aaCCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGc                             <  1:146692/149‑1 (MQ=255)
                                                                                                                         aTCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTagcaag                     >  2:1269739/1‑150 (MQ=255)
                                                                                                                             aaaCTCATTAATTTTAAAAGATTTTCAAAAGTTGGCACGCACCCTGCTATAACTCTGGAAAAAAAAAAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGa                 <  2:1260108/149‑1 (MQ=255)
                                                                                                                                   aTTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCa                                  >  1:833589/1‑127 (MQ=255)
                                                                                                                                   aTTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCa                                  <  2:833589/127‑1 (MQ=255)
                                                                                                                                    ttAATTTTAAAGAATTTTCAAAAGTTGGCAAGCACACTGCTATATCTATAATAAAAAAAAAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGAttttttt          <  2:116021/150‑1 (MQ=255)
                                                                                                                                     tAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGATTTTTTTg         <  1:1269739/150‑1 (MQ=255)
                                                                                                                                     tAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTATAGAAAAACAAAAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGATTTTTTTg         <  2:899729/150‑1 (MQ=255)
                                                                                                                                          ttAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGATTTTTTTGGAGTc    >  1:119684/1‑150 (MQ=255)
                                                                                                                                            aaaGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAACAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGATTTTTTTGGAGTCGa  >  2:1148493/1‑150 (MQ=255)
                                                                                                                                            aaaGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAAGAAACGGCAATGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGATTTTTTTGGAGTCGa  >  1:1124579/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGGTCACACGGGTGTCACATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGATTTTTTTGGAGTCGA  >  NC_002947/4018056‑4018345

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