Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,289,861:1 +A 8.9% coding (1128/2349 nt) PP_1126 → hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,289,8611.A 8.9% 161.6 / 9.8 45coding (1128/2349 nt)PP_1126hydrolase
Reads supporting (aligned to +/- strand):  ref base . (24/17);  new base A (2/2);  total (26/19)
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.83e-01

CCGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCGC  >  NC_002947/1289712‑1290005
                                                                                                                                                      |                                                                                                                                                
ccGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGcc                                                                                                                                                   >  2:331013/1‑150 (MQ=255)
      tATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCtcgtcg                                                                                                                                            >  2:676823/1‑150 (MQ=255)
             tGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGgc                                                                                                                                     >  2:656983/1‑150 (MQ=255)
                        gCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAg                                                                                                                          >  1:768380/1‑150 (MQ=255)
                            tGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAgg                                                                                                                      <  1:217645/150‑1 (MQ=255)
                                ggcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtc                                                                                                                                             <  2:153108/123‑1 (MQ=255)
                                ggcAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtc                                                                                                                                             >  1:153108/1‑123 (MQ=255)
                                      cggcggGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTAGGAGATTCTCAATTTCAAGCAGGGGCGCTGGtt                                                                                                            <  1:1107037/150‑1 (MQ=255)
                                                aCCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTAc                                                                                                  >  1:1011539/1‑150 (MQ=255)
                                                 gcTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTAcg                                                                                                 >  2:379800/2‑150 (MQ=255)
                                                        cGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCaa                                                                                          <  1:646086/150‑1 (MQ=255)
                                                         gCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAAc                                                                                         <  1:379800/150‑1 (MQ=255)
                                                         gCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAAc                                                                                         <  1:128840/150‑1 (MQ=255)
                                                             ggCGCGCAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGc                                                                                     >  1:1135940/1‑150 (MQ=255)
                                                                 cgAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAACCTCGGGTGCGGGTTCGACTTCTACGCCGCCAACCTGCAGCg                                                                                 >  2:473655/1‑150 (MQ=255)
                                                                  gAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGt                                                                                >  1:498775/1‑150 (MQ=255)
                                                                  gAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtc                                                                                                                                             >  2:826046/1‑89 (MQ=255)
                                                                  gAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGAGGCATTGCCAtcgtc                                                                                                                                             <  1:826046/89‑1 (MQ=255)
                                                                    aGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTAc                                                                              <  1:331013/150‑1 (MQ=255)
                                                                      gCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAgg                                                                                                                                       >  1:361855/1‑91 (MQ=255)
                                                                      gCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAgg                                                                                                                                       <  2:361855/91‑1 (MQ=255)
                                                                          gTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATgg                                                                        >  2:547478/1‑150 (MQ=255)
                                                                                    ttGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGTGAAAGCGACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGc                                                              <  2:866844/150‑1 (MQ=255)
                                                                                    ttGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATTGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCCACCTGCAGCGTACCTATGGCAAACCCTGc                                                              <  1:495469/150‑1 (MQ=255)
                                                                                       cacGCTGAGGGCGATGCACTGGGCGCTGTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCgg                                                           >  2:313342/3‑150 (MQ=255)
                                                                                          gCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGccc                                                        <  1:676823/150‑1 (MQ=255)
                                                                                                 gCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTgg                                                 >  2:1067601/1‑150 (MQ=255)
                                                                                                  cGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGa                                                <  1:656983/150‑1 (MQ=255)
                                                                                                  cGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGa                                                <  1:586318/150‑1 (MQ=255)
                                                                                                        aCTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGgcca                                          <  1:1067601/150‑1 (MQ=255)
                                                                                                         cTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGTCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTAc                                                                              >  1:853138/1‑113 (MQ=255)
                                                                                                         cTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTAc                                                                              <  2:853138/113‑1 (MQ=255)
                                                                                                              cGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCcagccagc                                    >  2:1386693/1‑150 (MQ=255)
                                                                                                               gCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATgg                                                                        <  1:819354/113‑1 (MQ=255)
                                                                                                               gCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATgg                                                                        >  2:819354/1‑113 (MQ=255)
                                                                                                                         gCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTgcgcgc                         >  2:1054102/1‑150 (MQ=255)
                                                                                                                                caacaaCGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGt                  >  1:1010931/1‑150 (MQ=255)
                                                                                                                                caacaaCGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGt                  >  2:933572/1‑150 (MQ=255)
                                                                                                                                    aacGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTg              >  1:786673/1‑150 (MQ=255)
                                                                                                                                      cgacgaTGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGa            >  2:608935/1‑150 (MQ=255)
                                                                                                                                       gacgaTGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGAt           >  1:122490/1‑150 (MQ=255)
                                                                                                                                           gtGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCg       <  2:122490/149‑1 (MQ=255)
                                                                                                                                               cAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCg   >  2:1334412/1‑150 (MQ=255)
                                                                                                                                                aaTGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCcagccagc                                    <  2:1180302/116‑1 (MQ=255)
                                                                                                                                                aaTGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCcagccagc                                    >  1:1180302/1‑116 (MQ=255)
                                                                                                                                                aaTGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGGCCGc  >  1:948893/1‑150 (MQ=255)
                                                                                                                                                aaTGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTCGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCGc  >  2:843977/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
CCGTTCTATTCCCTGGGGGTGAACGCCGTGGCGGCAGACGGCGGGCGTACCTATGTCGCCGGGCGCGAAGGCATGTTCAAGGCCTTGCCCGCTGAGGGCGATGCACTGGCCGCTTTCTATGGCGAAAGCAACAACGACGATGGCAATGCCTCGTCGCAAGGGCGCAATTTCAAGCAGGGGCGCTGGTTCGACTTCTACGCCGCCAACCTGCAGCGTACCTATGGCAAACCCTGCCGGCCCGCAGTGGAGGGCCAGCCAGCCGATTGCGCGCCGCCGGTGCTGGATGCCGACCGC  >  NC_002947/1289712‑1290005

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: