Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,033,329 T→C 88.7% L97P (CTG→CCG)  PP_0890 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,033,3290TC88.7% 110.7 / 5.0 44L97P (CTG→CCG) PP_0890hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (20/19);  total (23/21)
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.11e-01

CAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGTGCCTGCCCGCTCATGCCAATTGGCACCTGGATGGCGAAT  >  NC_002947/1033183‑1033465
                                                                                                                                                  |                                                                                                                                        
cAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTtgctgc                                                                                                                                        <  8:17007/149‑1 (MQ=255)
   ggTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCAt                                                                                                                                     >  5:224651/1‑149 (MQ=255)
       tACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGaa                                                                                                                                  <  5:88115/148‑1 (MQ=255)
         cTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATg                                                                                                                                >  8:89349/1‑148 (MQ=255)
            ccAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAgg                                                                                                                            >  8:27376/1‑149 (MQ=255)
                         aaGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTgg                                                                                                               >  5:272367/1‑149 (MQ=255)
                                     ccACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGttt                                                                                                    >  5:255317/1‑148 (MQ=255)
                                     ccACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGttt                                                                                                    >  5:250876/1‑148 (MQ=255)
                                     ccACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGttt                                                                                                    >  3:90776/1‑148 (MQ=255)
                                       aCTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCAc                                                                                                 >  2:324420/1‑149 (MQ=255)
                                       aCTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCAc                                                                                                 >  4:244111/1‑149 (MQ=255)
                                           aGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACga                                                                                               <  2:254133/147‑1 (MQ=255)
                                              ctgctgGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCt                                                                                          >  2:222341/1‑149 (MQ=255)
                                                 ctgGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGt                                                                                        <  5:234360/148‑1 (MQ=255)
                                                    gaggaAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATg                                                                                    <  4:90776/149‑1 (MQ=255)
                                                       gaAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAAGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGtt                                                                                  >  5:120729/1‑148 (MQ=255)
                                                            cagccagcAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGtt                                                                            <  4:81632/149‑1 (MQ=255)
                                                               ccagcAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGcc                                                                         <  5:253762/149‑1 (MQ=255)
                                                                   cagcagcGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGTCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGt                                                                     <  6:255317/149‑1 (MQ=255)
                                                                   cagcagcGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGt                                                                     <  6:250876/149‑1 (MQ=255)
                                                                       agcGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGt                                                                                                      >  1:157227/1‑112 (MQ=255)
                                                                       agcGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGt                                                                                                      <  2:157227/112‑1 (MQ=255)
                                                                               aCCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGcc                                                         >  1:33791/1‑149 (MQ=255)
                                                                               aCCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGcc                                                         <  3:244111/149‑1 (MQ=255)
                                                                               aCCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGcc                                                         <  1:324420/149‑1 (MQ=255)
                                                                               aCCGATGAAGGCATGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGATGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGcc                                                         >  6:40004/1‑149 (MQ=255)
                                                                                ccGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCCGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCt                                                        >  6:49270/1‑149 (MQ=255)
                                                                                      aaGGCCTGGGGCGGCGCAGTCAAGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTg                                                   <  6:120729/148‑1 (MQ=255)
                                                                                       aGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCt                                                 <  7:27376/149‑1 (MQ=255)
                                                                                        ggCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTg                                                <  6:272367/149‑1 (MQ=255)
                                                                                        ggCATGGGGCGGCGCAGTCAGGCGGCCTGGCAATGCGCGCGTGAGGCCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCATGGTTTAACGCTCTCTGTCATGTTCAAGTTGCCACGTCTGCTTTCAGCCTTATTGCTg                                                >  6:71017/1‑149 (MQ=255)
                                                                                               ggcggcGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCtgt                                         <  6:224651/149‑1 (MQ=255)
                                                                                                     gcAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGtgcctgc                                   <  2:33791/149‑1 (MQ=255)
                                                                                                         tCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGTGCCTGCCCGc                               >  6:271719/1‑149 (MQ=255)
                                                                                                         tCAGGCGGCCCGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGTGCCTGCCCGc                               <  5:49270/149‑1 (MQ=255)
                                                                                                           aGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGTGCCTGCCCGCt                              >  2:211679/1‑148 (MQ=255)
                                                                                                             gcggcCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGc                                                          <  4:248731/118‑1 (MQ=255)
                                                                                                             gcggcCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGc                                                          >  3:248731/1‑118 (MQ=255)
                                                                                                             gcggcCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGTGCCTGCCCGCTCAt                           >  5:26554/1‑149 (MQ=255)
                                                                                                             gcggcCCGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGTGCCTGCCCGCTCAt                           >  6:111059/1‑149 (MQ=255)
                                                                                                                       aaGGCGCGCGTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGTGCCTGCCCGCTCATGCCAATTgg                  >  3:304976/1‑148 (MQ=255)
                                                                                                                           cgcgcgTGAGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGTGCCTGCCCTCTCATGCCAATTGGCACCt             <  5:40004/149‑1 (MQ=255)
                                                                                                                                   aGGGCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGTGCCTGCCCGCTCATGCCAATTGGCACCTGGATGGc      <  7:89349/148‑1 (MQ=255)
                                                                                                                                      gCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGTGCCTGCCCGCTCATGCCAATTGGCACCTGGATGGCGAAt  >  6:310725/1‑149 (MQ=255)
                                                                                                                                      gCTGCCGCTTGCCGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGTGCCTGCCCGCTCATGCCAATTGGCACCTGGATGGCGAAt  <  1:222341/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                        
CAAGGTGTACTGCCAAGGGTACGTGAAGCAGTGGCTGCCACTGAGCCTGCTGGAGGAAAACAGCCAGCAGCAGCGAGTCACCGCTGAAGGCCTGGGGCGGCGCAGTCAGGCGGCCTGGCAAGGCGCGCGTGAGGGCTGCCGCTTGCTGCCATGAATGCAGGCGGGTAAGGTTGGCGGCCTGGTTTCACGAGATCTGTCATGTTCAAGTTGCCCCGTCTGCTGTCTGCCTTATTGCTGGCTCTGTGCCTGCCCGCTCATGCCAATTGGCACCTGGATGGCGAAT  >  NC_002947/1033183‑1033465

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