Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,139,001 T→A 100% I275F (ATC→TTC)  rne ← fused ribonucleaseE: endoribonuclease/RNA‑binding protein/RNA degradosome binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,139,0010TA82.9% 61.7 / 5.9 29I275F (ATC→TTC) rnefused ribonucleaseE: endoribonuclease/RNA‑binding protein/RNA degradosome binding protein
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base A (15/9);  total (19/10)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Not supported by required number of reads on both strands.

TTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGATCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCG  >  CP009273/1138862‑1139143
                                                                                                                                           |                                                                                                                                              
ttAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACtg                                                                                                                                             >  2:134615/1‑143 (MQ=255)
               tcGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGAGCGGGTa                                                                                                                              >  2:4380/1‑143 (MQ=255)
               tcGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCATCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGAGCTGGTa                                                                                                                              >  2:213196/1‑143 (MQ=255)
                        ccgccgcGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGATCTGTGACTCGATCTGGTAGTGGATGAACAgggg                                                                                                               >  1:158232/1‑149 (MQ=255)
                        ccgccgcGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGATCTGTGACTCGATCTGGTAGTCGCTGAACAGCgg                                                                                                               >  1:158226/1‑149 (MQ=255)
                                   cgcgcgTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCgg                                                                                                    >  1:7749/1‑149 (MQ=255)
                                                tggATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGAt                                                                                       <  1:251072/147‑1 (MQ=255)
                                                  gATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGAttt                                                                                     >  1:133461/1‑149 (MQ=255)
                                                    tGTCGATGGCCGTTACCGGCTCGGTGCTGTCGATAACAACGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGCCTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTtg                                                                                   <  1:134877/149‑1 (MQ=255)
                                                    tGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGGCTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTtg                                                                                   <  1:325281/149‑1 (MQ=255)
                                                                ttacagtaaccagtcTGGCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGGTTGATTTTGctgctg                                                                             >  2:86605/15‑143 (MQ=255)
                                                                     gCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCg                                                                  >  1:62916/1‑149 (MQ=255)
                                                                     gCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGGGTACAGTTTGATTTTGCTGCTGAAATCCGGGCg                                                                  >  1:292071/1‑149 (MQ=255)
                                                                       ttCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCg                                                                      >  2:89509/1‑143 (MQ=255)
                                                                         cGGTGCTGGCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGATCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCggg                                                                    <  2:15572/143‑1 (MQ=255)
                                                                                           tGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGa                                            <  1:152312/149‑1 (MQ=255)
                                                                                            ggAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGATCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTggc                                           >  1:108510/1‑147 (MQ=255)
                                                                                            ggAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGAc                                           <  1:149569/149‑1 (MQ=255)
                                                                                                    caGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGt                                         >  2:323992/1‑143 (MQ=255)
                                                                                                                gAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTc                       <  1:5063/149‑1 (MQ=255)
                                                                                                                gAACTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATATTGCTGCTGAAATCCGGGCGACCTACTGCAGCGCTATGCTGACGTGCCAGTTCGAGCACTTTc                       <  1:239504/149‑1 (MQ=255)
                                                                                                                 aaCTTCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCa                            <  2:19856/143‑1 (MQ=255)
                                                                                                                    ttCACGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGGGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGGGCACTTTCGGGt                   >  1:45034/1‑149 (MQ=255)
                                                                                                                       cggcgggggggCGGGCTCGGACTGTGGCTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCg                      <  2:119809/134‑1 (MQ=255)
                                                                                                                        cGCTGGAAGGCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGCGTTATc               >  1:83711/1‑149 (MQ=255)
                                                                                                                                 gCGGACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTGATCGAt            >  2:46269/1‑143 (MQ=255)
                                                                                                                                   ggACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGGGCACTTTCGGGCTATCGCTAAGGCttt    >  1:103146/1‑149 (MQ=255)
                                                                                                                                   ggACTCGAACTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGAttt    >  1:355866/1‑149 (MQ=255)
                                                                                                                                     aCTCGATCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTGTCGGGTTATCGATAAGGATTTCg  >  1:116546/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                              
TTAAACGCGGTTTCTTCGATATCGCCGCCGCGGGTCGCGCGTGCGGAGTTGATGTCGATGGCCGTTAACGCTTCGGTGCTGTCGATAACAATGGAACCACCAGACGGCAGACGAACTTCACGCTGGAAGGCGGACTCGATCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCG  >  CP009273/1138862‑1139143

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: