Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,644,994 A→G 15.2% intergenic (‑46/+20) PP_3211 ← / ← PP_3212 ABC transporter ATP‑binding protein/Rieske 2Fe‑2S family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,644,9940AG15.2% 70.5 / 6.3 33intergenic (‑46/+20)PP_3211/PP_3212ABC transporter ATP‑binding protein/Rieske 2Fe‑2S family protein
Reads supporting (aligned to +/- strand):  ref base A (14/14);  new base G (2/3);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

CCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACAC  >  NC_002947/3644858‑3645138
                                                                                                                                        |                                                                                                                                                
ccTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTggg                                                                                                                                     >  1:294212/1‑150 (MQ=255)
 cTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGc                                                                                                                                    <  1:454909/150‑1 (MQ=255)
       ggCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTgg                                                                                                                                      <  2:863898/142‑1 (MQ=255)
       ggCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTgg                                                                                                                                      >  1:863898/1‑142 (MQ=255)
         cGAACAAAGGGGTGTGCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCa                                                                                                                            <  2:1941/150‑1 (MQ=255)
            aCAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTg                                                                                                                         <  2:23763/150‑1 (MQ=255)
                     tgtCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGaa                                                                                                                <  1:526552/150‑1 (MQ=255)
                         cacaCCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACt                                                                                                            <  2:743473/150‑1 (MQ=255)
                         cacaCCGATGGTGTTCGGCAGGGCCTGCGGAAGCACGCCAGGCACGGATTTACGGTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGACGATGGGGAGGTGTCAGTGGCTGAGGAACACt                                                                                                            >  1:192959/1‑150 (MQ=255)
                            acCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCg                                                                                                         >  1:725114/1‑150 (MQ=255)
                            acCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTACGCCCCGGTGCGCCGGGGcat                                                                                                                                         <  2:537908/118‑3 (MQ=255)
                            acCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTACGCCCCGGTGCGCCGGGGcat                                                                                                                                         >  1:537908/1‑116 (MQ=255)
                                 tGGTGTTCGGCAGGGCCTGGGGCTGGACGTGGGAGTTGGATAGTGTCTGCATGAACATAGCGCGCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACGCTTCGCtctc                                                                                                    <  2:911073/150‑1 (MQ=255)
                                 tGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCtctc                                                                                                    <  2:725114/150‑1 (MQ=255)
                                         ggCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGGTGTCgtgttggg                                                                                            >  2:794047/1‑143 (MQ=255)
                                                           aCGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCagta                                                                          <  1:662258/150‑1 (MQ=255)
                                                                ggCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTg                                                                                                                                       <  1:334327/84‑1 (MQ=255)
                                                                ggCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTg                                                                                                                                       >  2:334327/1‑84 (MQ=255)
                                                                 gCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGc                                                                                                        <  1:329276/114‑1 (MQ=255)
                                                                 gCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGc                                                                                                        >  2:329276/1‑114 (MQ=255)
                                                                                  cATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCg                                                   >  2:986233/1‑150 (MQ=255)
                                                                                   aTCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCCGTAGTTGCCGACGGATAACCAACCCgg                                                  >  2:452564/1‑150 (MQ=255)
                                                                                    tCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCggg                                                 <  1:1015783/150‑1 (MQ=255)
                                                                                        cATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTg                                             >  2:50133/1‑150 (MQ=255)
                                                                                              gcTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAg                                                                                                                   <  2:139401/74‑1 (MQ=255)
                                                                                              gcTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAg                                                                                                                   >  1:139401/1‑74 (MQ=255)
                                                                                                              aTGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCa                       >  2:955846/1‑150 (MQ=255)
                                                                                                                gAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCAcc                     >  2:661614/1‑150 (MQ=255)
                                                                                                                         ccaacgCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGAt                                       >  2:754415/6‑123 (MQ=255)
                                                                                                                         ccaacgCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGAt                                       <  1:754415/118‑1 (MQ=255)
                                                                                                                              gCCCCGGTGCGCCGGGGCATTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTc       <  2:281687/150‑1 (MQ=255)
                                                                                                                                 ccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGgtgccggt    <  1:250521/150‑6 (MQ=255)
                                                                                                                                   ggCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGacac  >  1:978053/1‑150 (MQ=255)
                                                                                                                                        |                                                                                                                                                
CCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACAC  >  NC_002947/3644858‑3645138

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: