Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,572,293 C→T 10.7% E175K (GAG→AAG)  orn ← oligoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,572,2930CT10.7% 96.8 / 3.8 39E175K (GAG→AAG) ornoligoribonuclease
Reads supporting (aligned to +/- strand):  ref base C (19/15);  new base T (2/2);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

TCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGTTCCTTGAGGGTGGAGACATCCAGGTTGCGGTAGT  >  NC_002947/5572144‑5572436
                                                                                                                                                     |                                                                                                                                               
tCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTc                                                                                                                                                 >  1:75182/1‑150 (MQ=255)
               aTAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAg                                                                                                                                  >  1:463657/1‑150 (MQ=255)
                              aacaacATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGAtta                                                                                                                   >  1:1038741/1‑149 (MQ=255)
                               acaacaTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTNACCAGTNGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTcg                                                                                                                  >  1:978249/1‑150 (MQ=255)
                                  acaTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCgg                                                                                                               <  2:1058551/150‑1 (MQ=255)
                                          cGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCGGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTGGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTc                                                                                                       <  2:727282/150‑1 (MQ=255)
                                           gggTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTcc                                                                                                      >  2:28116/1‑150 (MQ=255)
                                                   gCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCTAg                                                                                              >  1:649867/1‑150 (MQ=255)
                                                      ggAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTg                                                                                           >  2:900707/1‑150 (MQ=255)
                                                       gAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTgg                                                                                          <  1:61337/150‑1 (MQ=255)
                                                            gcgcgcAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTTGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGa                                                                                                              <  2:892510/125‑1 (MQ=255)
                                                            gcgcgcAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTTGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGa                                                                                                              <  2:891170/125‑1 (MQ=255)
                                                            gcgcgcAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTTGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGa                                                                                                              >  1:892510/1‑125 (MQ=255)
                                                            gcgcgcAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTTGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGa                                                                                                              >  1:891170/1‑125 (MQ=255)
                                                                    tCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATCGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttct                                                                             <  2:539290/150‑1 (MQ=255)
                                                                    tCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttct                                                                             >  2:767890/1‑150 (MQ=255)
                                                                     cTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttctt                                                                            <  2:75182/150‑1 (MQ=255)
                                                                       aaCCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTACTGCCCTTCTTGa                                                                          >  1:601630/1‑150 (MQ=255)
                                                                        aCCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTAGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGaa                                                                         <  2:1038741/150‑1 (MQ=255)
                                                                          cAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGc                                                                       >  1:645251/1‑150 (MQ=255)
                                                                                  cAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTcacgc                                                               <  1:28116/150‑1 (MQ=255)
                                                                                     ggCNCCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTcacgcacg                                                            >  2:699080/1‑150 (MQ=255)
                                                                                        aCCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTcc                                                         >  2:864791/1‑150 (MQ=255)
                                                                                         ccAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCgg                                                                                                                             >  1:404828/1‑81 (MQ=255)
                                                                                         ccAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCgg                                                                                                                             <  2:404828/81‑1 (MQ=255)
                                                                                          cAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCgg                                                       >  1:1075425/1‑150 (MQ=255)
                                                                                          cAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCgg                                                       >  1:147574/1‑150 (MQ=255)
                                                                                                 cggCGCAGATATTGTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCAAGTCCGGCGCCCAg                                                >  2:1084781/2‑150 (MQ=255)
                                                                                                      cAGATACAATATCTGCGCCGTTTCGACTCACAGCTCGCTGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAgcgcgc                                           <  1:900707/150‑1 (MQ=255)
                                                                                                            cAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGCGCGCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGAGATCGGCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAg                                     <  2:242098/150‑1 (MQ=255)
                                                                                                                atCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGTTcc                                 <  2:1075425/150‑1 (MQ=255)
                                                                                                                     cgcCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGTTCCTTGAg                            <  2:978249/150‑1 (MQ=255)
                                                                                                                      gcCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGTTCCTTGAgg                           <  2:147574/150‑1 (MQ=255)
                                                                                                                                 tCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGTTCCTTGAGGGTGGAGACATc                >  1:439701/1‑150 (MQ=255)
                                                                                                                                   acaGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCCTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCTCCCAGCGCGCGGCCAGTTCCTTGAGGGTGGAGACATCCa              <  2:463657/150‑1 (MQ=255)
                                                                                                                                          gatgaAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGTTCCTTGAGGGTGGAGACATCCAGGTTGCg       >  2:1067099/1‑150 (MQ=255)
                                                                                                                                            tgaAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGTTCCTTGAGGGTGGAGACATCCAGGTTGCGGt     >  2:860018/1‑150 (MQ=255)
                                                                                                                                               aGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGGCACGCACGTCCGGCGCCCAGCGCGCGGCCAGTTCCTTGAGGGTGGAGACATCCAGGTTGCGGTAGt  <  2:645251/150‑1 (MQ=255)
                                                                                                                                                    tCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGTTCCTTGAGGGTGga                      >  1:954763/1‑125 (MQ=255)
                                                                                                                                                    tCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGTTCCTTGAGGGTGga                      <  2:954763/125‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGGCCAGTTCCTTGAGGGTGGAGACATCCAGGTTGCGGTAGT  >  NC_002947/5572144‑5572436

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