Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,646,064 G→T 9.1% intergenic (‑25/+14) PP_3212 ← / ← PP_3213 Rieske 2Fe‑2S family protein/ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,646,0640GT 9.1% 137.4 / 2.9 44intergenic (‑25/+14)PP_3212/PP_3213Rieske 2Fe‑2S family protein/ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (25/15);  new base T (2/2);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

CTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCAGACCCGCCTCGACCATGAAGTCGCGG  >  NC_002947/3645937‑3646201
                                                                                                                               |                                                                                                                                         
cTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCt                                                                                                                     <  1:671688/150‑1 (MQ=255)
       cagcagGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGAcc                                                                                                              >  2:238935/1‑150 (MQ=255)
           agGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGgt                                                                                                                                >  2:31216/1‑128 (MQ=255)
           agGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGgt                                                                                                                                <  1:31216/128‑1 (MQ=255)
            gggTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGAGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTc                                                                                                         >  1:519880/1‑150 (MQ=255)
             ggTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCg                                                                                                        >  2:789120/1‑150 (MQ=255)
               tGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGc                                                                                                      <  1:238935/150‑1 (MQ=255)
                gAATGGCTTGGGGCCGTCTTCTAGGGATGTCACAGGGAGCAGGGCATACCAGACACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCTCAGCGTTGAGCTTGTTATGCTGTTCAGCGGCTACCGACCTgaacggaa                                                                                                     <  2:278640/150‑8 (MQ=255)
                 aaTGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCt                                                                                                    >  2:970449/1‑150 (MQ=255)
                   tGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCa                                                                                                  <  1:970449/150‑1 (MQ=255)
                         ggggCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGt                                                                                            <  2:374992/150‑1 (MQ=255)
                                 cTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAg                                                                                    <  2:904643/150‑1 (MQ=255)
                                  ttCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGt                                                                                            <  2:469892/141‑1 (MQ=255)
                                  ttCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGt                                                                                            >  1:469892/1‑141 (MQ=255)
                                     gAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGc                                                                                                                   <  2:991753/115‑1 (MQ=255)
                                     gAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGc                                                                                                                   >  1:991753/1‑115 (MQ=255)
                                                tAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAgcg                                                                     >  2:250459/1‑150 (MQ=255)
                                                 aGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAgcgc                                                                    >  2:202972/1‑150 (MQ=255)
                                                     aGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATg                                                                                                                                  >  2:519471/1‑84 (MQ=255)
                                                     aGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATg                                                                                                                                  <  1:519471/84‑1 (MQ=255)
                                                          ggCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAg                                                           >  2:451365/1‑150 (MQ=255)
                                                                 cAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGATCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTc                                                                                 >  1:629618/1‑121 (MQ=255)
                                                                 cAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGATCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTc                                                                                 <  2:629618/121‑1 (MQ=255)
                                                                 cAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGATCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTc                                                                                 >  1:629725/1‑121 (MQ=255)
                                                                 cAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGATCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTc                                                                                 <  2:629725/121‑1 (MQ=255)
                                                                   gAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTgcgc                                                  >  2:74099/1‑150 (MQ=255)
                                                                                aCAGTTTGTTTGATGACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCgg                                     >  2:50224/1‑150 (MQ=255)
                                                                                aCAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCgg                                     <  1:250459/150‑1 (MQ=255)
                                                                                  aGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGa                                   >  1:1240098/1‑150 (MQ=255)
                                                                                  aGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACATGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGa                                   >  2:432497/1‑150 (MQ=255)
                                                                                  aGGTTGTTTGGTAACCAACATAATCAGGATTCATCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGa                                   >  1:1168898/1‑150 (MQ=255)
                                                                                     ttgttTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCg                                >  1:858357/1‑150 (MQ=255)
                                                                                       gttTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACATGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAg                              <  1:432497/150‑1 (MQ=255)
                                                                                         ttGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCa                            <  1:202972/150‑1 (MQ=255)
                                                                                          tGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCCGATCTTGCACAAAGCGCGTGGTGAAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCAg                           >  1:13071/1‑150 (MQ=255)
                                                                                            gTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCa                            >  2:398154/1‑147 (MQ=255)
                                                                                            gTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCa                            <  1:398154/147‑1 (MQ=255)
                                                                                            gTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCAGAc                         <  2:519880/150‑1 (MQ=255)
                                                                                                 cAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCAGACCCGcc                    >  2:672915/1‑150 (MQ=255)
                                                                                                        aTCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCAGACCCGCCTCGACCa             >  2:52609/1‑150 (MQ=255)
                                                                                                        aTCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCAGACCCGCCTCGACCa             >  1:1159178/1‑150 (MQ=255)
                                                                                                         tCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGACTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCAGACCCGCCTCGACCAt            >  1:565868/1‑150 (MQ=255)
                                                                                                             gATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCAGACCCGCCTCGACCATGAAg        >  2:964957/1‑150 (MQ=255)
                                                                                                                   tCAAAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCAGACCCGCCTCGACCATGAAGTCGCgg  >  2:689499/1‑150 (MQ=255)
                                                                                                                               gCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGa                               >  2:116851/1‑109 (MQ=255)
                                                                                                                               gCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGa                               <  1:116851/109‑1 (MQ=255)
                                                                                                                               |                                                                                                                                         
CTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGCTCGGCCTCAGGGAGTACCTTCAGGTCTTGCACAAAGCGCGTGGTGTAGGCGTTGCGCCAAGGTGCGTCGGGACCGAGCAGACCCGCCTCGACCATGAAGTCGCGG  >  NC_002947/3645937‑3646201

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