Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,419,548:1 +C 45.8% intergenic (‑37/+99) PP_1242 ← / ← PP_1244 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,419,5481.C45.8% ‑4.2 / 107.4 48intergenic (‑37/+99)PP_1242/PP_1244hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (8/18);  new base C (5/17);  total (13/35)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

CTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCA  >  NC_002947/1419401‑1419687
                                                                                                                                                     |                                                                                                                                            
caggcaGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGC                                                                                                                                             <  8:266686‑M1/147‑1 (MQ=255)
cTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGC                                                                                                                                             <  8:38618‑M1/149‑1 (MQ=255)
cTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGC                                                                                                                                             <  3:290633‑M1/149‑1 (MQ=255)
 tGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTATGGGTTGCCTGGGCTGGCc                                                                                                                                            >  7:7748‑M1/1‑148 (MQ=255)
        ggATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgc                                                                                                                                     >  4:249289‑M1/1‑141 (MQ=255)
          aTGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTGTGTTGTGTTGTCGGTGGTCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGATGGCctcttcgcgg                                                                                                                                   >  1:301618‑M1/1‑139 (MQ=255)
          aTGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgg                                                                                                                                   <  6:142350‑M1/149‑11 (MQ=255)
          aTGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGGTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGGTGCCTGGGCTGGCctcttcgcgg                                                                                                                                   <  3:389814‑M1/149‑11 (MQ=255)
              cGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcac                                                                                                                               >  5:365586‑M1/1‑135 (MQ=255)
              cGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcac                                                                                                                               >  1:380280‑M1/1‑135 (MQ=255)
                  ttCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcacgccc                                                                                                                           <  6:347662‑M1/149‑19 (MQ=255)
                       tGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcgtgcccgcgaa                                                                                                                      <  3:249289‑M1/149‑24 (MQ=255)
                       tGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcacgcccgctcc                                                                                                                      <  8:43730‑M1/149‑24 (MQ=255)
                          ggAACAGCGGGTGAGTGTTGTTTTCGCCTTTGTAGAAGAAGGGATTGCCGACTGTGGTGGGTTCTGTTGTTTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcacgcccgctcccac                                                                                                                   <  2:301618‑M1/149‑27 (MQ=255)
                           gAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcacgcccgctccca                                                                                                                    >  7:247186‑M1/1‑122 (MQ=255)
                                gCGGGTGAGTGTTGTTTTCGCCTTGTTAGAAGAAGGTCTTGCCGACTGTTGTGGGTTCTTTTGGGTTGGCGGTGGGCATTTTTTTGTTTCCCTGGATTCTGGGTTGCCTGTGCTGTCCTCTTCGCGGGCACGCCCGCTCCCACAGGTAc                                                                                                             <  4:155997/149‑1 (MQ=37)
                                gCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcacgcccgctcccacaggtac                                                                                                             <  4:265409‑M1/149‑33 (MQ=255)
                                 cGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcacgcccgctcccacaggtact                                                                                                            <  1:161396‑M1/149‑34 (MQ=255)
                                   ggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcacgcccgctcccacaggtactgt                                                                                                          <  6:365586‑M1/149‑36 (MQ=255)
                                       aGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCa                                                                                                       <  8:329029/148‑1 (MQ=255)
                                          gttgttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGc                                                                                                   <  6:245315/149‑1 (MQ=255)
                                           ttgttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTACCCCCGGTACTGTGCAGTggg                                                                                                  >  8:19500/1‑147 (MQ=255)
                                                     cTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgt                                                                                        <  2:315231/149‑1 (MQ=255)
                                                        ggTAGAAGCATGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgtg                                                                                       <  7:339624/147‑1 (MQ=255)
                                                        ggTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgtg                                                                                       <  7:276680/147‑1 (MQ=255)
                                                          ttgtaGCTGGGCTTGCCGTCTGTGGTGGGTTCTGTTGTTTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCAGGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGa                                                                                   <  4:234008/145‑1 (MQ=38)
                                                                 aGGTCTTGGCGTCTGTGGTGGGTTCTGTTGTGGTGTCGGTGGGCTTTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCGTCTTCGCGGGCACGCCCGCTCCCTCAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtgtg                                                                            <  3:137925/149‑1 (MQ=255)
                                                                 aGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtgtg                                                                            <  5:141460/149‑1 (MQ=255)
                                                                 aGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtgtg                                                                            >  7:156606/1‑149 (MQ=255)
                                                                 aGGTCTTGCCGACTGTGGGGGGTTCTGTTGTGTTGTCGGTGGTCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtgtg                                                                            <  4:331310/149‑1 (MQ=255)
                                                                  ggTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTTGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTgg                                                                           <  8:7748/149‑1 (MQ=255)
                                                                  ggTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcacgcccgctcccacagggtcaccacagccgtcacgcactgcacagtacctgtgg                                                                           <  3:278711‑M1/149‑67 (MQ=255)
                                                                    tCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGTCAGGCctcttcgcgggcacgcccgctcccacagggtcaccacagccgtcacgcactgcacagtacctgtggga                                                                         <  8:195875‑M1/149‑69 (MQ=255)
                                                                        gcggactGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcacgcccgctcccacagggtcaccaaagccgtgacgcactgcacagtacctgtgggagcg                                                                      <  7:99746‑M1/145‑72 (MQ=255)
                                                                        gCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGc                                                                       <  5:253313/147‑1 (MQ=255)
                                                                        gCCGACTGGGGTGGGTTCGTTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGc                                                                       <  8:247186/147‑1 (MQ=38)
                                                                         ccGTCTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCGTGGATTCTGGGGTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCggg                                                                    <  4:48677/149‑1 (MQ=38)
                                                                         ccGACTGGGGGGGGTTCTGTTGTGTTGTCGGGGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCgg                                                                     <  2:380280/148‑1 (MQ=38)
                                                                          cGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGc                                                                   <  4:197985/149‑1 (MQ=255)
                                                                           gACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcacgcccgctcccacagggtcaccacagccgtcacgcactgcacagtacctgtgggagc                                                                       <  4:266669‑M1/144‑71 (MQ=255)
                                                                           gACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcacgcccgctcccacagggtcaccacagccgtcacgcactgcacagtacctgtgggagc                                                                       >  3:266669‑M1/1‑74 (MQ=255)
                                                                                       ttCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttcgcgggcacgcccgctc                                                                                                                       >  4:118255‑M1/1‑62 (MQ=255)
                                                                                       tgcTGTTGTGTTGTGGGTGGGCATTGTTTTGTTTCCGTGGATTCTGGGTGGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCTGTGCGTGAGGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAaga                                                      <  6:134704/147‑1 (MQ=25)
                                                                                          tgttgtgttgTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtgtg                                                                            >  8:79162/1‑124 (MQ=255)
                                                                                          tgttgtgttgTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtgtg                                                                            <  7:79162/124‑1 (MQ=255)
                                                                                                                          ccTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCt                                                                                          <  6:247969/78‑1 (MQ=255)
                                                                                                                          ccTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCt                                                                                          >  5:247969/1‑78 (MQ=255)
                                                                                                                                           cTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCa  >  7:76959/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
CTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCA  >  NC_002947/1419401‑1419687

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: