Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,012,826 T→A 7.7% L413Q (CTG→CAG)  prfC → peptide chain release factor 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,012,8260TA 7.7% 149.4 / 2.9 52L413Q (CTG→CAG) prfCpeptide chain release factor 3
Reads supporting (aligned to +/- strand):  ref base T (22/26);  new base A (2/2);  total (24/28)
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.50e-01

ACAACCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGA  >  NC_002947/1012682‑1012970
                                                                                                                                                  |                                                                                                                                                
acaACCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTgcag                                                                                                                                              <  2:351140/149‑1 (MQ=255)
  aaCCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTgcagca                                                                                                                                            <  6:23287/149‑1 (MQ=255)
  aaCCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTgcag                                                                                                                                              >  1:40678/1‑147 (MQ=255)
  aaCCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTgcag                                                                                                                                              >  2:134849/1‑147 (MQ=255)
              aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCgaag                                                                                                                                 >  2:89979/1‑148 (MQ=255)
              aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCAGCTGCTGCAGGccttgac                                                                                                                                >  8:148538/1‑144 (MQ=255)
              aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCAGCTGCTGCAGGccttgac                                                                                                                                <  7:148538/149‑6 (MQ=255)
              aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCAGCTGCTGCAGGc                                                                                                                                      >  7:76439/1‑143 (MQ=255)
              aTCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCAGCTGCTGCAGGc                                                                                                                                      <  8:76439/143‑1 (MQ=255)
               tCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGGGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCgaag                                                                                                                                 >  5:9605/1‑147 (MQ=255)
                ccAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGaagaag                                                                                                                              >  2:197542/1‑149 (MQ=255)
                 cAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGaagaag                                                                                                                              <  5:294071/148‑1 (MQ=255)
                  aGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGt                                                                                                                            <  1:134849/149‑1 (MQ=255)
                   gATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTg                                                                                                                           <  3:357731/149‑1 (MQ=255)
                        gCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACTCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACt                                                                                                                      <  8:161213/149‑1 (MQ=255)
                                   aCCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGttctt                                                                                                           <  1:197542/149‑1 (MQ=255)
                                     cGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTcc                                                                                                         >  3:201746/1‑149 (MQ=255)
                                     cGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTcc                                                                                                         <  6:360997/149‑1 (MQ=255)
                                          gCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGTCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAg                                                                                                    <  8:132627/149‑1 (MQ=255)
                                           cGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAgc                                                                                                   >  7:78334/1‑149 (MQ=255)
                                                   tGGGCTTCACCGGTATCCCGCACTTCGTCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGAAACTGCGTCAAGGCCTGCAGCAGCTGGCCGaagaag                                                                                                                              >  1:16620/1‑115 (MQ=255)
                                                   tGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGaagaag                                                                                                                              <  2:16620/114‑1 (MQ=255)
                                                     ggCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAAc                                                                                           <  7:20116/147‑1 (MQ=255)
                                                         tCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAGAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCACCGAcat                                                                                      >  1:96456/1‑149 (MQ=255)
                                                         tCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGAcat                                                                                      >  6:156187/1‑148 (MQ=255)
                                                           aCCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGAcatcat                                                                                   >  8:176778/1‑149 (MQ=255)
                                                                    ccGCACTTCGCGCCGGGGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTg                                                                           <  6:9605/148‑1 (MQ=255)
                                                                    ccGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGc                                                                          >  8:171549/1‑149 (MQ=255)
                                                                    ccGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGc                                                                          >  5:191276/1‑149 (MQ=255)
                                                                        aCTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCACCGACATCATCCTgg                                                                              <  7:319907/141‑1 (MQ=255)
                                                                        aCTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTgg                                                                              >  8:319907/1‑141 (MQ=255)
                                                                            cgcgCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGcggtcggt                                                                   <  2:40678/148‑1 (MQ=255)
                                                                               cgcCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGtgtg                                                                <  7:171549/149‑1 (MQ=255)
                                                                                 ccGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGtgct                                                              <  4:293833/149‑1 (MQ=255)
                                                                                   ggAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGtgctgc                                                            <  5:156187/149‑1 (MQ=255)
                                                                                            ccgccgCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGAc                                                    <  2:96456/148‑1 (MQ=255)
                                                                                                         cTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTgg                                                                              >  6:378664/1‑109 (MQ=255)
                                                                                                         cTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCg                                      >  8:218717/1‑149 (MQ=255)
                                                                                                            aaaGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCtt                                   <  8:78334/149‑1 (MQ=255)
                                                                                                            aaaGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCt                                    <  4:218615/148‑1 (MQ=255)
                                                                                                                gACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGcggt                                                                       >  4:349455/1‑110 (MQ=255)
                                                                                                                gACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGcggt                                                                       <  3:349455/110‑1 (MQ=255)
                                                                                                                    cGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTgaaggaa                            <  7:218717/149‑1 (MQ=255)
                                                                                                                      cTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAaga                          <  4:201746/149‑1 (MQ=255)
                                                                                                                           aTCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTAc                      >  5:192279/1‑148 (MQ=255)
                                                                                                                            tCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTAc                      >  3:357985/1‑147 (MQ=255)
                                                                                                                              cAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGa                                                     <  7:45752/114‑1 (MQ=255)
                                                                                                                              cAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGa                                                     >  8:45752/1‑114 (MQ=255)
                                                                                                                               aaGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGt                 >  2:153827/1‑149 (MQ=255)
                                                                                                                                      tGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGtt                                                        <  6:368260/103‑1 (MQ=255)
                                                                                                                                      tGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGtt                                                        >  5:368260/1‑103 (MQ=255)
                                                                                                                                              ggCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGa  <  4:357985/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACAACCACGGCACCATCCAGATCGGCGACACCTTCACCGAAGGCGAAGCGCTGGGCTTCACCGGTATCCCGCACTTCGCGCCGGAGCTGTTCCGCCGCGTACGCCTGAAAGACCCGCTGAAATCCAAGCAACTGCGTCAAGGCCTGCAGCAGCTGGCCGAAGAAGGTGCTACTCAGGTGTTCTTCCCTGAGCGCAGCAACGACATCATCCTGGGTGCGGTCGGTGTGCTGCAGTTTGACGTGGTCGCCAGCCGCTTGAAGGAAGAGTACAAGGTGGAGTGTGCCTACGA  >  NC_002947/1012682‑1012970

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: