Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,278,938 A→G 9.1% F353S (TTC→TCC)  PP_1117 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,278,9380AG 9.1% 114.6 / 2.7 44F353S (TTC→TCC) PP_1117hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/21);  new base G (2/2);  total (21/23)
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.30e-01

CATCGGATACGGCCATTGAGACAGCCTGTGGGAAGAATGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCAAGCGTGTCAGCCATTACGACAATCGACTTTCCATGCAGCACGTCGAACAGCATTT  >  NC_002947/1278792‑1279085
                                                                                                                                                   |                                                                                                                                                    
cATCGGATACGGCCATTGAGACAGCCTGTGGGAAGAATGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAg                                                                                                                                                    <  1:73843/149‑1 (MQ=255)
       tACGGCCATTGAGACAGCCTGTGGGAAGAATGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCaa                                                                                                                                             <  7:75152/149‑1 (MQ=255)
             cATTGAGACAGCCTGTGGGAAGAATGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGt                                                                                                                                       <  1:268254/149‑1 (MQ=255)
               ttGAGACAGCCTGTGGGAAGAATGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAg                                                                                                                                     <  8:145877/149‑1 (MQ=255)
                  agaCAGCCTGTGGGAAGAATGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGc                                                                                                                                   <  5:60719/148‑1 (MQ=255)
                     cAGCCTGTGGGTATAATGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACCATGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAg                                                                                                                                 <  1:252182/148‑1 (MQ=255)
                      aGCCTGTGGGAAGAATGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAgg                                                                                                                                >  2:226576/1‑147 (MQ=255)
                        ccTGTGGGAAGAATGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAgg                                                                                                                                >  1:158466/1‑145 (MQ=255)
                          tgtgGGAAGAATGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAgg                                                                                                                                    >  2:171512/1‑139 (MQ=255)
                          tgtgGGAAGAATGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAgg                                                                                                                                    <  1:171512/139‑1 (MQ=255)
                                  gaaTGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTgcgc                                                                                                                  <  2:158466/149‑1 (MQ=255)
                                      gTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTg                                                                                                                     <  4:331922/142‑1 (MQ=255)
                                      gTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTg                                                                                                                     >  3:331922/1‑142 (MQ=255)
                                                            cATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGGAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGgtg                                                                                       >  6:203731/1‑149 (MQ=255)
                                                            cATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGGAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGgtg                                                                                       >  6:198753/1‑149 (MQ=255)
                                                             aTGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGgtg                                                                                       >  6:140337/1‑148 (MQ=255)
                                                                agagaTGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAg                                                                                    <  1:226576/148‑1 (MQ=255)
                                                                      gCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTAtg                                                                                                  >  7:195688/1‑128 (MQ=255)
                                                                      gCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTAtg                                                                                                  <  8:195688/128‑1 (MQ=255)
                                                                         agcagTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGcat                                                                                                           <  4:256841/116‑1 (MQ=255)
                                                                         agcagTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGcat                                                                                                           >  3:256841/1‑116 (MQ=255)
                                                                           cagTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCa                                                                        <  5:140337/149‑1 (MQ=255)
                                                                                     gTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTGTACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTa                                                                                                    >  2:209938/1‑111 (MQ=255)
                                                                                     gTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTGTACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTa                                                                                                    <  1:209938/111‑1 (MQ=255)
                                                                                     gTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGGAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCAcc                                                               <  5:203731/148‑1 (MQ=255)
                                                                                     gTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGGAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCAcc                                                               <  5:198753/148‑1 (MQ=255)
                                                                                     gTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCa                                                              >  7:326942/1‑149 (MQ=255)
                                                                                           agccagcGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCaa                                                                    >  6:9552/1‑137 (MQ=255)
                                                                                           agccagcGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCaa                                                                    <  5:9552/137‑1 (MQ=255)
                                                                                           agccagcGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCaa                                                        <  3:117450/149‑1 (MQ=255)
                                                                                           agccagcGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCaa                                                        >  7:103273/1‑149 (MQ=255)
                                                                                                        cTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCAcc                                                               <  4:68545/129‑1 (MQ=255)
                                                                                                        cTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCAcc                                                               >  3:68545/1‑129 (MQ=255)
                                                                                                          gtgGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCg                                                                                 >  1:239383/1‑109 (MQ=255)
                                                                                                          gtgGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCg                                                                                 <  2:239383/109‑1 (MQ=255)
                                                                                                          gtgGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCAAGCGTGTCAGCCATTa                                         >  3:39033/1‑149 (MQ=255)
                                                                                                           tgGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCAAGCGTGTCAg                                               >  3:53136/1‑142 (MQ=255)
                                                                                                           tgGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCAAGCGTGTCAg                                               <  4:53136/142‑1 (MQ=255)
                                                                                                                                          acGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCCGAGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCAAGCGTGTCAGCCATTACGACAATCGACTTTCCATGCAGCACGTCg            >  1:252516/1‑148 (MQ=255)
                                                                                                                                            gAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCTGTGTTCTCACCAATCACCATTTCAAGCGTGTCAGCCATTACGACAATCGACTTTCCATGCAGCACGTCGAACa        <  8:103273/149‑1 (MQ=255)
                                                                                                                                            gAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCAAGCGTGTCAGCCATTACGACAATCGACTTTCCATGCAGCACGTCGaa          >  4:340127/1‑147 (MQ=255)
                                                                                                                                                  gAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCAAGCGTGTCAGCCATTACGACAATCGACTTTCCa                        >  2:111498/1‑127 (MQ=255)
                                                                                                                                                  gAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCAAGCGTGTCAGCCATTACGACAATCGACTTTCCATGCAGCACGTCGAACAGCAttt  >  4:223987/1‑149 (MQ=255)
                                                                                                                                                  gAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGAATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCAAGCGTGTCAGCCATTACGACAATCGACTTTCCa                        <  1:111498/127‑1 (MQ=255)
                                                                                                                                                   aaGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCAAGCGTGTCAGCCATTACGACAATCGa                               >  8:141699/1‑119 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CATCGGATACGGCCATTGAGACAGCCTGTGGGAAGAATGTAGCGTGCAACACTGCTACACCATGAGAGATGCCAGCAGTTGCCTTGTAAGCAGCCAGCGACGTTCTGTGGAGTTTTGCACC‑TGCTCTTGGGTTATACACGAAGTGGAAGGACTCAACAGCGTAGGCAGGAAGATCATCTGCGCC‑AGCATCATCGTATGTGATGAAGGTGTAGGCGGTGTTCTCACCAATCACCATTTCAAGCGTGTCAGCCATTACGACAATCGACTTTCCATGCAGCACGTCGAACAGCATTT  >  NC_002947/1278792‑1279085

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