Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,104,742 T→C 8.1% N48D (AAC→GAC)  PP_1879 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,104,7420TC 8.1% 135.5 / 2.1 49N48D (AAC→GAC) PP_1879hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/22);  new base C (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATG  >  NC_002947/2104596‑2104888
                                                                                                                                                  |                                                                                                                                                   
gATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTGCGGGTGTGAATTTGCCCACGGTTGGATCTCCACCCGCAGGGGGGAATAGGTTGACCGGAACTCGACCGTGTTgg                                                                                                                                                   >  5:299019/1‑149 (MQ=255)
gATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTgg                                                                                                                                                   <  4:77677/149‑1 (MQ=255)
      cTTGCTGAAGGTCTTCTTGTCGACCAGTAAGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCCCGGGTGGACCTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTAGACCGTGTTGGGTTTcc                                                                                                                                             >  8:278973/1‑149 (MQ=255)
         gCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTccc                                                                                                                                            <  1:86725/147‑1 (MQ=255)
           tGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAAGAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCt                                                                                                                                        >  8:51579/1‑149 (MQ=255)
                gTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCaa                                                                                                                                    >  5:53110/1‑148 (MQ=255)
                        gTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTcca                                                                                                                           >  1:113033/1‑149 (MQ=255)
                        gTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTcc                                                                                                                            <  1:29688/148‑1 (MQ=255)
                           gACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTCGGGTTTCCCCCCTTCAACATCGTccacca                                                                                                                        >  7:55213/1‑149 (MQ=255)
                           gACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTCGGGTTTCCCCCCTTCAACATCGTccacca                                                                                                                        >  7:51099/1‑149 (MQ=255)
                            aCCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCcaccac                                                                                                                       >  5:64302/1‑149 (MQ=255)
                              cAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTCGGGTTTCCCCCCTTCAACATCGTCCACCACGc                                                                                                                     <  8:55213/149‑1 (MQ=255)
                              cAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTCGGGTTTCCCCCCTTCAACATCGTCCACCACGc                                                                                                                     <  8:51099/149‑1 (MQ=255)
                                gTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTc                                                                                                                             <  8:232404/139‑1 (MQ=255)
                                gTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTc                                                                                                                             >  7:232404/1‑139 (MQ=255)
                                           cAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTg                                                                                                                  >  4:38876/1‑139 (MQ=255)
                                           cAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTg                                                                                                                  <  3:38876/139‑1 (MQ=255)
                                            aGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggg                                                                                                       <  2:113033/149‑1 (MQ=255)
                                            aGTAGAGGAAAAGGTATTCCTCATCTATATCCGGGTGTTAATTTGCCCACGGGTGGAACTCCACCCGCAGTGGGTAATAGGTTGCCCGGAACTCGAACGTGGTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggg                                                                                                       <  6:299019/149‑1 (MQ=255)
                                             gTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCCGGGGGGAATAGGTTGCCCGGAACTCGACCGGGTTGGGTTTCCCCCCTTCACCAGAGTCCACCAGGCTTGGGTCGGTTgggg                                                                                                      >  6:168953/1‑149 (MQ=255)
                                                  ggAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGAc                                                                                                  >  5:250351/1‑148 (MQ=255)
                                                    aaaaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACaa                                                                                                <  6:53110/148‑1 (MQ=255)
                                                     aaaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGATaa                                                                                                >  1:57391/1‑147 (MQ=255)
                                                     aaaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAaga                                                                                              <  4:118874/149‑1 (MQ=255)
                                                                     tATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGt                                                                               >  7:116094/1‑148 (MQ=255)
                                                                      aTTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGtt                                                                              >  4:128657/1‑148 (MQ=255)
                                                                       tttCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCa                                                                            <  3:128657/149‑1 (MQ=255)
                                                                                  gAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGCCCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGt                                                                >  5:172847/1‑149 (MQ=255)
                                                                                   aaTTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGtt                                                               <  6:64302/149‑1 (MQ=255)
                                                                                    aTTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGAGTTTCCCCCCATCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTc                                                              >  8:235527/1‑149 (MQ=255)
                                                                                      ttGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTg                                                            >  2:87907/1‑149 (MQ=255)
                                                                                      ccccctgcGGGTGGAGTTCCACCCGTAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTg                                                            >  8:13917/8‑149 (MQ=255)
                                                                                        gCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGAGTTTCCCCCCATCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGc                                                           <  7:235527/148‑1 (MQ=255)
                                                                                             cgtggtgAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTAc                                                     <  8:269512/143‑1 (MQ=255)
                                                                                                   gAACTCGACCCGCATGGTGGAATAGGGTGCCCGGAACTCGCCCTTGTTGGGCTTCCCCCCTTCATCATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAgc                                               <  5:168953/149‑1 (MQ=255)
                                                                                                         cACCCGTAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGt                                         <  7:13917/149‑1 (MQ=255)
                                                                                                         cACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAgggg                                                                        >  7:104611/1‑119 (MQ=255)
                                                                                                         cACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAgggg                                                                        <  8:104611/119‑1 (MQ=255)
                                                                                                          aCCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGtt                                        >  4:210524/1‑149 (MQ=255)
                                                                                                              gCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGTTCtgt                                    <  8:116094/149‑1 (MQ=255)
                                                                                                                  gggggAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCACCATCGTCCACCACGCTGGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCa                                >  8:199294/1‑149 (MQ=255)
                                                                                                                    gggAATAGGTTGCCCTGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATc                              <  3:189917/149‑1 (MQ=255)
                                                                                                                    gggAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATc                              <  3:210524/149‑1 (MQ=255)
                                                                                                                      gAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGt                            <  6:250351/149‑1 (MQ=255)
                                                                                                                        aTAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTc                           >  6:154184/1‑148 (MQ=255)
                                                                                                                            ggtGCCCGGAACTCGACCGGGTTGGGTTTCCCCCCTTCAACCTCGTCCACCACGCTTGGTACGGCGGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCAtct                      <  6:172847/147‑1 (MQ=255)
                                                                                                                                          gACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTCTGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCAc           <  1:87907/147‑1 (MQ=255)
                                                                                                                                          gACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCa        >  6:219025/1‑149 (MQ=255)
                                                                                                                                            ccGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCAt       >  7:9498/1‑148 (MQ=255)
                                                                                                                                                 ttGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCTGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATg  <  5:219025/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGT‑TGGAGTTCTGCGGTTACCGGAGCGCTTGTTCTGTGTCATCGTCCATCTCCATGTGTCACGCATGAATG  >  NC_002947/2104596‑2104888

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