Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,419,549 +C intergenic (‑38/+98) 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.C95.8% 156.4 / 0.1 48intergenic (‑37/+99)PP_1242/PP_1244hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (23/23);  total (24/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 = 9.56e-01

CTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAA  >  NC_002947/1419401‑1419676
                                                                                                                                                      |                                                                                                                                  
cTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctggc                                                                                                                                   <  7:161734/149‑1 (MQ=255)
  ggCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctggcct                                                                                                                                 >  1:200422/1‑147 (MQ=255)
       gggATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTcg                                                                                                                            >  7:114203/1‑149 (MQ=255)
          aTGCCTGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCgg                                                                                                                         <  4:119615/149‑1 (MQ=255)
          aTGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCgg                                                                                                                         <  2:288009/149‑1 (MQ=255)
               gggTTCGAAGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCggg                                                                                                                        >  1:331189/1‑145 (MQ=255)
               gggTTCGAAGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCggg                                                                                                                        <  2:331189/145‑1 (MQ=255)
                         cGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTccc                                                                                                           >  7:182225/1‑148 (MQ=255)
                           gAACAGCGGGTGTGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGAGTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGTGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCcaca                                                                                                        <  5:23637/149‑1 (MQ=255)
                           gAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCcaca                                                                                                        <  8:182225/149‑1 (MQ=255)
                           gAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCca                                                                                                          >  8:48209/1‑147 (MQ=255)
                              cAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGt                                                                                                     <  7:48209/149‑1 (MQ=255)
                                  gggTGAGTGTTGTTTTCGCCTTGGTAGTAGCAGGTCTTGCGGACTGTGGTGGGGTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctggc                                                                                                                                   <  3:192784/115‑1 (MQ=255)
                                    gTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtgtg                                                                                               <  2:200422/149‑1 (MQ=255)
                                      gAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtgtg                                                                                               >  6:291416/1‑147 (MQ=255)
                                      gAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCa                                                                                             <  5:414446/149‑1 (MQ=255)
                                      gAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGAGGGCACGCCCGCTCCCACAGGTACTGTGCa                                                                                             >  3:286804/1‑149 (MQ=255)
                                       aGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTACCACAGGTACTGTGCAg                                                                                            >  3:178430/1‑149 (MQ=255)
                                         tgttgtTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCCGGCACGCCCGCTCCCACAGGTACTGTGCAGTg                                                                                          >  6:288926/1‑149 (MQ=255)
                                             gttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCAGTGCGt                                                                                       >  8:285315/1‑148 (MQ=255)
                                                     cTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgt                                                                              >  8:342334/1‑149 (MQ=255)
                                                          tAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCAc                                                                                                                     <  5:190650/105‑1 (MQ=255)
                                                          tAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCAc                                                                                                                     >  6:190650/1‑105 (MQ=255)
                                                            gatgCAGGTCTTGCGGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGGATTGTTTTGTTTCCCTGGTTTCTGGGTTGCCTGGGCTGGCCTCTTGGCTGGCACGCGCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTTGTGAc                                                                        <  4:286804/145‑1 (MQ=255)
                                                            gAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGAcc                                                                       >  1:91362/1‑149 (MQ=255)
                                                               gCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtg                                                                    >  8:100526/1‑149 (MQ=255)
                                                                          cGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTgg                                                                 >  4:273887/1‑141 (MQ=255)
                                                                          cGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTgg                                                                 <  3:273887/141‑1 (MQ=255)
                                                                           gACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCg                                                        <  7:342334/149‑1 (MQ=255)
                                                                            aCTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGATCGGGCGt                                                       <  5:288926/149‑1 (MQ=255)
                                                                              tgtgGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAgg                                                                                                      <  7:136355/100‑1 (MQ=255)
                                                                              tgtgGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAgg                                                                                                      >  8:136355/1‑100 (MQ=255)
                                                                                  gtggGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtgtg                                                                                               >  2:273907/1‑103 (MQ=255)
                                                                                       ttCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtgt                                                                   >  5:36054/1‑126 (MQ=255)
                                                                                       ttCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtgt                                                                   <  6:36054/126‑1 (MQ=255)
                                                                                         cttttgtgttgtCGGTGGGCTTTGTTTTTTTTCCGGGTTTTCTGGGTTGCCTGGGCGGGCCTCTTCGGGGGCACGCCGGCTCCCACTGTTAGTGTGCTTTGGGTGAGGGCTGTGGTGACCGTGTGGTTGGGGGCGTGCCCGCGAagag                                          <  4:178430/145‑1 (MQ=255)
                                                                                          tgttgtgttgTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGc                                        >  7:302617/1‑149 (MQ=255)
                                                                                           gttgtgttgtGGGTGGGATTTTTTTTTTTCCCGTGTATTGTGGGTTGCCTGGGCGGCCCTCTGCGCGGCCACGCCCGCTCCCACTGTTAGTGTGCTGGGCGTTACGGCTGGGGTTACCGTGTGGTTGGGGGCGTGCCCGCGAAGAGGc                                        <  2:91362/148‑1 (MQ=255)
                                                                                              gtgttgtCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACCCCCCCGACCACAGGTACTGTGGAGTGGGTGGCGGCTGTTGTGGCCCGGTGGGAGCGGGCGTGGCCCCGAAGAGGCCTGa                                    >  6:352324/1‑149 (MQ=255)
                                                                                                    tCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtgt                                                                   <  6:23527/113‑1 (MQ=255)
                                                                                                    tCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtgt                                                                   >  5:23527/1‑113 (MQ=255)
                                                                                                             aTTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGa                                                                         <  1:169261/98‑1 (MQ=255)
                                                                                                             aTTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGa                                                                         >  2:169261/1‑98 (MQ=255)
                                                                                                             aTTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAgg                     >  5:207783/1‑149 (MQ=255)
                                                                                                                  tttgtttCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCCCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGgtgt                  >  2:295898/1‑147 (MQ=255)
                                                                                                                   ttgtttCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGa                <  7:285315/148‑1 (MQ=255)
                                                                                                                    tgtttCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAg               <  7:100526/148‑1 (MQ=255)
                                                                                                                           cTGGATTCTGGGTTGCCTGGGCTGGCATCTTCGCGGGCACGCCCGCTCCCCCAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCAGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGGGAGGCTCtgtg       >  5:72159/1‑149 (MQ=37)
                                                                                                                               gATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTgg      <  5:291416/147‑1 (MQ=255)
                                                                                                                                  tCTGGGTTGCCTGGGGTGGCCTGTTCGCGGTCACGCCCGCTCCCACTGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCGTGTGGTTGGGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTaa  <  8:114203/148‑1 (MQ=255)
                                                                                                                                   cTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCTCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTaa  <  1:295898/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                  
CTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAA  >  NC_002947/1419401‑1419676

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: