Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 774,988 A→T 8.0% N54Y (AAC→TAC)  PP_0666 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947774,9880AT 8.0% 145.5 / 3.5 50N54Y (AAC→TAC) PP_0666hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base T (2/2);  total (28/22)
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.86e-01

CAGGCGCGCGGATGTCACTGGCGGTCAGTTACAGGGGCCTGTTCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGCGGGCCTGCGTGGGTGCCGACCTGGGGGCTGGC  >  NC_002947/774848‑775134
                                                                                                                                            |                                                                                                                                                  
cAGGCGCGCGGATGTCACTGGCGGTCAGTTACAGGGGCCTGTTCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCgg                                                                                                                                             <  2:259078/148‑1 (MQ=255)
                      ggTCAGTTACAGGGGCCTGTTCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGtt                                                                                                                       >  4:162210/1‑148 (MQ=255)
                       gTCAGTTACATGGGCCTGTTCGAAACGGCAGGGATCGTGGGCGAGGAACTGCAACTGGATGTTCAGGTGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCa                                                                                                                     <  8:235703/149‑1 (MQ=255)
                        tCAGTTACAGGGGGCTGTTCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCTACGTCGGCAAGGCGCAGCTGACGCGGGTTc                                                                                                                      >  8:79918/1‑147 (MQ=255)
                        tCAGTTACAGGGGCCTGTTCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTTATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAg                                                                                                                    <  2:346675/149‑1 (MQ=255)
                        tCAGTTACAGGGGCCTGTTCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCTACGTCGGCAAGGCGCAGCTGACGCGGGTTc                                                                                                                      >  8:84777/1‑147 (MQ=255)
                        tCAGTTACAGGGGCCTGTTCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCTACGTCGGCAAGGCGCAGCTGACGCGGGTTc                                                                                                                      <  7:84777/147‑1 (MQ=255)
                        tCAGTTACAGGGGCCTGTTCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCTACGTCGGCAAGGCGCAGCTGACGCGGGTTc                                                                                                                      <  7:79918/147‑1 (MQ=255)
                        tCAGTTACAGGGGCCTGTTCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAg                                                                                                                    <  2:66501/149‑1 (MQ=255)
                            ttACAGGGGCCTGTTCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGgcgc                                                                                                                                     >  5:152823/1‑128 (MQ=255)
                            ttACAGGGGCCTGTTCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGgcgc                                                                                                                                     <  6:152823/128‑1 (MQ=255)
                                         ttCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTa                                                                                                   >  4:294294/1‑149 (MQ=255)
                                            gAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTAcc                                                                                                 <  3:294294/148‑1 (MQ=255)
                                             aaaCGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAg                                                                                                                    <  7:209634/128‑1 (MQ=255)
                                             aaaCGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAg                                                                                                                    >  8:209634/1‑128 (MQ=255)
                                             aaaCGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTgg                                                                                                              >  6:8779/1‑134 (MQ=255)
                                                 ggcaggcaTCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGCCACATCTGAACCCGCGTCAGCTGCGCCTTGCCGa                                                                                                      <  7:64759/138‑1 (MQ=25)
                                                 ggcaggcaTCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGCCACATCTGAACCCGCGTCAGCTGCGCCTTGCCGa                                                                                                      >  8:64759/1‑138 (MQ=25)
                                                    aggcaTCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTc                                                                                         >  1:261424/1‑148 (MQ=255)
                                                                    gACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAgt                                                                        <  5:333969/149‑1 (MQ=255)
                                                                       cTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTa                                                                     >  4:193576/1‑149 (MQ=255)
                                                                             cAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGc                                                               >  8:212060/1‑149 (MQ=255)
                                                                                gATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTg                                                            >  3:145544/1‑149 (MQ=255)
                                                                                  tgtgCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGt                                                          >  2:325150/1‑149 (MQ=255)
                                                                                           gCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCt                                                 >  4:187361/1‑149 (MQ=255)
                                                                                              gTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCtgtg                                              <  3:162210/149‑1 (MQ=255)
                                                                                              gTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCtgtg                                              <  3:193576/149‑1 (MQ=255)
                                                                                                     cAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGtt                                                         <  1:57610/131‑1 (MQ=255)
                                                                                                     cAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGtt                                                         >  2:57610/1‑131 (MQ=255)
                                                                                                       ggCGCTGAGCGTGATCGACGGCCTGATTGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCgg                                     >  1:191347/1‑149 (MQ=255)
                                                                                                        gcgcTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCt                                                 <  3:176586/136‑1 (MQ=255)
                                                                                                        gcgcTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCt                                                 >  4:176586/1‑136 (MQ=255)
                                                                                                              aGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGCggg                               <  7:212060/148‑1 (MQ=255)
                                                                                                                      cGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGCGGGCCTGCGTgg                      >  1:202159/1‑149 (MQ=255)
                                                                                                                       gACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCtgt                                               <  4:207138/123‑1 (MQ=255)
                                                                                                                       gACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCtgt                                               >  3:207138/1‑123 (MQ=255)
                                                                                                                         cGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGCGGGCCTGCGTGGGTg                   <  2:261424/149‑1 (MQ=255)
                                                                                                                            ccTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCtgt                                               <  4:27556/118‑1 (MQ=255)
                                                                                                                            ccTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCtgt                                               >  3:27556/1‑118 (MQ=255)
                                                                                                                            ccTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGCGGGCCTGCGTGGGTGcc                 >  2:96219/1‑148 (MQ=255)
                                                                                                                              tGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGCGGGCCTGCGTGGGTGCCGAc              >  1:186936/1‑149 (MQ=255)
                                                                                                                               gATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGc                                  <  1:126062/128‑1 (MQ=255)
                                                                                                                               gATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGc                                  >  2:126062/1‑128 (MQ=255)
                                                                                                                               gATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGCGGGCCTGCGTGGGTGCCGAcc             >  6:88945/1‑149 (MQ=255)
                                                                                                                                   gTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGCGGGCCTGCGTGGGTGCCGACCTggg         >  6:331400/1‑149 (MQ=255)
                                                                                                                                     gCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCaaaa                                         >  8:144522/1‑115 (MQ=255)
                                                                                                                                     gCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCaaaa                                         <  7:144522/115‑1 (MQ=255)
                                                                                                                                      cAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGCGGGCCTGCGTGGGTGCCGACCTggggg       <  7:21198/148‑1 (MQ=255)
                                                                                                                                       aaGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGCGGGCCTGCGTGGGTGCCGACCTGGGggct     >  4:68792/1‑149 (MQ=255)
                                                                                                                                           cAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGCGGGCCTGCGTGGGTGCCGACCTGGGggctggc  >  1:148812/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
CAGGCGCGCGGATGTCACTGGCGGTCAGTTACAGGGGCCTGTTCGAAACGGCAGGCATCGTGGCCGAGGACCTGCAACAGGATGTGCAGGGGCAGTTGCGTCAGGCGCTGAGCGTGATCGACGGCCTGATGGTGCAAGCCAACGTCGGCAAGGCGCAGCTGACGCGGGTTCAGATGTGGCTTGCCGACTACCGCCACTTCGACCTGGTCAACGAAGTGTACGACGCCTGGTTGCAAGGCTGTGCAAAACCGGTGCGGGCCTGCGTGGGTGCCGACCTGGGGGCTGGC  >  NC_002947/774848‑775134

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