Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,391,521 G→C 8.3% W244S (TGG→TCG) ‡ PP_2994 → FMN‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,391,5210GC 8.3% 145.2 / 4.9 48W244S (TGG→TCG) ‡PP_2994FMN‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (26/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 = 8.73e-01

GGAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTCGCCAAAA  >  NC_002947/3391375‑3391670
                                                                                                                                                   |                                                                                                                                                     
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)
                                                                                                                                                   ggCTCAACGACCTGGGCGTAGACAT‑GTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTCGCCaaaa  <  6:119030/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGAAAACCGGGCGCGCCTGTTGCTCGAAATCGTCAAGGCGGTCCGGGCCGTGGTTTCTGTGGATTTTGCAGTGGCGGTCAAACTCAACTCCGCCGATTTCCAGCGCGGAGGGTTCACCACCGAAGACGCGAAAAAG‑GTCGTTGAGTGGCTCAACGACCTGGGCGTAGACATGGTGGAGTTATCGGGAGGCAGCTACGAAGTCCCCGCCATGCAGGGTGACGCGCGCGATGGTCGCACCCTGGCTCGCGAAGCTTATTTCGTCGAGTTCGCCCGTGATCTCCAAACCGTCGCCAAAA  >  NC_002947/3391375‑3391670

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