Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,015,758 G→C 9.3% I21M (ATC→ATG)  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,7580GC 9.3% 127.8 / 4.4 43I21M (ATC→ATG) 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 = 9.76e-01

CGCCGCTGCAGGCCATCGCCCAGTGGCGGGCGGCCGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTGGTCAACCCCAC  >  NC_002947/6015609‑6015907
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                       <  2:230507/94‑1 (MQ=255)
                                                                                        gggCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTc                                                                                                                       >  1:230507/1‑94 (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)
                                                                                                                                                     gatatggGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTGGTCAACCCCAc  <  1:105838/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCGCTGCAGGCCATCGCCCAGTGGCGGGCGGCCGAGGGAAAGTGCGGTATAGATGGTCAGGCGCCGTTCGGGCAGCTCGCGTACGCGGGCGTAAAGCGCGTTGACGAAGGCGTTGGGTTTCCCCAGGCCCAAAGGCAGGCCCATATGGATATGGGCTGGCAGCCGGGACAAGACCTGCTCGACAGCCTGGTCGATGGAGCAGAGGTGCATGGGGCCTCCTGAAACATCCGTGGGTTCAGGGGTTGGACATGAGCTCAAGCGGGTTGGTTGCCTAGGGCAAACTCTGGTCAACCCCAC  >  NC_002947/6015609‑6015907

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: ‑