Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,337,383 G→A 10.2% G33G (GGG→GGA)  PP_1163 → FAD‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,337,3830GA10.2% 102.1 / 3.5 39G33G (GGG→GGA) PP_1163FAD‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (20/15);  new base A (2/2);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

CGGGCTGCAAACCCACGCCTGAAACCCCGCCATCGCGAGGAACGCCGCCATGCAACAATCTCCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAA  >  NC_002947/1337236‑1337519
                                                                                                                                                   |                                                                                                                                        
cGGTCTGCAAACCCACGCCTGAAACCCCTCCATCGCGAGGAACGCCGCCATTCATCAATCTCCTGATCACCGCCCCTCGCCCTGGCATCCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCggg                                                                                                                                          <  1:120542/148‑1 (MQ=255)
         aaCCCACGCCTGAAACCCCGCCATCGCGAGGAACGCCGCCATGCAACAATCTCCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCa                                                                                                                                <  7:229262/149‑1 (MQ=255)
             cACGCCTGAAACCCCGCCATCGCGAGGAACGCCGCCATGCAACAATCTCCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCg                                                                                                                            >  8:99904/1‑149 (MQ=255)
                         cccGCCATCGCGAGGAACGCCGCCATGCAACAATCTCCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAATAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGcc                                                                                                                >  2:195674/1‑149 (MQ=255)
                          ccGCCATCGCGAGGAACGCCGCCATGCAACAATCTCCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGAc                                                                                                                        <  4:46078/140‑1 (MQ=255)
                          ccGCCATCGCGAGGAACGCCGCCATGCAACAATCTCCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGAc                                                                                                                        >  3:46078/1‑140 (MQ=255)
                            gCCATCGCGAGGAACGCCGCCATGCAACAATCTCCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGa                                                                                                             <  7:99904/149‑1 (MQ=255)
                                   cgAGGAACGCCGCCATGCAACAATCTCCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCAc                                                                                                      >  2:21728/1‑149 (MQ=255)
                                         aCGCCGCCATGCAACAATCTCCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACg                                                                                                >  5:319500/1‑149 (MQ=255)
                                                     aacaaTCTCCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGACAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAg                                                                                    >  6:37782/1‑149 (MQ=255)
                                                     aacaaTCTCCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGCTGTAGCAGAACGCATGGAGGCGTTCGGACAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAg                                                                                    >  4:25174/1‑149 (MQ=255)
                                                           ctcCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCg                                                                              >  8:235031/1‑149 (MQ=255)
                                                              cTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTc                                                                           <  6:319500/149‑1 (MQ=255)
                                                               tGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCa                                                                          >  4:11399/1‑149 (MQ=255)
                                                                 aTCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGACAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCa                                                                          <  3:25174/147‑1 (MQ=255)
                                                                 aTCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGACAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCa                                                                          <  5:37782/147‑1 (MQ=255)
                                                                      cGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTcgc                                                                   >  4:103827/1‑149 (MQ=255)
                                                                              gCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAgtg                                                           <  7:235031/149‑1 (MQ=255)
                                                                                ccTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGtgtg                                                         >  2:304866/1‑149 (MQ=255)
                                                                                      aTGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGcc                                                   >  3:154133/1‑149 (MQ=255)
                                                                                      aTGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGcc                                                   >  4:253590/1‑149 (MQ=255)
                                                                                          gggTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAg                                                >  7:89170/1‑148 (MQ=255)
                                                                                            gTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCACCGCCAGTGTGGATGCCCAGGGt                                             >  8:100221/1‑149 (MQ=255)
                                                                                               aaaaGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCg                                           <  1:21728/148‑1 (MQ=255)
                                                                                                aaaGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGc                                         >  3:108708/1‑149 (MQ=255)
                                                                                                     cccTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTgg                                    <  1:304866/149‑1 (MQ=255)
                                                                                                      ccTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTggg                                   >  4:284017/1‑149 (MQ=255)
                                                                                                          cAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAAGGCCGACCCTGGGcatc                               <  6:177860/149‑4 (MQ=255)
                                                                                                            ggAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACgc                             >  8:279879/1‑149 (MQ=255)
                                                                                                                aaGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACgctgct                         >  8:209181/1‑149 (MQ=255)
                                                                                                                 aGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCACCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACGctgctg                        <  7:100221/149‑1 (MQ=255)
                                                                                                                   gTCGGGGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACGCTGCTgg                       >  7:291675/1‑148 (MQ=255)
                                                                                                                      ggTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACGCTGCTGGAAgg                   <  3:284017/149‑1 (MQ=255)
                                                                                                                        tgtAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGc                 <  3:11399/149‑1 (MQ=255)
                                                                                                                        tgtAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGc                 <  4:154133/149‑1 (MQ=255)
                                                                                                                                       ggAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAgg                                               >  7:96151/1‑104 (MQ=255)
                                                                                                                                       ggAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAgg                                               <  8:96151/104‑1 (MQ=255)
                                                                                                                                        gAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCGGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTaa  <  3:103827/148‑1 (MQ=255)
                                                                                                                                        gAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGCCa               <  2:58760/135‑1 (MQ=255)
                                                                                                                                        gAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGCCa               >  1:58760/1‑135 (MQ=255)
                                                                                                                                                   gCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACgctgc                          <  6:94009/113‑1 (MQ=255)
                                                                                                                                                   gCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACgctgc                          >  5:94009/1‑113 (MQ=255)
                                                                                                                                                   |                                                                                                                                        
CGGGCTGCAAACCCACGCCTGAAACCCCGCCATCGCGAGGAACGCCGCCATGCAACAATCTCCTGATCACCGCCCCTCGCCCTGGCATGCGGGTGAAAAGACCCTGCAGGAAAAGGTCGGTGTAGCAGAACGCATGGAGGCGTTCGGGCAAAAGGTCATTCGTGACCACATGCCCGACCAGCACCGAACGTTCTATCACCAGTTACCGTTCATGGTCGCCGCCAGTGTGGATGCCCAGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAA  >  NC_002947/1337236‑1337519

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