Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,597,416 G→T 7.8% E67* (GAA→TAA)  PP_2274 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,597,4160GT 7.8% 153.4 / 4.2 51E67* (GAA→TAA) PP_2274hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/27);  new base T (2/2);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

CCTTTCAAAATCCAATCGGACTTCGCTCGCTTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGAGTATCCGGCACGTCCGCTTGATGCGGAACTGGGCGCTGA  >  NC_002947/2597269‑2597562
                                                                                                                                                   |                                                                                                                                                  
ccTTTCAAAATCCAATCGGACTTCGCTCGCTTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGa                                                                                                                                                   <  8:214995/149‑1 (MQ=255)
     cAAAATCCAATCGGACTTCGCTCGCTTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAAc                                                                                                                                               >  3:42560/1‑148 (MQ=255)
          tCCAATCGGACTTCGCTCGCTTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTAc                                                                                                                                         <  5:44438/149‑1 (MQ=255)
             aaTCGGACTTCGCTCGCTTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCtg                                                                                                                                      >  1:93494/1‑149 (MQ=255)
              aTCGGACTTCGCTCGCTTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCtgt                                                                                                                                     <  2:352543/149‑1 (MQ=255)
                       cgctcgctTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCAt                                                                                                                             <  4:340293/148‑1 (MQ=255)
                        gctcgctTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCAtca                                                                                                                           <  2:263974/149‑1 (MQ=255)
                          tcgctTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTAc                                                                                                                                         <  8:339799/133‑1 (MQ=255)
                          tcgctTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTAc                                                                                                                                         >  7:339799/1‑133 (MQ=255)
                            gctTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCAtcactt                                                                                                                        <  8:334845/148‑1 (MQ=255)
                             ctTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGTAAAACGATTAc                                                                                                                                         >  4:151404/1‑130 (MQ=255)
                             ctTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGTAAAACGATTAc                                                                                                                                         >  4:145812/1‑130 (MQ=255)
                             ctTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGTAAAACGATTAc                                                                                                                                         <  3:151404/130‑1 (MQ=255)
                             ctTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGTAAAACGATTAc                                                                                                                                         <  3:145812/130‑1 (MQ=255)
                             ctTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCAtcacttca                                                                                                                      >  7:203651/1‑149 (MQ=255)
                                cGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCAtcacttcac                                                                                                                     <  5:55956/147‑1 (MQ=255)
                                 gAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGCCCCCAtcacttccct                                                                                                                    >  6:358104/1‑147 (MQ=255)
                                 gAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTaa                                                                                                                  <  4:345795/149‑1 (MQ=255)
                                    gCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGCCCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGt                                                                                                                >  8:26617/1‑148 (MQ=255)
                                       cAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCa                                                                                                              >  3:190187/1‑147 (MQ=255)
                                         ggCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCagaga                                                                                                          <  4:190187/149‑1 (MQ=255)
                                           cAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAgg                                                                                                        <  3:313191/149‑1 (MQ=255)
                                           cAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAgg                                                                                                        <  5:83813/149‑1 (MQ=255)
                                                     gctggcatcgctggGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATg                                                                                              <  3:314828/149‑1 (MQ=255)
                                                      ctggcatcgctggGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGt                                                                                             >  6:381402/1‑149 (MQ=255)
                                                         gcatcgctggGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGtt                                                                                            <  5:381402/147‑1 (MQ=255)
                                                         gcatcgctggGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGtt                                                                                            <  8:203651/147‑1 (MQ=255)
                                                                  ggCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTg                                                                                   <  7:26617/147‑1 (MQ=255)
                                                                   gCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCCCTCCCATGTTGAAACCCTGGac                                                                                >  7:46931/1‑149 (MQ=255)
                                                                       gCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGaccat                                                                             <  4:42560/148‑1 (MQ=255)
                                                                         tCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTACCATGTTGAAACCCTGGaccatcac                                                                          >  5:339431/1‑149 (MQ=255)
                                                                                       ggTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACg                                                            <  2:93494/149‑1 (MQ=255)
                                                                                        gTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACgg                                                           >  1:86158/1‑149 (MQ=255)
                                                                                         tCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGc                                                          <  3:211195/149‑1 (MQ=255)
                                                                                              gCCAGTATGTTCGCAAGTGGCTTGTCACCGGCACCGGCCTTGACATTCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCTTCACCATCGACGTGAACGGCACCaa                                                     <  5:80225/149‑1 (MQ=255)
                                                                                               ccAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAgg                                                                                                        <  7:134544/97‑1 (MQ=255)
                                                                                               ccAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAgg                                                                                                        >  8:134544/1‑97 (MQ=255)
                                                                                                cAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAg                                                    >  7:24325/1‑148 (MQ=255)
                                                                                                cAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGc                                                   >  3:304306/1‑149 (MQ=255)
                                                                                                   tATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACg                                                >  2:303019/1‑149 (MQ=255)
                                                                                                          gCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGa                                         <  5:139468/149‑1 (MQ=255)
                                                                                                          gCAAGTGGCGTGTCACCCGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTgg                                          <  2:360014/148‑1 (MQ=255)
                                                                                                                          cggcaccggcCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGAGTATCCGGCACGTCCg                         >  3:13837/1‑149 (MQ=255)
                                                                                                                           ggcaccggcCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGAGTATCCGGCACGTCCGc                        <  2:86158/149‑1 (MQ=255)
                                                                                                                             caccggcCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGAGTATCCGGCACGTCCGCtt                      >  5:67104/1‑149 (MQ=255)
                                                                                                                             caccggcCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGAGTATCCGGCACGTCCGCtt                      >  5:72462/1‑149 (MQ=255)
                                                                                                                               ccggcCTTGACCTGCTAAGTGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTTGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGAGTATCCGGCACGTCCGCTTGa                    <  2:290568/149‑1 (MQ=255)
                                                                                                                                   ccTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGAGTATCCGGCACGTCCGCTTGATGCg                <  4:13837/149‑1 (MQ=255)
                                                                                                                                   ccTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGGCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGAGTATCCGGCACGTCCGCTTGATGCg                <  8:202444/149‑1 (MQ=255)
                                                                                                                                        aCCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGAGTATCCGGCACGTCCGCTTGATGCGGAACt           >  7:320988/1‑149 (MQ=255)
                                                                                                                                          cTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGAGTATCCGGCACGTCCGCTTGATGCGGAACTgg         >  3:125525/1‑149 (MQ=255)
                                                                                                                                                 gggAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCGCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGAGTATCCGGCACGTCCGCTTGATGCGGAACTGGGCGCTGa  <  6:339431/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCTTTCAAAATCCAATCGGACTTCGCTCGCTTCGAAGCCCAGGCAGTAGCATCGCTGGCATCGCTGGGCTTGCTCAGCACCCTGGAAGGTCACGGCCAGTATGGTCGCAAGTGGCGTGTCACCGGCACCGGCCTTGACCTGCTAAGGGAAAACGATTACCTGTGACCCCATCACTTCACTAAGTCAGAGAGGCACTCCCATGTTGAAACCCTGGACCATCACCATCGACGTGAACGGCACCAAGCACGCATTGGAGTATCCGGCACGTCCGCTTGATGCGGAACTGGGCGCTGA  >  NC_002947/2597269‑2597562

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