Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,391,520 T→C 8.5% W244R (TGG→CGG) ‡ PP_2994 → FMN‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,391,5200TC 8.5% 128.1 / 2.9 48W244R (TGG→CGG) ‡PP_2994FMN‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (25/18);  new base C (2/2);  total (28/20)
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.84e-01

TCTGGAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTC  >  NC_002947/3391372‑3391663
                                                                                                                                                     |                                                                                                                                               
tctGGAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGt                                                                                                                                                 >  3:56374/1‑149 (MQ=255)
   ggAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGc                                                                                                                                              >  8:211324/1‑149 (MQ=255)
    gAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCt                                                                                                                                             >  4:319571/1‑149 (MQ=255)
    gAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTACAGCGCGGAGGGTTCCCCACCGAAGACGCGAAAAAC‑GTCGTTGAGTGGc                                                                                                                                              >  6:44348/1‑148 (MQ=255)
            ggCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGAcc                                                                                                                                     <  1:64900/149‑1 (MQ=255)
            ggCGCGCCTGTTGCGCGAACTCGTCAAGGCGGGCCGGGACGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGCCACTCAACGAcc                                                                                                                                     <  4:247838/149‑1 (MQ=255)
                cgcCTGTTGCTCTAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTGATCTCTGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCTAAAAAG‑GTCGTTTAGTGGCTCAACGACCTgg                                                                                                                                  <  4:157353/148‑1 (MQ=255)
                   cTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGt                                                                                                                              <  1:235999/149‑1 (MQ=255)
                      ttGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTa                                                                                                                             <  3:319571/147‑1 (MQ=255)
                          tCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACAtg                                                                                                                       >  1:224337/1‑149 (MQ=255)
                            gAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACAtgg                                                                                                                      <  7:211324/148‑1 (MQ=255)
                              aaTCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACAtggtg                                                                                                                    >  5:134748/1‑148 (MQ=255)
                                  gTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGtt                                                                                                               >  5:329759/1‑149 (MQ=255)
                                      aGGCGGTCCGGGCCGTTGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCg                                                                                                           <  6:134748/149‑1 (MQ=255)
                                        gCGGTCCGGGCCGTGGTATCTCTGGATTATGCAGTGGCGGACAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑CTCGTAGAGAGGCTCAACGACCTGGGCGAAGACATGGTCGAGTTATCgc                                                                                                          >  2:47806/1‑147 (MQ=255)
                                           gTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAgg                                                                                                      >  7:314205/1‑149 (MQ=255)
                                                    gTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACg                                                                                              >  5:123546/1‑148 (MQ=255)
                                                     tGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCTCCGATTTCCAGCGCGGAGGGTTCACAACCGAATAAGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAg                                                                                                   >  7:274722/1‑142 (MQ=255)
                                                     tGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGaa                                                                                            <  5:44348/149‑1 (MQ=255)
                                                     tGGTTTATGTGGATTTTGCAGTTGTGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGTAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTATACATGGTGGAGTTATCGGGAGGCAg                                                                                                   <  8:274722/142‑1 (MQ=255)
                                                        tttCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGAGGAGTTATCGGGAGGCAGCTACGAAGTc                                                                                         >  8:53218/1‑149 (MQ=255)
                                                         ttCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTcc                                                                                        >  6:271594/1‑149 (MQ=255)
                                                             gtgGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGc                                                                                    >  1:40678/1‑149 (MQ=255)
                                                                gATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCAt                                                                                 >  6:154953/1‑149 (MQ=255)
                                                                     tGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCTAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAg                                                                             <  4:56374/148‑1 (MQ=255)
                                                                                tCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GGCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGAcgcgcgc                                                                 <  5:271594/149‑1 (MQ=255)
                                                                                  aaaCTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTCGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGa                                                               >  6:251199/1‑149 (MQ=255)
                                                                                   aactcaactcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGAt                                                              >  1:137329/1‑149 (MQ=255)
                                                                                      tcaactcCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGt                                                           >  2:344251/1‑149 (MQ=255)
                                                                                                        aGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGaa                                         <  5:154953/149‑1 (MQ=255)
                                                                                                          cgcgGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTCTCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGc                                       <  8:314205/149‑1 (MQ=255)
                                                                                                            cgGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCt                                      >  1:318021/1‑148 (MQ=255)
                                                                                                             ggAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTa                                    <  3:91477/149‑1 (MQ=255)
                                                                                                               aGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGTCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTAtt                                  <  6:329759/149‑1 (MQ=255)
                                                                                                               aGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTAtt                                  <  2:40678/149‑1 (MQ=255)
                                                                                                                   ttCACCACCGAAGACGCGAAAAAGTGTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTccg                               >  7:213340/1‑149 (MQ=255)
                                                                                                                       ccaccGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAg                          >  5:24979/1‑149 (MQ=255)
                                                                                                                             aaGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGcc                    >  3:319757/1‑149 (MQ=255)
                                                                                                                               gACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCg                  <  6:123546/149‑1 (MQ=255)
                                                                                                                                    gAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAtc             >  6:188397/1‑149 (MQ=255)
                                                                                                                                      cccaG‑GTCGTTGAGCCACTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGAtctc           >  5:324753/4‑149 (MQ=255)
                                                                                                                                          g‑gTCGTTGAGCCACTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCaaa       >  5:300740/1‑149 (MQ=255)
                                                                                                                                            gTCGTTGAGTGGCTCAACGTCGTGGGGGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCaaa       <  2:137329/148‑1 (MQ=255)
                                                                                                                                            gTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAAc      >  6:169127/1‑149 (MQ=255)
                                                                                                                                             tCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAAcc     >  1:182851/1‑149 (MQ=255)
                                                                                                                                             tCGTTGAGCCACTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAAcc     <  6:324753/149‑1 (MQ=255)
                                                                                                                                               gTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGt   >  6:156177/1‑149 (MQ=255)
                                                                                                                                               gTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATCTCGTCGAGTTCGCCCGTGATCTCCAAACCGt   >  6:314822/1‑149 (MQ=255)
                                                                                                                                                ttGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTc  <  4:319757/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TCTGGAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTC  >  NC_002947/3391372‑3391663

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