Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,788,503 G→C 8.7% Y44* (TAC→TAG)  PP_3350 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,788,5030GC 8.7% 140.9 / 3.9 46Y44* (TAC→TAG) PP_3350hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base C (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01

GCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCA  >  NC_002947/3788365‑3788649
                                                                                                                                           |                                                                                                                                                  
gCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGa                                                                                                                                          >  3:320177/1‑149 (MQ=255)
  cAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTa                                                                                                                                                  <  8:242521/139‑1 (MQ=255)
  cAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTa                                                                                                                                                  >  7:242521/1‑139 (MQ=255)
   aaaGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATcc                                                                                                                                       <  7:8048/149‑1 (MQ=255)
      gaagGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGc                                                                                                                                    >  3:9412/1‑149 (MQ=255)
            tttGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAgggg                                                                                                                              >  1:329431/1‑149 (MQ=255)
            tttGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAgggg                                                                                                                              >  7:301335/1‑149 (MQ=255)
                               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                                                                 >  3:210885/1‑149 (MQ=255)
                                                                        tGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAg                                                                 >  4:160955/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                                                                                                      <  7:6907/108‑1 (MQ=255)
                                                                            aGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGctct                                                                                                      >  2:78161/1‑108 (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:119500/1‑149 (MQ=255)
                                                                                             gCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGggcag                                            <  7:11984/149‑1 (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                                             >  3:44005/1‑104 (MQ=255)
                                                                                                                                         g‑gTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGggca                                             <  4:44005/104‑1 (MQ=255)
                                                                                                                                         g‑gTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGt                                 >  3:317730/1‑116 (MQ=255)
                                                                                                                                         g‑gTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGt                                 <  4:317730/116‑1 (MQ=255)
                                                                                                                                         g‑gTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGt                                 <  4:312672/116‑1 (MQ=255)
                                                                                                                                         g‑gTAGTCCTCGATCCAGCGAGGGG‑CCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGt                                 >  3:312672/1‑116 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
GCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCG‑GTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCA  >  NC_002947/3788365‑3788649

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: ‑