Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,271,650 G→A 7.8% I145I (ATC→ATT)  PP_2869 ← FMN‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,271,6500GA 7.8% 148.7 / 2.2 51I145I (ATC→ATT) PP_2869FMN‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base A (2/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

ACCGTTGTAGCCATCACTGCGTACGTTCACCCGCGGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGTCTGCCGACTCCATGATCTCCCGCCCAGGGTGAAACAGTTGCGACAGCACCAC  >  NC_002947/3271503‑3271790
                                                                                                                                                   |                                                                                                                                            
aCCGTTGTAGCCATCACTGCGTACGTTCACCCGCGGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGa                                                                                                                                             >  8:129824/1‑149 (MQ=255)
         gCCATCACTGCGTACGTTCACCCGCGGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGcc                                                                                                                                     <  7:137479/148‑1 (MQ=255)
             tCACTGCGTACGTTCACCCGCGGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTc                                                                                                                                <  1:97324/149‑1 (MQ=255)
                       cGTTCACCCGCGGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGt                                                                                                                      <  7:129824/149‑1 (MQ=255)
                         ttCACCCGCGGGTTGTGAAACTGCGCGGGCTGGTAGCCTTGGCTGGCGACCACCTCTTCGCCATCCAGGCCGGCCTGGTGGAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTgg                                                                                                                    <  4:22021/149‑1 (MQ=255)
                         ttCACCCGCGGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTgg                                                                                                                    <  4:142541/149‑1 (MQ=255)
                          tCACCCGCGGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCAATCATTGCCTgg                                                                                                                                  >  1:126498/1‑134 (MQ=255)
                          tCACCCGCGGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCAATCATTGCCTgg                                                                                                                                  <  2:130679/134‑1 (MQ=255)
                          tCACCCGCGGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCAATCATTGCCTgg                                                                                                                                  <  2:126498/134‑1 (MQ=255)
                          tCACCCGCGGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCAATCATTGCCTgg                                                                                                                                  >  1:130679/1‑134 (MQ=255)
                                gcgGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCAtt                                                                                                                                        <  1:76103/122‑1 (MQ=255)
                                gcgGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCAtt                                                                                                                                        >  2:76103/1‑122 (MQ=255)
                                     ttGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGaa                                                                                                        <  5:267108/149‑1 (MQ=255)
                                                     gaaGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTTCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGGTCATTGTGCAcc                                                                                        >  1:60887/3‑149 (MQ=255)
                                                             ccATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGc                                                                                >  5:44088/1‑149 (MQ=255)
                                                                tggctggcCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGcc                                                                               >  5:137783/1‑147 (MQ=255)
                                                                tggctggcCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGa                                                                             <  4:147615/149‑1 (MQ=255)
                                                                 ggctggcCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGAt                                                                                                                                            >  2:35267/1‑85 (MQ=255)
                                                                 ggctggcCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGAt                                                                                                                                            <  1:35267/85‑1 (MQ=255)
                                                                 ggctggcCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAg                                                                            >  4:189450/1‑149 (MQ=255)
                                                                   ctggcCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTa                                                                          <  1:268224/149‑1 (MQ=255)
                                                                   ctggcCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTa                                                                          <  8:108258/149‑1 (MQ=255)
                                                                   ctggcCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCCGAAGCGCTCATTGGGCACCGCAGACGCCGAg                                                                            >  3:48140/1‑147 (MQ=255)
                                                                    tggcCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAg                                                                         <  6:44088/149‑1 (MQ=255)
                                                                        caccacCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTa                                                                     >  5:250413/1‑149 (MQ=255)
                                                                               tCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAgcgc                                                                                                    <  6:123432/111‑1 (MQ=255)
                                                                               tCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAgcgc                                                                                                    >  5:123432/1‑111 (MQ=255)
                                                                               tCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGa                                                                             <  8:154151/134‑1 (MQ=255)
                                                                               tCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGa                                                                             >  7:154151/1‑134 (MQ=255)
                                                                               tCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCaa                                                              >  1:115198/1‑149 (MQ=255)
                                                                                cGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCaa                                                              >  4:244070/1‑148 (MQ=255)
                                                                                   cGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAgg                                                          <  6:137783/149‑1 (MQ=255)
                                                                                    gCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGCCGCCGAGTAGGCTACCGCCAACAGGc                                                         >  2:176897/1‑149 (MQ=255)
                                                                                     ccATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACACGGAAGCGCTCATTGGACACCGCAGACGCCGAGTAGGCTACCGCCAACAGGcc                                                        >  7:89599/1‑149 (MQ=255)
                                                                                       aTCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGt                                                      >  8:107143/1‑149 (MQ=255)
                                                                                        tCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGTc                                                     >  2:91732/1‑149 (MQ=255)
                                                                                            ggccggccTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTa                                                                     <  6:254299/129‑1 (MQ=255)
                                                                                            ggccggccTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTa                                                                     >  5:254299/1‑129 (MQ=255)
                                                                                                                        gcTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTg                                                                                              <  2:257956/76‑1 (MQ=255)
                                                                                                                            gagTAGCCCTGGACGATTTAGTCGATCATTGCCTGGTTCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGTCTGCCGACTCCATGATCTCCCGCCCAGGGTGAAACAg                 >  4:107211/3‑149 (MQ=255)
                                                                                                                             cGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGTCTGCCGACTCCATGATCTCCCGCCCAGGGTGAAACAGt                >  1:13324/1‑149 (MQ=255)
                                                                                                                              gTAGCCCTGCACGATTTCGTCGATCATTGCCTGGGCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTTGGCTACCGCCCACAGTCCGTGTGCCGACTCCATGATCTCCCGCCCAGGGTGAAACAGtt               <  3:15609/149‑1 (MQ=255)
                                                                                                                               tAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGTCTGCCGACTCCATGATCTCCCGCCCAGGGTGAAACAGTTg              <  3:189450/149‑1 (MQ=255)
                                                                                                                                   ccTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGTCTGCCGACTCCATGATCTCCCGCCCAGGGTGAAACAGTTGCGa           >  7:233395/1‑148 (MQ=255)
                                                                                                                                   ccTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGTCTGCCGACTCCATGATCTCCCGCCCAGGGTGAAACAGTTGCGAc          <  2:13324/149‑1 (MQ=255)
                                                                                                                                     tGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGTCTGCCGACTCCATGATCTCCCGCCCAGGGTGAAACAGTTGCGACAg        >  1:257488/1‑149 (MQ=255)
                                                                                                                                       gACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAg                                                                            >  2:218071/1‑79 (MQ=255)
                                                                                                                                       gACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAg                                                                            <  1:218071/79‑1 (MQ=255)
                                                                                                                                          gATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGTCTGCCGACTCCATg                                        <  6:37437/112‑1 (MQ=255)
                                                                                                                                          gATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGTCTGCCGACTCCATg                                        >  5:37437/1‑112 (MQ=255)
                                                                                                                                           aTTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGTCTGCCGACTCCATGATCTCCCGCCCAGGGTGAAACAGTTGCGACAGcaccac  <  2:115198/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
ACCGTTGTAGCCATCACTGCGTACGTTCACCCGCGGGTTGAGAAACTGCGCGGGCAGGTAGCCATGGCTGGCCACCACCTCGACGCCATCCAGGCCGGCCTGGTGCAGGCGGCGGGCGGCGCTGGCGTAGCCCTGGACGATTTCGTCGATCATTGCCTGGTCCAGCGCCCGTGGCATCACCCGGAAGCGCTCATTGGGCACCGCAGACGCCGAGTAGGCTACCGCCAACAGGCCGTCTGCCGACTCCATGATCTCCCGCCCAGGGTGAAACAGTTGCGACAGCACCAC  >  NC_002947/3271503‑3271790

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