Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,615,482 C→A 28.5% E177* (GAA→TAA)  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,4820CA28.5% 19.9 / 3.8 31E177* (GAA→TAA) U712_18635Peptide chain release factor 1
Reads supporting (aligned to +/- strand):  ref base C (8/10);  new base A (2/7);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCCTTCTCATGGATGTCGACCTCTACCTCTTCTGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAATACA  >  CP006881/3615353‑3615613
                                                                                                                                 |                                                                                                                                   
tatCCTTCTCATGGATGTCGACCTCTACCTCTTCTGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTGTTGATAttt                                                                                                                               >  3:159248/1‑136 (MQ=255)
  tCCTTCTCATGGATGTCGACCTCTACCTCTTCTGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGt                                                                                                                           <  7:69472/138‑1 (MQ=255)
                 tCGACCTCAACCTCTTCTGCCTCTGGAAGACATGCCAGTGTTGCTGTGGATGTATGGAGGCGTCCGCCTGTTTCTGTTTCCGGAACACGCTGTAGGCGGTGGGCGCCGTTTTAATATTTTAGTTTAGAATAcgcg                                                                                                               <  4:159248/135‑1 (MQ=255)
                           cctctTCTGCCTCCGGAAGACATGCCACTGTTGCTGTGGATGTATGGAGGCGGCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTATAATACGCGCCGCTTCCTGt                                                                                                    <  6:90600/136‑1 (MQ=255)
                           cctctTCTGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGt                                                                                                    <  5:176598/136‑1 (MQ=255)
                                tctGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATAttt                                                                                                                               >  8:49822/1‑104 (MQ=255)
                                tctGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATAttt                                                                                                                               <  7:49822/104‑1 (MQ=255)
                                                  gCCACTGTTGCAGTGGATGTATGGATGCGTACGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAgatgat                                                                                     <  6:392262/128‑1 (MQ=255)
                                                  gCCACTGTTGCAGTGGATGTATGGATGCGTACGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAgatgat                                                                                     >  5:392262/1‑128 (MQ=255)
                                                           gCAGTGGATGTTTGGATGCGTCGGCATGATTATGTTTCCGGAAGAGCCTTTACGGGGGGGGCGCATTTTTAATATTTGAGTTTACAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTc                                                                 <  7:79402/139‑1 (MQ=255)
                                                              gTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTGATATTTGAGTTTTTAATACGCGCCGCTGACTGTAATAATAAAGATGATCTCTTTGTAGCCGCCCGGTCCgg                                                              >  6:238808/1‑139 (MQ=255)
                                                              gTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTGATATTTGAGTTTAGAATACGCGCCCCGTCCTGTAATAATAAAGATGGTGTCTTTGTAGCCGCCCGGTCCgg                                                              >  5:95909/1‑139 (MQ=255)
                                                              gTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTTGAAAACGCGCCCCGTCCTGTAATCATAAAGATTATATATTTGTAGCCGCCCGGTCCgg                                                              >  8:79402/1‑139 (MQ=255)
                                                                   tGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTGTTCATTTTTGAGTTTAGAATACGCCGCGCGTCCTTTTATAATAAATATTTTTAAGTTTTTGCCGCccggtccggtagca                                                         >  5:426734/1‑135 (MQ=255)
                                                                        ggttGCGGCCGCCTTATTCTGTTTCCGGACCAGGGTGTCCGGGGGGGGCGCGGTTTTATTATTTGAGTTAAGAATCCGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACAt                                                        <  3:277431/132‑1 (MQ=255)
                                                                        ggATGCGGCCGCCTGATTCGGTTTCCGCAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTGGAGTTTAGAATCCGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCg                                                     <  1:49179/138‑1 (MQ=255)
                                                                           tGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTAATATTTGAGTTTAGAATACGCGCCGCTTCCTGTTATCATAAAGATGATCGATTTGTAGCCGCCCGTTGCGGTTTCCTTCGcgct                                                 >  1:57840/1‑136 (MQ=255)
                                                                              gTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTTGAATACGCGCCGCGTCCTGTTATCATAAAGATGATCTGTTTTTTGCCGCCCGGTCCGGGTACAATCGGTTGCAt                                              >  5:89815/1‑139 (MQ=255)
                                                                               tCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTGATATTTGAGTTTATAATAAGCGCCGCGGACTGTAATCATAAATATTATCTCTTTGTTTCCCCCCGTTCCGGTTTCATTAGCTTGCATg                                             >  5:15758/1‑139 (MQ=255)
                                                                                 ggcgtGATTGTGTTTCCGCACCCGGCTGTCCGGGGGGGGCGCGTTTTTAATTTTTGTTTTTAGAATCCGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGAc                                           <  1:423190/135‑1 (MQ=255)
                                                                                 cGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATTTTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAgatgat                                                                                     >  6:202083/1‑97 (MQ=255)
                                                                                 cGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAgatgat                                                                                     <  5:202083/97‑1 (MQ=255)
                                                                                           tGTTTCCGGAACACGCTGTCCGCGGTGGGCGCCGTTTTAATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGtt                                    <  3:398161/136‑1 (MQ=255)
                                                                                                   gAACAGGCTGTACGCGGTGGGCGACTTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTg                         <  6:143848/139‑1 (MQ=255)
                                                                                                             tACGCGGTGGGCGCCGTTTTAATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTGTTTAGCCGCCCGGTCCGGGTACATTCGCGTCCATGACCTCTGTTTTGCAGCCCTGAAGCGGAGCa               >  1:141894/1‑139 (MQ=255)
                                                                                                                  ggTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACg          <  6:95909/139‑1 (MQ=255)
                                                                                                                     gggCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACggg        >  6:406863/1‑138 (MQ=255)
                                                                                                                       gcgcTGTTTTAATATTTGAGTTAAGAATACGAGCCGCTTCCTGTAATCATAAAGTTGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAAt     <  6:426734/139‑1 (MQ=255)
                                                                                                                       gcgcGGTTTTAATTTTTTAGTTTATAAACCGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAAt     <  6:15758/139‑1 (MQ=255)
                                                                                                                       gcgcCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGaa      >  5:360651/1‑138 (MQ=255)
                                                                                                                          ccGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAATACa  <  6:89815/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                   
TATCCTTCTCATGGATGTCGACCTCTACCTCTTCTGCCTCAGGAAGACATGCCACTGTTGCAGTGGATGTATGGATGCGTCCGCCTGATTCTGTTTCCGGAACACGCTGTACGCGGTGGGCGCCGTTTTCATATTTGAGTTTAGAATACGCGCCGCTTCCTGTAATCATAAAGATGATCTCTTTGTAGCCGCCCGTTCCGGTTACATTCGCTTCCATGACCTCTGTTTTCCAGCCCTGAAGCTCAGCATAACGGGAATACA  >  CP006881/3615353‑3615613

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: