Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 159,931 A→T 7.3% L153Q (CTG→CAG)  PP_0149 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947159,9310AT 7.3% 156.1 / 2.7 55L153Q (CTG→CAG) PP_0149hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base T (2/2);  total (26/29)
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.55e-01

GCCATTGTGCTGTTCGCGCTCGCCTCGTCGGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGCTCTTCACGCACCGTGCCATGTGGGCGGTCGGTGC  >  NC_002947/159792‑160075
                                                                                                                                           |                                                                                                                                                
gCCATTGTGCTGTTCGCGCTCGCCTCGTCGGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCtt                                                                                                                                         <  4:7974/149‑1 (MQ=255)
  cATTGTGCTGTTCGCGCTCGCCTCGTCGGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGc                                                                                                                                       <  8:271795/149‑1 (MQ=255)
    ttGTGCTGTTCGCGCTCGCCTCGTCGGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCGGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTg                                                                                                                                     <  7:301416/149‑1 (MQ=255)
        gCTGTTCGCGCTCGCCTCGTCGGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGt                                                                                                                                 >  3:244917/1‑149 (MQ=255)
            ttCGCGCTCGCCTCGTCGGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTg                                                                                                                              >  2:41551/1‑148 (MQ=255)
                 gcTCGCCTCGTCGGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGt                                                                                                                        <  6:82772/149‑1 (MQ=255)
                  cTCGCCTCGTCGGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGCGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTg                                                                                                                       <  7:116936/149‑1 (MQ=255)
                   tCGCCTCGTCGGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCGGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTg                                                                                                                       >  1:201285/1‑148 (MQ=255)
                          gtcgGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAgc                                                                                                               <  4:287427/149‑1 (MQ=255)
                           tcgGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAgcg                                                                                                              >  2:237265/1‑149 (MQ=255)
                           tcgGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTCACCAGAGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAgcg                                                                                                              >  4:233684/1‑149 (MQ=255)
                               cccACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCAtt                                                                                                          >  8:183401/1‑149 (MQ=255)
                                   cGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCg                                                                                                      <  2:99181/149‑1 (MQ=255)
                                          ctgcCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTgg                                                                                               <  1:88275/149‑1 (MQ=255)
                                           tggcACCCGTCCAGGTCGAGGAGGCGTGGGTGTTCGCTGTGTGTGAAGCGGCGCTCGGCCTGCACGCGGCATTCTTCCGGCAGGACGACCTTGCCAAGCGTGCCTTTCTGGCGTTTTAAGGGTGCAATCAGCGCATTTTCTGCGCTgg                                                                                               <  4:66418/145‑1 (MQ=255)
                                           tgcCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTggg                                                                                              >  2:158401/1‑149 (MQ=255)
                                                             aGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTc                                                                            <  3:233684/149‑1 (MQ=255)
                                                             aGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTc                                                                            <  4:103913/149‑1 (MQ=255)
                                                                  gCGTGGGTCGTCGCTCGGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTaga                                                                       <  2:201285/149‑1 (MQ=255)
                                                                  gCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTa                                                                         >  5:15074/1‑147 (MQ=255)
                                                                    gTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAgaga                                                                     <  1:41551/149‑1 (MQ=255)
                                                                     tGGGTCGTCGCTCAGGGTGAAGCGGCGCTGGGCCGGCCCCCGGCAGTATTCCGGGCGGCCGACCTTGGCCAGCGTGGCTTTGGGGCGGTTTCCGCGGTGAATAAAGGCATTTTGGGGGGGGGGCCACCCGCGGGGCCGGTCGGagtggt                                                                    >  6:245696/1‑145 (MQ=255)
                                                                         tcgtcgCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGt                                                                <  7:102442/149‑1 (MQ=255)
                                                                               cTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTa                                                                         >  8:167313/1‑134 (MQ=255)
                                                                               cTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTa                                                                         <  7:167313/134‑1 (MQ=255)
                                                                                tCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCgg                                                                                    <  6:35620/122‑1 (MQ=255)
                                                                                tCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCgg                                                                                    >  5:35620/1‑122 (MQ=255)
                                                                                tCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGt                                                                <  4:258670/142‑1 (MQ=255)
                                                                                tCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGt                                                                >  3:258670/1‑142 (MQ=255)
                                                                                  aGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCACAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGAATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTg                                                       >  6:150369/1‑149 (MQ=255)
                                                                                   gggTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTgca                                                     <  4:244917/149‑1 (MQ=255)
                                                                                     gTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTg                                                                                                <  6:18184/105‑1 (MQ=255)
                                                                                     gTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTg                                                                                                >  5:18184/1‑105 (MQ=255)
                                                                                        aaGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCAc                                                 >  3:94108/1‑149 (MQ=255)
                                                                                         aGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTg                                                               >  4:118923/1‑134 (MQ=255)
                                                                                         aGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTg                                                               <  3:118923/134‑1 (MQ=255)
                                                                                              cgcTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGc                                           <  1:158401/149‑1 (MQ=255)
                                                                                              cgcTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGc                                           <  4:94108/149‑1 (MQ=255)
                                                                                               gcTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGc                                           >  2:174387/1‑148 (MQ=255)
                                                                                               gcTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGGGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGCGATGGGTGAGCTCGTGCAGCACGGTGGc                                           >  2:181581/1‑148 (MQ=255)
                                                                                                      cTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGct                                   >  6:302565/1‑149 (MQ=255)
                                                                                                       tGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGctc                                  >  1:128164/1‑149 (MQ=255)
                                                                                                         cccGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGTCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGCTCtt                                <  6:15074/149‑1 (MQ=255)
                                                                                                          ccGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCGGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGCTCTTc                               <  7:183401/149‑1 (MQ=255)
                                                                                                           cgcgGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGCTCTTCa                              >  6:79766/1‑149 (MQ=255)
                                                                                                            gcgGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGCTCTTCAc                             >  5:332443/1‑149 (MQ=255)
                                                                                                                     ttcCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGCTCTTCACGc                           >  7:104302/1‑142 (MQ=255)
                                                                                                                     ttcCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGCTCTTCACGc                           <  8:104302/142‑1 (MQ=255)
                                                                                                                         ggCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGCTCTTCACGCACCGTGCCAtg                <  5:245696/149‑1 (MQ=255)
                                                                                                                            aGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGCTCTTCACGCACCGTGCCATGTgg             >  4:112371/1‑149 (MQ=255)
                                                                                                                                cGACCTTGCCCTGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTgcag                                                    >  6:33164/1‑106 (MQ=255)
                                                                                                                                cGACCTTGCCCTGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTgcag                                                    >  6:27897/1‑106 (MQ=255)
                                                                                                                                cGACCTTGCCCTGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTgcag                                                    <  5:27897/106‑1 (MQ=255)
                                                                                                                                cGACCTTGCCCTGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTgcag                                                    <  5:33164/106‑1 (MQ=255)
                                                                                                                                        cccAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGCTCTTCACGCACCGTGCCATGTGGGCGGTCGGTgc  <  5:79766/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                
GCCATTGTGCTGTTCGCGCTCGCCTCGTCGGCCCACGTACTGCTGCCACCCGGCCAGGTCGAGCAGGCGTGGGTCGTCGCTCAGGGTGAAGCGGCGCTCGGCCTGCCCGCGGCATTCTTCCGGCAGGCCGACCTTGCCCAGCGTGCCTTGCTGGCGTTTGCAGCGTGCAATCAGCGCATTTTCGGCGCTGGGCCACAGGCGGGCACGGTCGTAGAGATGGGTGAGCTCGTGCAGCACGGTGGCCAGCAGCTCTTCACGCACCGTGCCATGTGGGCGGTCGGTGC  >  NC_002947/159792‑160075

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: