Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,342,483 G→T 11.4% A229A (GCC→GCA)  rnb ← ribonuclease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,342,4830GT11.4% 132.7 / 10.4 53A229A (GCC→GCA) rnbribonuclease II
Reads supporting (aligned to +/- strand):  ref base G (19/28);  new base T (3/3);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CGCGCAGTGAGCAGAGATCGTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATGGTGACAAAATCCAGCGCGGTCAGATCTTCACGAACCAGACCTTCA  >  CP009273/1342347‑1342618
                                                                                                                                        |                                                                                                                                       
cgcgcAGTGAGCAGAGATCGTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCCGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTGTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCa                                                                                                                                      >  1:577922/1‑140 (MQ=255)
cgcgcAGTGAGCAGAGATCGTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCa                                                                                                                                      >  1:599127/1‑140 (MQ=255)
                 tCGTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGATGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAAGCGCCACAATCAGCTg                                                                                                                     >  1:312969/1‑140 (MQ=255)
                 tCGTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTg                                                                                                                     >  1:130298/1‑140 (MQ=255)
                 tCGTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTg                                                                                                                     >  1:111612/1‑140 (MQ=255)
                 tCGTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTg                                                                                                                     >  1:617144/1‑140 (MQ=255)
                 tCGTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTg                                                                                                                     >  1:234732/1‑140 (MQ=255)
                 tCGTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTATGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTg                                                                                                                     >  1:677657/1‑140 (MQ=255)
                 tCGTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCAGGCAGAGAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTg                                                                                                                     >  1:72161/1‑140 (MQ=255)
                   gTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGac                                                                                                                   >  1:81556/1‑139 (MQ=255)
                   gTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGaa                                                                                                                   >  1:507709/1‑140 (MQ=255)
                   gTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGaa                                                                                                                   >  1:430643/1‑140 (MQ=255)
                            tttgcgAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCAt                                                                                                          <  1:633473/137‑1 (MQ=255)
                            tttgcgAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCAt                                                                                                          <  1:376169/137‑1 (MQ=255)
                                                tGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGaaa                                                                                      <  2:139987/140‑1 (MQ=255)
                                                tGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGaaa                                                                                      <  2:641686/140‑1 (MQ=255)
                                                tGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGaaa                                                                                      <  2:407369/140‑1 (MQ=255)
                                                  ttGAATCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCTGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAg                                                                                    >  1:162076/1‑140 (MQ=255)
                                                  ttGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCTGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAg                                                                                    >  1:603634/1‑140 (MQ=255)
                                                  ttGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCTGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAg                                                                                    >  1:162123/1‑140 (MQ=255)
                                                  ttGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGcc                                                                                               >  1:512338/1‑129 (MQ=255)
                                                  ttGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGcc                                                                                               <  2:512338/129‑1 (MQ=255)
                                                     aaGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGc                                                                                 <  2:546745/140‑1 (MQ=255)
                                                     aaGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGc                                                                                 <  2:432112/140‑1 (MQ=255)
                                                     aaGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGc                                                                                 <  2:428092/140‑1 (MQ=255)
                                                     aaGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGc                                                                                 <  2:31944/140‑1 (MQ=255)
                                                       gCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGt                                                                               <  2:205445/140‑1 (MQ=255)
                                                       gCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGt                                                                               <  2:181056/140‑1 (MQ=255)
                                                           ggCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATc                                                                           <  2:165231/140‑1 (MQ=255)
                                                           ggCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATc                                                                           >  1:688619/1‑140 (MQ=255)
                                                            gCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATcc                                                                          >  2:632365/1‑140 (MQ=255)
                                                             cAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCTGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCa                                                                         <  2:162123/140‑1 (MQ=255)
                                                             cAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCTGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCa                                                                         <  2:162076/140‑1 (MQ=255)
                                                             cAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCTGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCa                                                                         <  2:603634/140‑1 (MQ=255)
                                                                    ttGGTGATTGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGACTTAAGCAATCCTCGCGGTTGGATCGGCAATAGCCACAATCAGCTGAAGTTTGTCATCCGGCAAGGCCTTAGCGAAAAGGGCGTCATCCATAtctta                                                                  <  1:181553/140‑2 (MQ=255)
                                                                    ttGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATAtcttc                                                                  <  1:297011/140‑1 (MQ=255)
                                                                    ttGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATAtcttc                                                                  <  1:162298/140‑1 (MQ=255)
                                                                    ttGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATAtcttc                                                                  <  1:181858/140‑1 (MQ=255)
                                                                    ttGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATAtcttc                                                                  <  1:584753/140‑1 (MQ=255)
                                                                    ttGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATAtcttc                                                                  <  1:208769/140‑1 (MQ=255)
                                                                    ttGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATAtcttc                                                                  <  1:406281/140‑1 (MQ=255)
                                                                    gtGGTGAATGCGCGAATTTTCGCGGCTTTGTCCATCTTGCTGCATTCAGCAATCCACGCGGTTGGAGCGGCAATAGCCACAATCAGCTGAAGTTTGGCTTCCGGCAACGCCTTAGCGAAACGGGCGTCATCCATAtctac                                                                  <  1:288720/139‑3 (MQ=255)
                                                                                        cgcgGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATgg                                              <  2:192467/140‑1 (MQ=255)
                                                                                        cgcgGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATgg                                              <  2:30/140‑1 (MQ=255)
                                                                                        cgcgGATTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATgg                                              <  2:265173/140‑1 (MQ=255)
                                                                                             cTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATGGTGACa                                         <  2:617144/140‑1 (MQ=255)
                                                                                             cTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATGGTGACa                                         <  2:130298/140‑1 (MQ=255)
                                                                                             cTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATGGTGACa                                         <  2:234732/140‑1 (MQ=255)
                                                                                                                        tCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATGGTGACAAAATCCAGCGCGGTCAGATCTTCACGa              >  1:209731/1‑140 (MQ=255)
                                                                                                                                    gATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATGGTGACAAAATCCAGCGCGGTCAGATCTTCACGAACCAGACCTTCa  >  1:314632/1‑140 (MQ=255)
                                                                                                                                    gATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATGGTGACAAAATCCAGCGCGGTCAGATCTTCACGAACCAGACCTTCa  >  1:210177/1‑140 (MQ=255)
                                                                                                                                       cGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATgg                                              >  2:133119/1‑93 (MQ=255)
                                                                                                                                       cGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATgg                                              <  1:133119/93‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
CGCGCAGTGAGCAGAGATCGTCAGAAAGCTCGCGAGGCAGCATAGGGATGTTGAAGCCAGGCAGATAGTTGGTGAATGCGCGAATTTTCGCGGCTTTGTCCAGCTTGCTGCCTTCAGCAATCCACGCGGTTGGATCGGCAATCGCCACAATCAGCTGAAGTTTGTCATCCGGCAACGCCTTAGCGAAAAGGGCGTCATCCATATCTTCTGTGCTGGCACTGTCGATGGTGACAAAATCCAGCGCGGTCAGATCTTCACGAACCAGACCTTCA  >  CP009273/1342347‑1342618

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