Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,312,732 T→A 11.1% intergenic (+30/‑132) PP_2912 → / → hemB hypothetical protein/5‑aminolevulinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,312,7320TA11.1% 126.4 / 7.4 46intergenic (+30/‑132)PP_2912/hemBhypothetical protein/5‑aminolevulinate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (20/20);  new base A (3/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

TGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTCC  >  NC_002947/3312590‑3312879
                                                                                                                                              |                                                                                                                                                    
tGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCcgcg                                                                                                                                              <  2:837176/150‑1 (MQ=255)
      agaTCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGC‑GGGTAACCCCGCGAAGAgg                                                                                                                                        <  2:1296151/150‑1 (MQ=255)
           gACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAac                                                                                                                                   <  1:1270642/150‑1 (MQ=255)
                          gacgacGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACgg                                                                                                                    >  2:1172493/1‑150 (MQ=255)
                             gacgCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCg                                                                                                                 >  2:554673/1‑150 (MQ=255)
                                     tcgCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCAtt                                                                                                         <  1:554673/150‑1 (MQ=255)
                                                     gccgaggccgagCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtc                                                                                         <  1:589917/150‑1 (MQ=255)
                                                                     ggccggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCttt                                                                         >  2:381763/1‑150 (MQ=255)
                                                                       ccggccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTc                                                                       >  2:674361/1‑150 (MQ=255)
                                                                          gccAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACt                                                                    <  2:828320/150‑1 (MQ=255)
                                                                               ggCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCtt                                                                          <  1:790930/139‑1 (MQ=255)
                                                                               ggCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTATCGTCATTTAGCCtt                                                                          >  2:790930/1‑139 (MQ=255)
                                                                                   gaggagGTAACGCCGCCCTTGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCa                                                                                   <  2:231135/126‑1 (MQ=255)
                                                                                   gaggagGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCa                                                                                   >  1:231135/1‑126 (MQ=255)
                                                                                      gaggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAg                                                        >  2:967421/1‑150 (MQ=255)
                                                                                       aggTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAgg                                                       <  1:1172493/150‑1 (MQ=255)
                                                                                           aaCGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATa                                                   >  1:1416332/1‑150 (MQ=255)
                                                                                            aCGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATaa                                                  >  2:879771/1‑150 (MQ=255)
                                                                                              gccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAgg                                                <  1:381763/150‑1 (MQ=255)
                                                                                              gccgccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTGA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAgg                                                >  1:982815/1‑150 (MQ=255)
                                                                                                 gccCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGtat                                             >  2:810088/1‑150 (MQ=255)
                                                                                                  cccTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTCGCCTTTTCTCTGGCGCGCGACAGGCATAAGGtata                                            <  2:237446/150‑1 (MQ=255)
                                                                                                          gCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGGCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTGCACTGGCGCGCGGCAGCCATAAGGTATGCGCCATTc                                    >  1:498528/1‑150 (MQ=255)
                                                                                                             gtagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                       >  1:427792/1‑96 (MQ=255)
                                                                                                             gtagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                       <  2:427792/96‑1 (MQ=255)
                                                                                                              tagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGgc                                                                 >  1:997029/1‑117 (MQ=255)
                                                                                                              tagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGgc                                                                 <  2:997029/117‑1 (MQ=255)
                                                                                                              tagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGtt                                     >  1:641748/1‑145 (MQ=255)
                                                                                                              tagtTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGtt                                     <  2:641748/145‑1 (MQ=255)
                                                                                                                   tATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCt                           <  1:967421/150‑1 (MQ=255)
                                                                                                                        cTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGc                      >  2:868956/1‑150 (MQ=255)
                                                                                                                              gggctctTCGC‑GGGTAACCCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                       >  1:1425105/4‑79 (MQ=25)
                                                                                                                              ggCCTCTTCGCGGGTAAC‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATtctc                                                                                       <  2:1425105/78‑1 (MQ=25)
                                                                                                                               gCCTCTTCGCGGGTAAC‑CCCGCGAAGAGGCAAACACAGGCAAACCACGGCAGTTGCCATTGACGCCAGACACATTCTCGTAATTTAGCCTTTTAAGTGGCGCGCGACAGGCGTAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATg               >  2:360532/1‑150 (MQ=255)
                                                                                                                               gCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTgg                                                                  <  2:803376/99‑1 (MQ=255)
                                                                                                                               gCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTgg                                                                  >  1:803376/1‑99 (MQ=255)
                                                                                                                                 ctctTCGCGGGTAAC‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTc             >  2:523212/1‑150 (MQ=255)
                                                                                                                                  tctTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCAt                >  2:600829/1‑146 (MQ=255)
                                                                                                                                  tctTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCAt                <  1:600829/146‑1 (MQ=255)
                                                                                                                                    ttCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACa                                                         >  2:77119/1‑103 (MQ=255)
                                                                                                                                    ttCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACa                                                         <  1:77119/103‑1 (MQ=255)
                                                                                                                                    ttCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaa          <  2:1426982/150‑1 (MQ=255)
                                                                                                                                     tCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTccaac         <  1:868956/150‑1 (MQ=255)
                                                                                                                                        cgGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGc                                                                                                             >  1:1289070/1‑47 (MQ=255)
                                                                                                                                        cgGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGc                                                                                                             <  2:1289070/47‑1 (MQ=255)
                                                                                                                                            gTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTca  >  1:765109/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TGCGTGAGATCGACGCCGCCCGGATGGACGACGCTCGTCGCCGTGCAGCACGCGCCGAGGCCGAGCATGGGCCGGCCAAGGCAGAGGAGGTAACGCCGCCCTGGCTGCAGTAGTTTATTCCTGTACCGCCTCTTCGCGGGGTTA‑CCCGCGAAGAGGCCAACACAGGCAAACCACGGCCGTTGCCATTGACGCAAGACACATTCTCGTCATTTAGCCTTTTCACTGGCGCGCGACAGGCATAAGGTATACGCCGTTCAAATTCCCTGGAGCTGTCATGTCCAACCAATTCC  >  NC_002947/3312590‑3312879

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