Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,391,841 A→T 8.0% K351* (AAA→TAA)  PP_2994 → FMN‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,391,8410AT 8.0% 143.7 / 3.5 50K351* (AAA→TAA) PP_2994FMN‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (24/22);  new base T (2/2);  total (26/24)
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.06e-01

GTCGGCGCCCGGTGGCCGAAAGCGTCGTGCAGAGCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGTTGCGCGCGCTGATCCAGCAACAGCTCTTCATGGCCTGTCGCGCTCG  >  NC_002947/3391698‑3391986
                                                                                                                                               |                                                                                                                                                 
gTCGGCGCCCGGTGGCCGAAAGCGTCGTGCAGAGCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAAc                                                                                                                                              <  7:38237/149‑1 (MQ=255)
         cGGTGGCCGAAAGCGTCGTGCAGAGCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATc                                                                                                                                     <  8:84011/149‑1 (MQ=255)
         cGGTGGCCGAAAGCGTCGTGCAGAGCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATc                                                                                                                                     >  6:222561/1‑149 (MQ=255)
                cGAAAGCGTCGTGCAGAGCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCAACGCATCCt                                                                                                                              >  3:155757/1‑149 (MQ=255)
                 gAAAGCGTCGTGCAGAGCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCtt                                                                                                                             >  2:186954/1‑149 (MQ=255)
                  aaaGCGTCGTGCAGAGCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTg                                                                                                                            <  5:84811/149‑1 (MQ=255)
                          gtGCAGAGCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACa                                                                                                                     <  5:222561/148‑1 (MQ=255)
                            gCAGAGCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAAc                                                                                                                                              <  8:104033/121‑1 (MQ=255)
                            gCAGAGCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAAc                                                                                                                                              >  7:104033/1‑121 (MQ=255)
                             cAGAGCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATgg                                                                                                                  >  2:261183/1‑148 (MQ=255)
                               gagCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCag                                                                                                               >  3:176127/1‑149 (MQ=255)
                                          aCATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTg                                                                                                                            <  3:111303/125‑1 (MQ=255)
                                          aCATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTg                                                                                                                            >  4:111303/1‑125 (MQ=255)
                                           cATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCAACGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTc                                                                                                   <  4:155757/149‑1 (MQ=255)
                                            aTGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAAc                                                                                                                                              >  5:271325/1‑105 (MQ=255)
                                            aTGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAAc                                                                                                                                              <  6:271325/105‑1 (MQ=255)
                                             tggtggGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCaa                                                                                                 <  1:186954/149‑1 (MQ=255)
                                                             ggCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGTAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGcc                                                                                 >  4:219828/1‑149 (MQ=255)
                                                             ggCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGTAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGcc                                                                                 >  4:223514/1‑149 (MQ=255)
                                                               cGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCgg                                                                               <  4:176127/149‑1 (MQ=255)
                                                                      cATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTAGGCGAACATGGAAGTAGTGAAGTTTCAACTGCGCAAACTCCGCCGGGCTcaa                                                                         >  8:100439/1‑148 (MQ=255)
                                                                        tCGATCCGCACTTGCCTCGTGAGTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAg                                                                      >  7:167545/1‑149 (MQ=255)
                                                                        tCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAg                                                                      >  2:47164/1‑149 (MQ=255)
                                                                          gATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGTAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGc                                                                     <  3:219828/148‑1 (MQ=255)
                                                                          gATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGTAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGc                                                                     <  3:223514/148‑1 (MQ=255)
                                                                            tCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCTATTGCCACCGATCACCTGGAAAAACAAAGCCATAGCATCCTTGGCGAACATGGCAGTACTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGccc                                                                   >  1:195925/1‑148 (MQ=255)
                                                                              cGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGcccgac                                                                <  1:47164/149‑1 (MQ=255)
                                                                               gCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGcccgacc                                                               >  8:186659/1‑149 (MQ=255)
                                                                                  cTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGcccgaccc                                                              >  5:69092/1‑147 (MQ=255)
                                                                                     gCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCg                                                         >  1:5377/1‑149 (MQ=255)
                                                                                       cTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGtt                                                       >  8:111792/1‑149 (MQ=255)
                                                                                          gTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGcc                                                                                 >  8:163493/1‑120 (MQ=255)
                                                                                          gTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGcc                                                                                 <  7:163493/120‑1 (MQ=255)
                                                                                          gTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTc                                                     <  2:5377/148‑1 (MQ=255)
                                                                                          gTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCg                                                    >  2:65297/1‑149 (MQ=255)
                                                                                                 gTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGtt                                               <  6:69092/147‑1 (MQ=255)
                                                                                                 gTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGtt                                               <  1:261183/147‑1 (MQ=255)
                                                                                                      aGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTgcgc                                                                                           >  7:233220/1‑98 (MQ=255)
                                                                                                      aGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTgcgc                                                                                           <  8:233220/98‑1 (MQ=255)
                                                                                                      aGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGt                                                <  7:29320/141‑1 (MQ=255)
                                                                                                      aGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGt                                                >  8:29320/1‑141 (MQ=255)
                                                                                                      aGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGTTgcgcgcg                                        <  8:218903/149‑1 (MQ=255)
                                                                                                                  aCGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGTTGCGCGCGCTGATCCAGCa                             >  5:86068/1‑148 (MQ=255)
                                                                                                                         tCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGTTGCGCGCGCTGATCCAGCAACAGctct                     >  8:11975/1‑149 (MQ=255)
                                                                                                                           aaTTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGTTGCGCGCGCTGATCCAGCAACAGCTCTTc                   <  1:65297/149‑1 (MQ=255)
                                                                                                                                     gATGACCTGGAAAAACAAAGCCCTCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCGGTTGCGCGCGCTGATCCAGCAACAGCTCTTCATGGCCTGTc         <  7:100439/149‑1 (MQ=255)
                                                                                                                                         aCCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGTTGCGCGCGCTGATCCAGCAACAGCTCTTCATGGCCTGTCgcgc     <  8:167545/149‑1 (MQ=255)
                                                                                                                                          ccTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGTTGCGCGCGCTGATc                                  <  4:161987/119‑1 (MQ=255)
                                                                                                                                          ccTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGTTGCGCGCGCTGATc                                  >  3:161987/1‑119 (MQ=255)
                                                                                                                                            tGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGTTGCGCGCGCTGATCCAGCAACAGCTCTTCATGGCCTGTCGCGCTCg  >  5:221552/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GTCGGCGCCCGGTGGCCGAAAGCGTCGTGCAGAGCGGTGTGGACATGGTGGGTATCGGAACGGCGTTAGCCATCGATCCGCACTTGCCTCGTGACTGGTTGCAGGGCAAGGACAACGCCCCTCAATTGCCACCGATCACCTGGAAAAACAAAGCCATCGCATCCTTGGCGAACATGGCAGTAGTGAAGTTTCAACTGCGCAAACTCAGCCGGGCTAAAAAGCCCGACCCGAGCGTTTCGCCGTTGCGCGCGCTGATCCAGCAACAGCTCTTCATGGCCTGTCGCGCTCG  >  NC_002947/3391698‑3391986

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