Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,788,526 Δ1 bp 100% coding (109/1392 nt) PP_3350 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,788,5260C.92.5% 204.3 / 8.7 53coding (109/1392 nt)PP_3350hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base . (25/24);  total (27/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 = 9.58e-01
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTT  >  NC_002947/3788396‑3788670
                                                                                                                                   |                                                                                                                                                
cGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCccggccgg                                                                                                                               <  2:33480/149‑1 (MQ=255)
 ggTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCcggccggc                                                                                                                              <  4:9412/149‑1 (MQ=255)
       gggCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGctctct                                                                                                                         >  6:76714/1‑148 (MQ=255)
        ggCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTc                                                                                                                       <  2:329431/149‑1 (MQ=255)
              aGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCgg                                                                                                                 <  4:320177/149‑1 (MQ=255)
                                       tctGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGca                                                                                         >  6:92089/1‑148 (MQ=255)
                                         tGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAATGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCg                                                                                                 <  2:175084/138‑1 (MQ=255)
                                         tGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAGGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCg                                                                                                 >  1:175084/1‑138 (MQ=255)
                                         tGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAg                                                                                      >  4:160955/1‑149 (MQ=255)
                                         tGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAg                                                                                      >  3:210885/1‑149 (MQ=255)
                                           aCAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGcaca                                                                                       <  5:92089/146‑1 (MQ=255)
                                             aGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGCGTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGctct                                                                                                                           <  1:78161/109‑1 (MQ=255)
                                             aGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGctct                                                                                                                           >  8:6907/1‑108 (MQ=255)
                                             aGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGctct                                                                                                                           >  2:78161/1‑108 (MQ=255)
                                             aGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGctct                                                                                                                           <  7:6907/108‑1 (MQ=255)
                                                      tAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑CTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCAt                                                                          >  2:158326/1‑149 (MQ=255)
                                                      tAGCGCACGACATCGAAAGGCGAAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑CTAGTCCTCGATCCACCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCAGGCCGCGTAGCCGGACTGGGACAGCGAAAGCCCCAt                                                                          >  8:189641/1‑149 (MQ=255)
                                                           cAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAgggg                                                                    <  8:301335/149‑1 (MQ=255)
                                                              gCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGggcag                                                                 >  8:11984/1‑149 (MQ=255)
                                                              gCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGggcag                                                                 <  7:11984/149‑1 (MQ=255)
                                                              gCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGggcag                                                                 >  7:119500/1‑149 (MQ=255)
                                                                cATCGAAAGGGTCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑CTAGTCCTCTATCCTGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCTCAGCGAAAGCACCATGAGGggcagg                                                                <  7:189641/149‑1 (MQ=255)
                                                                cATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑CTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGggcagg                                                                <  1:158326/149‑1 (MQ=255)
                                                                    gAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATa                                                            <  8:119500/148‑1 (MQ=255)
                                                                     aaaGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATAcc                                                          >  7:296586/1‑149 (MQ=255)
                                                                        ggggCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATAc                                                           >  2:338280/1‑145 (MQ=255)
                                                                        ggggCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATAc                                                           <  1:338280/145‑1 (MQ=255)
                                                                         gggCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCt                                                                                                              >  2:77629/1‑93 (MQ=255)
                                                                         gggCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCt                                                                                                              <  1:77629/93‑1 (MQ=255)
                                                                             aGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTccc                                                                                                                                      >  4:344949/1‑65 (MQ=255)
                                                                             aGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTccc                                                                                                                                      <  3:344949/65‑1 (MQ=255)
                                                                               tGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAAc                                                 <  8:296586/148‑1 (MQ=255)
                                                                                 gCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACgg                                               <  3:160955/148‑1 (MQ=255)
                                                                                      ggCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACAATGAGGGGCAGGCATACCCCGTTCAACGGCTTGa                                          <  5:76714/148‑1 (MQ=255)
                                                                                                tcgAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAAtt                               >  1:67034/1‑149 (MQ=255)
                                                                                                tcgAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAAtt                               >  5:7227/1‑149 (MQ=255)
                                                                                                    gaAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCAAGGCCGGCTCTCTTAAAGAGGGCTGGCCGCGGAGGCGGAATGGCACAGCGAAAGCACCATCAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAAATTacc                           >  7:91259/1‑149 (MQ=255)
                                                                                                     aaaGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCt                          <  4:210885/149‑1 (MQ=255)
                                                                                                       aGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGAGAAAGCACCATGAGTGGCAGGCATACCCCGTTCAACGGCTTGAAAATCAGAATTTCCCtt                         >  8:231380/1‑148 (MQ=255)
                                                                                                         gg‑gTAGTCCTCGATCCAGCGAGGGG‑CCCTGTCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAATCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCa                       <  7:231380/147‑1 (MQ=255)
                                                                                                          g‑gTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGggca                                                                  <  4:44005/104‑1 (MQ=255)
                                                                                                          g‑gTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGggca                                                                  >  3:44005/1‑104 (MQ=255)
                                                                                                          g‑gTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGt                                                      >  3:317730/1‑116 (MQ=255)
                                                                                                          g‑gTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGt                                                      >  3:312672/1‑116 (MQ=255)
                                                                                                          g‑gTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGt                                                      <  4:312672/116‑1 (MQ=255)
                                                                                                          g‑gTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGt                                                      <  4:317730/116‑1 (MQ=255)
                                                                                                                tCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGttt                >  1:100889/1‑149 (MQ=255)
                                                                                                                  cTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGtttt               >  5:36726/1‑148 (MQ=255)
                                                                                                                  cTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGtttt               >  7:238410/1‑148 (MQ=255)
                                                                                                                   tCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTta             >  6:338926/1‑149 (MQ=255)
                                                                                                                     gATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTtat            >  5:18580/1‑148 (MQ=255)
                                                                                                                      aTCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATaa          <  2:67034/149‑1 (MQ=255)
                                                                                                                           gCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCGAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAgg                                                                      <  8:70845/84‑1 (MQ=255)
                                                                                                                               aggggCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCtt  <  6:7227/148‑1 (MQ=255)
                                                                                                                               aggggCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCt   >  4:110485/2‑148 (MQ=255)
                                                                                                                                   |                                                                                                                                                
CGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTT  >  NC_002947/3788396‑3788670

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