Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,615,489 G→T 23.2% L174L (CTC→CTA U712_18635 ← Peptide chain release factor 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,615,4890GT23.2% 43.3 / 3.0 36L174L (CTC→CTAU712_18635Peptide chain release factor 1
Reads supporting (aligned to +/- strand):  ref base G (12/15);  new base T (6/3);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGGATGTCGACCTCTACCTCTTCTGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAATACATTC  >  CP006881/3615362‑3615616
                                                                                                                               |                                                                                                                               
cATGGATGTCGACCTCTACCTCTTCTGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAAt                                                                                                                       >  4:188092/1‑138 (MQ=255)
           acctctacctctTCTGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTGTTGATATTTTAGTTTACAATACGCGCCGCGTc                                                                                                           >  1:52891/1‑139 (MQ=255)
               ctacctgtTCTGCCTCTGTAAGACATCCCAGTGTGGCTGGGTATGTATGGATGCGTCCGCCTTATTCTGTTTCCGGAACCCGCTGTACGGGGTGGGCGCGGTTTTAATATTTGTGTTTAGAATACGCGCCGCTTCCTGt                                                                                                       <  6:338396/139‑1 (MQ=255)
               ctacctctTCTGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTg                                                                                                        >  8:64003/1‑138 (MQ=255)
                  cctctTCTGCCTCAGGAAGACATGCCACTGTTGCAGTGGATTTATGGATGCGGCCGCCTGATGATGTTTCCGGAACACGCGGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGt                                                                                                       <  8:7651/136‑1 (MQ=255)
                        ctGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTTAGTTTAGAATACGCGCCGGTTCCTGTAATCATaaa                                                                                              >  8:362288/1‑139 (MQ=255)
                        ctGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTGTTCATATTTTAGTTTAGAATACCCGCCGCGTCCGGTAATCATaa                                                                                               >  5:338396/1‑138 (MQ=255)
                              ccgcaaGACATGCCACTGTTGCTGTGGTTGTTTGGAGGGGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGGGGTGGGCGCGGTTTTAATATTTGAGTTTAGAATCCGCGCCGCTTCCTGTAATCATAAAgatgat                                                                                        <  6:158842/135‑1 (MQ=255)
                               aGGACGACATGCAAGTGTTGCTGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACAGGCTGTACGCGGTGGGGGCCGTTTTCATATTTGAGTTTAGATTCCGCGCCGCTTCCTGTAATCATAAAgatgat                                                                                        <  7:362288/138‑1 (MQ=255)
                                     aCATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTg                                                                                 >  1:307513/1‑139 (MQ=255)
                                         gCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCAGTTATCATAAAGATGATCTCTTTGTTgc                                                                             >  8:50826/1‑139 (MQ=255)
                                                              tGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCg                                                        <  7:64003/139‑1 (MQ=255)
                                                               ggATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTTTACGCGGTGGGCGCGGTTTTAATATTTGAGTTTAGAAAACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCg                                                        <  3:143742/138‑1 (MQ=255)
                                                                gATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTAATATTTTAGTTTAGAAAAAGCCCCGCTTCCTGTAATCATAAAGAAGAAATCTTTGTAGCCGCCCGGTCCGGGTAAATTCGCt                                                      >  5:274231/1‑139 (MQ=255)
                                                                  tGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTATAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTATCCGCCCCTTCCGGTTACATTCGCTTc                                                    >  3:153887/1‑139 (MQ=255)
                                                                  tGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTGTTGATATTTTTGTTTAGAATACGCGCCGCGTACGGTAATAATAAAGATGATCTGTTTGTTGCCGCCCGTTCCGGGTCAAATCGCTTc                                                    >  4:105146/1‑139 (MQ=255)
                                                                        cGCCTGATTCTGTTTCGGGAACACGCTGTACGCGGTGGGCGCGGTTTTAATATTTGAGTTTATAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGAc                                              <  2:307513/139‑1 (MQ=255)
                                                                          ccTGATTGTGTTTCGGGAACACGCTGTACGGGTTGGGCGGGGTTTTAATTTTTTAGTTTATATTACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACct                                            <  6:274231/139‑1 (MQ=255)
                                                                                   ttttCGGTAACACGCTGTACGGGGTGGTGGCCGTTTTAATTTTTTAGTTTAGAATCCGCGACGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTc                                    <  2:52891/137‑1 (MQ=255)
                                                                                    tttCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTTAGTTTAGAATACGCGCCGCTTCCTGGAATCATAAAGATGATGTCTTTTTTGCCGCCCGTTCCGGTTTCATTCGCGTCCCTGACCGCTGTTTTc                                    >  5:254649/1‑137 (MQ=255)
                                                                                      tCCGGAACAGGGTGTCCGGGGTGGGCGAGTTTTTAATATTTGTTTTTATATTACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAg                                 <  6:254649/138‑1 (MQ=255)
                                                                                             caggCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAg                         <  7:50826/136‑1 (MQ=255)
                                                                                               acaggTACGCGGGGGGCGCCGTTTTCATATTTGCGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCt                       <  7:24295/135‑1 (MQ=255)
                                                                                                    tACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTTATCATAAAGATGATCACTTTGTAGCCGCCCCTTCCGCTTAAATTAGCTTCCAAGACCTCTGTTTTCCAGCCCTGAAGCTCAGc                   >  4:168450/1‑138 (MQ=255)
                                                                                                          gTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCGTCCGGTAATCATAAAGATGATCTCTTTGTTGCCGCCCGGTCCGGGTACATTCGCGTCCATGACCTCTGTTTTCCCGCCCTGAAGCTCAGCCTAACgg            >  3:21004/1‑139 (MQ=255)
                                                                                                            gggCGCCGTTTTCATATTTGAGTTTATAATACGAGCCGCTTCCTGTAATAATAAATATGATCTAGTTGTAGCCTCCCGTTCCGGTTACATTCGCTTACATGACATATTTTTTCCAGCCCTGAAGCTCATCATAActgt           >  6:176578/1‑135 (MQ=255)
                                                                                                            gggCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACggg           >  4:106570/1‑138 (MQ=255)
                                                                                                              gcgcGGTTTTAATATTTTTTTTTTTAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAAt        <  4:153887/139‑1 (MQ=255)
                                                                                                              gcgcCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATTATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCAt                 >  3:137862/1‑130 (MQ=255)
                                                                                                              gcgcCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCAt                 <  4:137862/130‑1 (MQ=255)
                                                                                                              gcgcCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAAt        <  3:188092/139‑1 (MQ=255)
                                                                                                              gcgcCGTTTTCATATTTGAGTTAATAATCCGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAAt        <  1:245024/139‑1 (MQ=255)
                                                                                                              gcgcCGTTTTCATATTTGAGTTAAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAAt        <  3:105146/139‑1 (MQ=255)
                                                                                                              gcgcCGTTTTAAAATTTGAGTTTTGAATAACCGCCCCGTCCGGTTATCATAAAAATGATCTCTTTGTAGCCGCCCCGTCCGGTTTCATTCGCTTGCATTACCGCTGGTTTGCCGCCCGGAAGCTCACCAAAACGGGagg        >  8:24295/1‑137 (MQ=255)
                                                                                                                 cggTTTTAATTTTTGCTTTTAGAATCCGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAATACa     <  4:21004/137‑1 (MQ=255)
                                                                                                                    ttttCATATTTGAGTTTAGAATACGCGCCGCTTACTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAATACATTc  >  3:156400/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                               
CATGGATGTCGACCTCTACCTCTTCTGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAATACATTC  >  CP006881/3615362‑3615616

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