Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,015,757 G→C 9.3% H22D (CAT→GAT)  PP_5266 ← acetyl‑CoA hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,015,7570GC 9.3% 125.6 / 4.6 43H22D (CAT→GAT) PP_5266acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (16/23);  new base C (2/2);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

CGCCGCTGCAGGCCATCGCCCAGTGGCGGGCGGCCGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTGGTCA  >  NC_002947/6015609‑6015900
                                                                                                                                                    |                                                                                                                                               
cgccgcTGCAGGCCATCGCCCAGTGGCGGGCGGCCGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCatatgg                                                                                                                                                <  1:628356/150‑1 (MQ=255)
          ggCCATCGCCCAGTGGCGGGCGGCCGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTg                                                                                                                                      >  1:620071/1‑150 (MQ=255)
           gCCATCGCCCAGTGGCGGGCGGCCGAGGGAAAGTGCGGGATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTgg                                                                                                                                     <  2:177154/150‑1 (MQ=255)
              aTCGCCCAGTGGCGGGCGGCCGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGACTCGGGTTCCCCCCGGCCCAAAGGCA‑GCCCATATCCATATGGGcctgcctt                                                                                                                                 <  2:726939/150‑8 (MQ=255)
                              cggcCGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGc                                                                                                                  <  1:519250/150‑1 (MQ=255)
                                 ccGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCt                                                                                                                 <  2:356411/148‑1 (MQ=255)
                                 ccGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCt                                                                                                                 >  1:356411/1‑148 (MQ=255)
                                   gAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGAc                                                                                                             >  2:112289/1‑150 (MQ=255)
                                    aGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACa                                                                                                            >  1:105416/1‑150 (MQ=255)
                                               ggTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGc                                                                                                          >  2:532494/1‑141 (MQ=255)
                                               ggTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGc                                                                                                          <  1:532494/141‑1 (MQ=255)
                                               ggTATAGATGGTCAGGCGCCGGTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGAt                                                                                                 <  1:112289/150‑1 (MQ=255)
                                                  ataGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGa                                                                                              <  2:620071/150‑1 (MQ=255)
                                                        ggTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCagag                                                                                        >  2:392581/1‑150 (MQ=255)
                                                              gcgcCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCgg                                                                                                                              <  2:192027/106‑1 (MQ=255)
                                                              gcgcCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCgg                                                                                                                              >  1:192027/1‑106 (MQ=255)
                                                               tgcCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATg                                                                                 <  1:187536/149‑1 (MQ=255)
                                                               cgcCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATg                                                                                 <  1:498543/150‑1 (MQ=255)
                                                               cgcCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATg                                                                                 <  1:372827/150‑1 (MQ=255)
                                                                 ccGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATggg                                                                               <  1:392581/150‑1 (MQ=255)
                                                                   gTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATgg                                                                                <  1:223127/147‑1 (MQ=255)
                                                                   gTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATgg                                                                                >  2:223127/1‑147 (MQ=255)
                                                                      cGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCTGTCGCAAAGCCAGGCCCAGATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGcctc                                                                          <  2:289842/150‑1 (MQ=255)
                                                                                     cgcgGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTgg                                                           <  1:535866/150‑1 (MQ=255)
                                                                                        gggCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTc                                                                                                                >  1:230507/1‑94 (MQ=255)
                                                                                        gggCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTc                                                                                                                <  2:230507/94‑1 (MQ=255)
                                                                                            gTAAAGCGCGTTGACGAAGGCGTTGGGTTTCTCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCTTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGaaa                                                                    <  1:609557/134‑1 (MQ=255)
                                                                                            gTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGaaa                                                                    >  2:609557/1‑134 (MQ=255)
                                                                                                                            ccAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGt                                                             >  2:165510/1‑109 (MQ=255)
                                                                                                                            ccAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGt                                                             <  1:165510/109‑1 (MQ=255)
                                                                                                                            ccAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGcc                    >  2:281232/1‑150 (MQ=255)
                                                                                                                               ggCGCAAACGCAGGCCCACATGGCTTCGTTCTCGCATCCGGGACAAGACCTTCTCGACAGCCTGGTCGATGTAGCAGAGGTGCATGGGGCGTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAg                 <  2:147737/150‑1 (MQ=255)
                                                                                                                                  ccAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCa                                                      <  2:103811/110‑1 (MQ=255)
                                                                                                                                  ccAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCa                                                      >  1:103811/1‑110 (MQ=255)
                                                                                                                                   cAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCa             >  2:343348/1‑150 (MQ=255)
                                                                                                                                   aaaaGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCa             <  2:375686/149‑1 (MQ=255)
                                                                                                                                    aaaGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCaa            >  2:180128/1‑150 (MQ=255)
                                                                                                                                          cagcCCATATCCATATTGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCgg                              <  1:659287/124‑1 (MQ=255)
                                                                                                                                          cagcCCATATCCATATTGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCgg                              >  2:659287/3‑126 (MQ=255)
                                                                                                                                          aggcCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTg      <  2:684476/150‑1 (MQ=255)
                                                                                                                                           ggcCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGTGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTgg     <  1:281232/150‑1 (MQ=255)
                                                                                                                                            gcCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCTTGCGTTCAGGGGTTGGACATGAGCCCACGCGGGTTGGTTGCCTAGGGCAAACTCTGGt    <  2:105416/150‑1 (MQ=255)
                                                                                                                                            gcCCATATCCATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTGAGGTGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTGGt    >  1:727177/1‑150 (MQ=255)
                                                                                                                                              ccATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTGGTCa  >  2:105838/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CGCCGCTGCAGGCCATCGCCCAGTGGCGGGCGGCCGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTGGTCA  >  NC_002947/6015609‑6015900

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