Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,706,587 T→A 9.8% L202Q (CTG→CAG)  PP_5009 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,706,5870TA 9.8% 166.8 / 8.6 61L202Q (CTG→CAG) PP_5009hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base A (4/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

TCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGT  >  NC_002947/5706443‑5706731
                                                                                                                                                |                                                                                                                                                
tCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCt                                                                                                                                              >  1:60562/1‑149 (MQ=255)
      aaGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAg                                                                                                                                        <  7:77838/149‑1 (MQ=255)
      aaGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAg                                                                                                                                        <  2:60562/149‑1 (MQ=255)
                 gTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAgg                                                                                                                                       >  5:147344/1‑139 (MQ=255)
                 gTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAgg                                                                                                                                       <  6:147344/139‑1 (MQ=255)
                 gTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAg                                                                                                                              >  2:318666/1‑148 (MQ=255)
                    tcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCAc                                                                                                                          >  4:132459/1‑149 (MQ=255)
                        gcATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTc                                                                                                                      <  6:354518/149‑1 (MQ=255)
                              tctGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACg                                                                                                                >  7:60753/1‑149 (MQ=255)
                                  gctggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAATGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGa                                                                                                              <  5:291180/147‑1 (MQ=255)
                                  gctggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATg                                                                                                            >  1:107991/1‑149 (MQ=255)
                                        gCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCgggg                                                                                                      <  5:39443/149‑1 (MQ=255)
                                        gCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCgggg                                                                                                      <  5:243239/149‑1 (MQ=255)
                                              cATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCa                                                                                                                           <  3:361284/122‑1 (MQ=255)
                                              cATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCa                                                                                                                           >  4:361284/1‑122 (MQ=255)
                                              cATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGcacaca                                                                                                >  1:187450/1‑149 (MQ=255)
                                                      cTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTAcc                                                                                        <  5:63599/149‑1 (MQ=255)
                                                         cATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTa                                                                                     <  8:90258/149‑1 (MQ=255)
                                                              cTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGa                                                                                <  2:187450/149‑1 (MQ=255)
                                                                 gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGa                                                                                                                                         >  6:3751/1‑89 (MQ=255)
                                                                 gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGa                                                                                <  6:246490/146‑1 (MQ=255)
                                                                 gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGa                                                                                >  5:246490/1‑146 (MQ=255)
                                                                 gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGAt                                                                               <  2:269215/147‑1 (MQ=255)
                                                                       ctccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCt                                                                       >  6:370684/1‑149 (MQ=255)
                                                                       ctccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCAGGCTGTCgtg                                                                                                                                        >  4:102770/1‑82 (MQ=255)
                                                                       ctccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCAGGCTGTCgtg                                                                                                                                        <  3:102770/84‑3 (MQ=255)
                                                                        tccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCt                                                                       >  3:169445/1‑148 (MQ=255)
                                                                           ctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtcat                                                                    >  2:319401/1‑148 (MQ=255)
                                                                            tcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAAc                                                                                                                 <  1:299992/102‑1 (MQ=255)
                                                                            tcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAAc                                                                                                                 >  2:299992/1‑102 (MQ=255)
                                                                              gCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATcac                                                                 <  8:60753/148‑1 (MQ=255)
                                                                               ccAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATcacca                                                               <  5:370684/149‑1 (MQ=255)
                                                                                         gCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCacacac                                                                                               >  7:174439/1‑107 (MQ=255)
                                                                                         gCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCacacac                                                                                               <  8:174439/107‑1 (MQ=255)
                                                                                          cAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaac                                                    <  7:221002/149‑1 (MQ=255)
                                                                                           aGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaacc                                                   <  1:318666/149‑1 (MQ=255)
                                                                                           aGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaacc                                                   <  2:154460/149‑1 (MQ=255)
                                                                                           aGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaa                                                     <  7:116664/147‑1 (MQ=255)
                                                                                           aGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaa                                                     >  8:116664/1‑147 (MQ=255)
                                                                                                   cAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTg                                           >  5:65920/1‑149 (MQ=255)
                                                                                                      aaGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGtt                                        <  3:132459/149‑1 (MQ=255)
                                                                                                       aGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTc                                       <  5:163932/149‑1 (MQ=255)
                                                                                                            cAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgc                                    <  1:319401/147‑1 (MQ=255)
                                                                                                                    gCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGc                           >  3:42006/1‑148 (MQ=255)
                                                                                                                      ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTAcc                                                                                        <  5:370960/85‑1 (MQ=255)
                                                                                                                      ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTAcc                                                                                        >  6:370960/1‑85 (MQ=255)
                                                                                                                      ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCt                        >  6:240371/1‑149 (MQ=255)
                                                                                                                        tGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGc                      >  1:160024/1‑149 (MQ=255)
                                                                                                                         gCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcg                     >  5:268432/1‑149 (MQ=255)
                                                                                                                          cTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcgg                    <  6:268432/149‑1 (MQ=255)
                                                                                                                          cTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGc                      <  2:160024/147‑1 (MQ=255)
                                                                                                                             gACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcggttc                 <  6:65920/149‑1 (MQ=255)
                                                                                                                                  ggTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCggttcggttc            <  3:192339/149‑1 (MQ=255)
                                                                                                                                    tACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgt          >  7:48693/1‑149 (MQ=255)
                                                                                                                                      ctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtcg        >  3:95653/3‑149 (MQ=255)
                                                                                                                                      ctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtcg        >  1:253353/3‑149 (MQ=255)
                                                                                                                                      ctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtcg        >  8:91219/3‑149 (MQ=255)
                                                                                                                                          aCAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCAt    <  7:91219/149‑1 (MQ=255)
                                                                                                                                           cAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGa                                                                                                              >  7:131711/1‑42 (MQ=255)
                                                                                                                                           cAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGa                                                                                                              <  8:131711/42‑1 (MQ=255)
                                                                                                                                            aGCCTGGCTGTCGAGGCGATACAAGTCACTCTAGGAAAGGATGGCGGGGCACACACCGTTACCCTACCGGATCTAATGCTCATCAACACAAGTGGCAACCGCAACCTGCTTCTGCAATGCAAGCCCTGAGCTTCGGTTCGACGCCATGt  >  3:227564/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGT  >  NC_002947/5706443‑5706731

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