Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,508,789 A→T 17.3% M109K (ATG→AAG) ‡ PP_2203 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,508,7890AT17.3% 82.4 / 10.1 45M109K (ATG→AAG) ‡PP_2203hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/14);  new base T (3/5);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

GGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGC  >  NC_002947/2508642‑2508932
                                                                                                                                                   |                                                                                                                                               
ggcggcACGTAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGGATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCAt                                                                                                                                                <  2:221178/149‑1 (MQ=255)
       ccAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCGCCCGGc                                                                                                                                         <  1:309444/149‑1 (MQ=255)
        cAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCa                                                                                                                                        >  5:34671/1‑149 (MQ=255)
         aaaGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCCCCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCACGGCAc                                                                                                                                       >  1:1819/1‑149 (MQ=255)
               gcTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCGCCCGGCat                                                                                                                                       >  6:174921/1‑142 (MQ=255)
               gcTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCGCCCGGCat                                                                                                                                       <  5:174921/143‑2 (MQ=255)
                   ataaaCCGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGCCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGcc                                                                                                                              <  3:181167/146‑1 (MQ=255)
                   agaAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGccc                                                                                                                             >  3:1161/1‑149 (MQ=255)
                         aGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAgg                                                                                                                       >  3:272080/1‑149 (MQ=255)
                         aGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGAGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGa                                                                                                                         <  3:71276/147‑1 (MQ=255)
                               ccTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCGCCCGGCat                                                                                                                                       >  1:40842/1‑126 (MQ=39)
                               ccTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCGCCCGGCat                                                                                                                                       <  2:40842/127‑2 (MQ=39)
                                          cATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCcag                                                                                                      >  3:296544/1‑149 (MQ=255)
                                                           cACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCGGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTcgg                                                                                     <  4:86594/149‑1 (MQ=255)
                                                               gtgtCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTcggccgg                                                                                 >  5:308363/1‑149 (MQ=255)
                                                                        ccATCCGGGGCCACCCCGATCCCCTGGTCGGCATCTTGGGGGTCCCCCATCTCCATGCCGGGCCACCCCTTGGGATTCCCCGCCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGcc                                                                        <  4:206501/149‑1 (MQ=255)
                                                                           tCCGGGGCCCCCCCGATGACCTGGTAGCCATCCTTGGGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATCCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTc                                                                     <  1:284702/149‑1 (MQ=255)
                                                                                 gTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTcgacga                                                                >  3:231392/1‑148 (MQ=255)
                                                                                 ctCACCCCGATGCCGTGTTAGGCTTCCTTGGTGCCACCCATCTCCTTGCGGGGCTCGCCCTTGGGGGTGCCGGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGaa                                                               <  4:185252/148‑1 (MQ=255)
                                                                                        cGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGCCGATCTTgccg                                                        >  6:75722/1‑149 (MQ=255)
                                                                                            cagcTGGTAGGCAGCCTTGTGGCCACCCATCTCAAGCCGGGGGGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGa                                                    <  5:105906/146‑1 (MQ=255)
                                                                                             aCCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGAt                                                   >  6:158664/1‑149 (MQ=255)
                                                                                              ccTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATc                                                  >  7:56107/1‑149 (MQ=255)
                                                                                               cTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCa                                                 >  7:177021/1‑149 (MQ=255)
                                                                                                       cATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCACATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGACAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTAGACCAACTTTCCGCCGATCACCCCGgtg                                         >  7:268871/1‑149 (MQ=255)
                                                                                                           cTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTgg                                      <  5:308170/148‑1 (MQ=255)
                                                                                                              gtgGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCt                                                                       <  1:77801/112‑1 (MQ=255)
                                                                                                              gtgGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCt                                                                       >  2:77801/1‑112 (MQ=255)
                                                                                                                 gTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCacgacg                               >  3:262295/1‑149 (MQ=255)
                                                                                                                 gTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCacgacg                               >  5:244259/1‑149 (MQ=255)
                                                                                                                  tCACCCATCTCCATGCCGGGCGAGCCCTTGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGc                              <  2:8463/149‑1 (MQ=255)
                                                                                                                      ccATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAgg                          >  6:101364/1‑149 (MQ=255)
                                                                                                                             cttGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGt                       <  4:231392/143‑1 (MQ=255)
                                                                                                                             cATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGgcg                   >  7:156427/1‑149 (MQ=255)
                                                                                                                               tGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGgcgc                  >  2:304701/1‑148 (MQ=255)
                                                                                                                                   gggCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGc             >  2:118144/1‑149 (MQ=255)
                                                                                                                                   ggcgtgccccAGGGCCTTGCCGGCCACCGCCAGGCCCTTTAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGTCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGc             <  8:268871/142‑1 (MQ=255)
                                                                                                                                      cGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGa          <  4:296544/149‑1 (MQ=255)
                                                                                                                                       gAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAg         <  8:177021/149‑1 (MQ=255)
                                                                                                                                        aGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCTCgaccaat        >  1:76295/1‑143 (MQ=255)
                                                                                                                                          cccATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCg      >  5:231774/1‑149 (MQ=255)
                                                                                                                                           ccATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGGGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGt     >  5:265574/1‑149 (MQ=255)
                                                                                                                                           ccATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGt     >  2:215153/1‑149 (MQ=255)
                                                                                                                                           ccATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGt     <  1:215153/149‑1 (MQ=255)
                                                                                                                                              tGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGc  <  5:158664/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGC  >  NC_002947/2508642‑2508932

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑