Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,975,166 A→C 11.1% intergenic (+87/+20) PP_3504 → / ← PP_3505 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,975,1660AC11.1% 144.2 / 11.6 54intergenic (+87/+20)PP_3504/PP_3505hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base C (3/3);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

CGAATTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATC  >  NC_002947/3975028‑3975312
                                                                                                                                          |                                                                                                                                                  
cGAATTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACACTTTGTCGCaa                                                                                                                                          >  4:173268/1‑149 (MQ=255)
cGAATTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTc                                                                                                                                              <  8:37381/145‑1 (MQ=255)
cGAATTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTc                                                                                                                                              >  7:37381/1‑145 (MQ=255)
cGAATTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCaa                                                                                                                                          <  4:330494/149‑1 (MQ=255)
  aaTTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACACTTTGTCGCaatt                                                                                                                                        <  3:173268/149‑3 (MQ=255)
   aTTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACACTTTGTCGCaattg                                                                                                                                       <  3:263765/149‑4 (MQ=255)
      gtgtgAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGc                                                                                                                                     <  4:143686/148‑1 (MQ=255)
      gtgtgAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGc                                                                                                                                     <  8:300610/148‑1 (MQ=255)
       tgtgAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGc                                                                                                                                     >  4:352113/1‑147 (MQ=255)
       tgtgAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCa                                                                                                                                   <  4:115817/149‑1 (MQ=255)
        gtgAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCa                                                                                                                                   <  3:69969/148‑1 (MQ=255)
           aaaaGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGtt                                                                                                                                >  4:71640/1‑148 (MQ=255)
           aaaaGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGtt                                                                                                                                <  3:352113/148‑1 (MQ=255)
              aGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAAt                                                                                                                            <  3:71640/149‑1 (MQ=255)
                        ccTGCAACCGTGGTCTATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCt                                                                                                                  <  1:297913/149‑1 (MQ=255)
                        ccTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCt                                                                                                                  >  3:161296/1‑149 (MQ=255)
                         cTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACACCTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCATCCTGCCAGTTCAAGCCATCACCCt                                                                                                                  >  6:41551/1‑148 (MQ=255)
                         cTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCt                                                                                                                  >  5:273666/1‑148 (MQ=255)
                          tGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTg                                                                                                                 >  3:17077/1‑148 (MQ=255)
                                cGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAAttgtt                                                                                                          <  2:30903/149‑1 (MQ=255)
                                           ccGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTgg                                                                                               <  6:273666/149‑1 (MQ=255)
                                                 tAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAAt                                                                                           >  1:262402/1‑147 (MQ=255)
                                                                caacaaGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCt                                                                          >  8:62512/1‑149 (MQ=255)
                                                                 aacaaGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCGACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTg                                                                         >  7:227464/1‑149 (MQ=255)
                                                                  acaaGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGc                                                                        <  3:13993/149‑1 (MQ=255)
                                                                    aaGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGc                                                                        >  7:236751/1‑147 (MQ=255)
                                                                     aGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAAcc                                                                                                                                        >  7:50730/1‑82 (MQ=255)
                                                                     aGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAAcc                                                                                                                                        <  8:50730/82‑1 (MQ=255)
                                                                      gtgcgtgcATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTATGCCAGa                                                                    >  3:195175/1‑149 (MQ=255)
                                                                       tgcgtgcATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCGACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGAc                                                                   <  8:227464/149‑1 (MQ=255)
                                                                                 gTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTgcgcgc                                                          <  2:262402/148‑1 (MQ=255)
                                                                                        aCTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAgg                                                  >  1:44534/1‑149 (MQ=255)
                                                                                            gTACAAGTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCAAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCa                                              >  1:193740/1‑149 (MQ=255)
                                                                                             tACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTg                                                               <  1:301289/131‑1 (MQ=255)
                                                                                             tACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTg                                                               >  2:301289/1‑131 (MQ=255)
                                                                                              aCAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCAc                                             <  3:282121/148‑1 (MQ=255)
                                                                                                aaCTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACggg                                          <  4:161296/149‑1 (MQ=255)
                                                                                                   ttCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAAttgtt                                                                                                          <  8:296243/82‑1 (MQ=255)
                                                                                                   ttCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAAttgtt                                                                                                          >  7:296243/1‑82 (MQ=255)
                                                                                                        tgttgctgttgTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGgc                                                 >  1:186249/1‑134 (MQ=255)
                                                                                                        tgttgctgttgTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGgc                                                 <  2:186249/134‑1 (MQ=255)
                                                                                                            gctgttgTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCa                               >  1:197415/1‑148 (MQ=255)
                                                                                                             ctgttgTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCat                              <  7:62512/148‑1 (MQ=255)
                                                                                                                ttgtGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATc                           <  2:197415/148‑1 (MQ=255)
                                                                                                                    ggCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTGATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGAcc                       >  1:326474/1‑148 (MQ=255)
                                                                                                                       aaCTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACcaac                    >  1:87174/1‑148 (MQ=255)
                                                                                                                           tggcaggttGCGACACTTTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGt               >  2:338879/8‑149 (MQ=255)
                                                                                                                             cccAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTgg             >  8:206386/1‑149 (MQ=255)
                                                                                                                                aaTTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGgcg           >  7:86703/1‑148 (MQ=255)
                                                                                                                                 gttGCGACACTTTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCTATATCGACca                      <  2:108753/135‑1 (MQ=255)
                                                                                                                                 gttGCGACACTTTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCTATATCGACca                      >  1:108753/2‑136 (MQ=255)
                                                                                                                                     cGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACgg      >  2:243354/1‑148 (MQ=255)
                                                                                                                                      gACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCa    >  3:199481/1‑149 (MQ=255)
                                                                                                                                        cAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATc  <  4:17077/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
CGAATTGTGTGAAAAGGGCTACGACCTGCAACCGTGGTCGATGCCGAAGTAATATCGGCCAGCACAACAAGTGCGTGCATGGTCACGCACTTGTACAACTTCTCTGTTGCTGTTGTGGCAACTTTCCCAATTGCGACAAAGTGTCGCAACCTGCCAGTTCAATCCATCACCCTGAATTGTTGCAAGTGCTGGTAATCAGCCTGCAATTGTTCTGCCAGACGCTGCGCGCGTCCCAGGCGCACGGGCGCCAGTTCCATATCGACCAACAACGTGGCGCACGGCATC  >  NC_002947/3975028‑3975312

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: