Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,909,758 C→A 100% E12* (GAA→TAA)  proQ ← RNA chaperone, putative ProP translation regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,909,7580CA96.4% 198.9 / 0.0 56E12* (GAA→TAA) proQRNA chaperone, putative ProP translation regulator
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (23/31);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCT  >  CP009273/1909614‑1909902
                                                                                                                                                |                                                                                                                                                
cGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTAc                                                                                                                                             >  1:151471/1‑150 (MQ=255)
cGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTAc                                                                                                                                             <  1:960885/150‑1 (MQ=255)
cGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTAc                                                                                                                                             <  1:384482/150‑1 (MQ=255)
       cAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTAtt                                                                                                                                      <  1:9270/150‑1 (MQ=255)
        aaTTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTAtt                                                                                                                                      >  2:633227/1‑149 (MQ=255)
          ttGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCaa                                                                                                                                   <  2:798346/150‑1 (MQ=255)
              gTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTa                                                                                                                               <  1:753300/150‑1 (MQ=255)
                 ttGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAgg                                                                                                                             <  2:984723/149‑1 (MQ=255)
                  tGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGt                                                                                                                            <  2:232956/149‑1 (MQ=255)
                   gCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGtt                                                                                                                           <  2:292812/149‑1 (MQ=255)
                    cTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTg                                                                                                                          <  2:247950/149‑1 (MQ=255)
                       aGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGAt                                                                                                                        <  1:619896/148‑1 (MQ=255)
                       aGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGAt                                                                                                                        <  1:239088/148‑1 (MQ=255)
                        ggTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGAtttt                                                                                                                     <  2:448815/150‑1 (MQ=255)
                            cATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCa                                                                                                                  <  1:546979/149‑1 (MQ=255)
                             aTTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCAt                                                                                                                 >  2:200688/1‑149 (MQ=255)
                             aTTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATg                                                                                                                <  2:476736/150‑1 (MQ=255)
                                  cccAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAAtt                                                                                                           >  1:453207/1‑150 (MQ=255)
                                  cccAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAAt                                                                                                            <  1:490482/149‑1 (MQ=255)
                                     aGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTcc                                                                                                        >  2:103021/1‑150 (MQ=255)
                                       cAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCt                                                                                                       >  1:589918/1‑149 (MQ=255)
                                               tCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAAc                                                                                              >  2:163169/1‑150 (MQ=255)
                                                   ccAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGAc                                                                                          <  2:467879/150‑1 (MQ=255)
                                                    cAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACg                                                                                         <  1:498653/150‑1 (MQ=255)
                                                     aaaTCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACg                                                                                         >  2:271052/1‑149 (MQ=255)
                                                      aaTCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACg                                                                                         >  1:767468/1‑148 (MQ=255)
                                                             aaaaaTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCcaaca                                                                                <  1:633227/150‑1 (MQ=255)
                                                              aaaaTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCcaacaa                                                                               <  1:295982/150‑1 (MQ=255)
                                                               aaaTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAg                                                                              <  1:914708/150‑1 (MQ=255)
                                                                    ccGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAg                                                                         >  1:794392/1‑150 (MQ=255)
                                                                        ttttCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCAt                                                                     >  2:737568/1‑150 (MQ=255)
                                                                            cAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGaa                                                                  >  1:699688/1‑149 (MQ=255)
                                                                                ggACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGc                                                             >  2:332919/1‑150 (MQ=255)
                                                                                      gcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCAt                                                        >  2:241050/1‑149 (MQ=255)
                                                                                       cTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATta                                                      >  2:912069/1‑150 (MQ=255)
                                                                                        ttCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATta                                                      >  2:658775/1‑149 (MQ=255)
                                                                                        ttCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATta                                                      >  1:466429/1‑149 (MQ=255)
                                                                                          cACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATtat                                                     <  2:97301/148‑1 (MQ=255)
                                                                                          cACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATaa                                                   >  1:1012537/1‑150 (MQ=255)
                                                                                           aCCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATaa                                                   <  1:482052/149‑1 (MQ=255)
                                                                                           aCCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAAt                                                  <  2:453207/150‑1 (MQ=255)
                                                                                                 cGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCt                                            <  1:737568/150‑1 (MQ=255)
                                                                                                 cGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCt                                            <  1:103021/150‑1 (MQ=255)
                                                                                                        aaaCAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTa                                     <  2:897169/150‑1 (MQ=255)
                                                                                                            aGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATg                                 >  2:175973/1‑150 (MQ=255)
                                                                                                                     aaCGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTAt                        <  2:876321/150‑1 (MQ=255)
                                                                                                                       cGTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATcc                      <  2:671897/150‑1 (MQ=255)
                                                                                                                        gTTCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCg                     >  1:558853/1‑150 (MQ=255)
                                                                                                                         ttCGGCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGt                    <  1:260986/150‑1 (MQ=255)
                                                                                                                            ggCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTg                  <  2:122067/149‑1 (MQ=255)
                                                                                                                            ggCCAGAAACGCGATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGa                 >  1:132745/1‑150 (MQ=255)
                                                                                                                                   aaCGCGATTGCTTATTTACTGCTATTCAACTTAGGTTGAATTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCt          <  2:70061/150‑1 (MQ=255)
                                                                                                                                     cgcgATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGc        >  1:973496/1‑150 (MQ=255)
                                                                                                                                       cgATTACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGa      >  2:845252/1‑150 (MQ=255)
                                                                                                                                          ttACTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGc   >  2:137574/1‑150 (MQ=255)
                                                                                                                                            aCTTATTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCt  <  1:175973/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGGATCGCAATTGCGTTTTGCTCAGGTTCATTTCCCCAGCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCATTATAATGACGCTATCAGTAAATGCTACGTTATCCGTTGATTATCCTGCGACGCT  >  CP009273/1909614‑1909902

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