Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,817,899 G→T 9.4% G199C (GGT→TGT)  PP_4242 → conserved protein of unknown function with Rhs repeat‑associated core domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,817,8990GT 9.4% 167.8 / 3.0 57G199C (GGT→TGT) PP_4242conserved protein of unknown function with Rhs repeat‑associated core domain
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base T (2/4);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGGCATCAGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAATATTCGAGTAC  >  NC_002947/4817751‑4818046
                                                                                                                                                    |                                                                                                                                                   
tggtggCATCATGAGCTCTACGTTGTCTTTAGAGTCTTGCGTTCAGCCCCATCGACATCAAATTGTTGATCATTGGCAGGGTTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAg                                                                                                                                                     <  3:78132/149‑1 (MQ=255)
   tggCATCAGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGt                                                                                                                                                   >  6:66562/1‑148 (MQ=255)
     gCATCAGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTg                                                                                                                                                  >  5:190852/1‑147 (MQ=255)
        tCAGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACCTTTGTGAAGCt                                                                                                                                             >  1:281302/1‑149 (MQ=255)
          aGGATCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAATCTTCACCTTTGTGAAGCttg                                                                                                                                           <  5:360271/149‑3 (MQ=255)
           ggAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCgt                                                                                                                                          <  7:190310/149‑1 (MQ=255)
           ggAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCgt                                                                                                                                          <  1:118299/149‑1 (MQ=255)
           ggAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCgt                                                                                                                                          <  5:66562/149‑1 (MQ=255)
            gAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACCTTTGTGAAGCttggc                                                                                                                                         <  2:281302/149‑5 (MQ=255)
              gCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACCTTTGTGAAGCttggcgg                                                                                                                                       <  5:362431/149‑7 (MQ=255)
                  tACGTGGTCTTTTGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTg                                                                                                                                   >  8:138351/1‑149 (MQ=255)
                  tACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTg                                                                                                                                   >  3:245220/1‑149 (MQ=255)
                   aCGTGGTCTTTTGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTg                                                                                                                                   <  7:138351/148‑1 (MQ=255)
                    cGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTAttt                                                                                                                                    <  5:188837/146‑1 (MQ=255)
                    cGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgt                                                                                                                                 <  7:163556/149‑1 (MQ=255)
                       ggTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgtgtg                                                                                                                              >  8:313005/1‑149 (MQ=255)
                       ggTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgt                                                                                                                                 >  3:57978/1‑146 (MQ=255)
                       ggTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgt                                                                                                                                 <  4:57978/146‑1 (MQ=255)
                              agagTGTGTCGTTGAGCTCCATCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGAAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCGGCCAAGCTTCACAAATGTGAAGCTCGTGTATTTGgtgt                                                                                                                               <  7:158830/141‑1 (MQ=255)
                              agagTCAGGAGTTCAGCCCCAGCGACATAAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTAAGGCGTGGCCGAGCGCTATGCGAATGCAGACGCCAAGATTCACAAATGTGAAGCTCGTGTATTTGatgtgagatag                                                                                                                        >  8:158830/1‑142 (MQ=255)
                                agTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCaa                                                                                                                      <  6:190852/148‑1 (MQ=255)
                                 gTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAgg                                                                                                                    <  7:313005/149‑1 (MQ=255)
                                    tGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTAt                                                                                                                                      <  7:112647/128‑1 (MQ=255)
                                    tGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTAt                                                                                                                                      >  8:112647/1‑128 (MQ=255)
                                    tGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGAtt                                                                                                                 >  1:192298/1‑149 (MQ=255)
                                      gCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTc                                                                                                               >  1:14623/1‑149 (MQ=255)
                                         ttCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGAt                                                                                                            >  3:1775/1‑149 (MQ=255)
                                                 cAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTtgt                                                                                                    <  2:14623/149‑1 (MQ=255)
                                                           aaaTTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgtgt                                                                                                                               >  1:294631/1‑112 (MQ=255)
                                                           aaaTTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgtgt                                                                                                                               <  2:294631/112‑1 (MQ=255)
                                                             aTTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCaa                                                                                        <  1:55981/149‑1 (MQ=255)
                                                                gttgATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGAGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCTATGTGCTTGTGGCTGCGGGCAAGa                                                                                      >  2:50339/1‑148 (MQ=255)
                                                                            cAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTg                                                                                                                                                  >  5:201632/1‑76 (MQ=255)
                                                                            cAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTg                                                                                                                                                  <  6:201632/76‑1 (MQ=255)
                                                                                ggTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATatgggat                                                                     <  4:1775/149‑1 (MQ=255)
                                                                                 gTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgt                                                                                                                                 <  5:362971/88‑1 (MQ=255)
                                                                                 gTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGgt                                                                                                                                 >  6:362971/1‑88 (MQ=255)
                                                                                  ttttattttaCACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCaa                                                                                        <  6:166491/128‑1 (MQ=255)
                                                                                       ttttaCACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGaaa                                                               >  4:109554/1‑148 (MQ=255)
                                                                                         ttaCACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTgg                                                                                                  <  4:342320/111‑1 (MQ=255)
                                                                                         ttaCACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTgg                                                                                                  >  3:342320/1‑111 (MQ=255)
                                                                                            cacGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGTGAAAAATGAg                                                         <  1:50339/149‑1 (MQ=255)
                                                                                                    tATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAg                                                 >  7:328193/1‑149 (MQ=255)
                                                                                                           tGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGc                                                       >  3:186085/1‑136 (MQ=255)
                                                                                                           tGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGc                                                       <  4:186085/136‑1 (MQ=255)
                                                                                                                    gcTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAAtt                                 <  4:245220/149‑1 (MQ=255)
                                                                                                                      tATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAAtt                                 <  8:299481/147‑1 (MQ=255)
                                                                                                                      tATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTga                               >  7:299481/1‑149 (MQ=255)
                                                                                                                      tATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGCCACCCAGGGAAAGTTATTTAATTg                                >  6:359728/1‑148 (MQ=255)
                                                                                                                                 aGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCAt                    >  7:162509/1‑149 (MQ=255)
                                                                                                                                 aGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCAt                    >  2:145854/1‑149 (MQ=255)
                                                                                                                                      ccAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTAc                >  7:58393/1‑148 (MQ=255)
                                                                                                                                       cAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTa              <  3:109554/149‑1 (MQ=255)
                                                                                                                                        aaGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAAcc                                                   <  1:337594/111‑1 (MQ=255)
                                                                                                                                        aaGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAAcc                                                   >  2:337594/1‑111 (MQ=255)
                                                                                                                                         aGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAat            <  2:192298/149‑1 (MQ=255)
                                                                                                                                                   aGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAATATTCGAGTac  >  3:124799/1‑149 (MQ=255)
                                                                                                                                                   aGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAATATTCGAGTac  <  1:145854/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGTGGCATCAGGAGCTCTACGTGGTCTTTAGAGTCTGGCGTTCAGCCCCAGCGACATCAAATTGTTGATCATGGGCAGGGGTTTTATTTTACACGCAACTATGGCGTGGCCGAGCGCTATGCGAATGCAGCCGCCAAGCTTCACAAAGGTGAAGCTCGTGTATTTGGTGTGTACATCAAGGATTTCGATGTGCTTGTGGCTGCGGGCAAGATCAAAAAAATATGGGATGGGAAAAATGAGGCAACCAGGGAAAGTTATTTAATTGAAGAACATGCATTTACTAATATTCGAGTAC  >  NC_002947/4817751‑4818046

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