Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 647,430 C→T R27H (CGC→CAC)  citG ← 2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913647,4300CT88.6% 98.5 / 4.6 36R27H (CGC→CAC) citG2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (15/16);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGCGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGG  >  NC_000913/647212‑647653
                                                                                                                                                                                                                          |                                                                                                                                                                                                                               
cTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGt                                                                                                                                                                                                               >  3:354696/1‑237 (MQ=255)
   cGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCa                                                                                                                                                                                                          <  3:473823/239‑1 (MQ=255)
            tACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGCGCCAGCCCAGCAGGGCGTACTCa                                                                                                                                                                                                          <  1:185494/230‑1 (MQ=255)
                  ttGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGccagcccagc                                                                                                                                                                                                                      <  4:220013/212‑1 (MQ=255)
                  ttGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGccagcccagc                                                                                                                                                                                                                      >  3:220013/1‑212 (MQ=255)
                  ttGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGccag                                                                                                                                                                                                                            >  3:278661/1‑206 (MQ=255)
                      gCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGCGCCAGCCCACCCGGGAGAACGCATCGATTAATg                                                                                                                                                                                                >  2:185494/1‑230 (MQ=255)
                         tAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCt                                                                                                                                                                                    >  3:450012/1‑239 (MQ=255)
                                        cGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTg                                                                                                                                                                     <  3:158476/239‑1 (MQ=255)
                                                        tCCGCACTCCAGGCACAAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGa                                                                                                                                                     <  4:363920/239‑1 (MQ=255)
                                                          cGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACa                                                                                                                                                   <  3:245888/239‑1 (MQ=255)
                                                             aCTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCg                                                                                                                                                >  3:356360/1‑239 (MQ=255)
                                                                     cACCAAATTCAATGAACCGGGGTACCCAGCCCGGAATCGCCAGCGCGCTGCGGTGGAAATCCTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTcag                                                                                                                                           <  4:356360/236‑1 (MQ=255)
                                                                                aTGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTTTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGGGCGTCCAttc                                                                                                                             >  4:246095/1‑237 (MQ=255)
                                                                                 tGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCg                                                                                                                            <  3:246095/239‑1 (MQ=255)
                                                                                    aaCGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCa                                                                                                                                                                                        <  3:281179/176‑1 (MQ=255)
                                                                                    aaCGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCa                                                                                                                                                                                        >  4:281179/1‑176 (MQ=255)
                                                                                                                      gCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCAcgc                                                                                       <  2:64985/239‑1 (MQ=255)
                                                                                                                                                     tgCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATg                                                                                                                                                                                                >  3:357027/1‑103 (MQ=255)
                                                                                                                                                     tgCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATg                                                                                                                                                                                                <  4:357027/103‑1 (MQ=255)
                                                                                                                                                       cgcACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAg                                                       >  1:200355/1‑238 (MQ=255)
                                                                                                                                                          aCCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACa                                                                                                                                                   >  2:187957/1‑143 (MQ=255)
                                                                                                                                                          aCCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACa                                                                                                                                                   <  1:187957/143‑1 (MQ=255)
                                                                                                                                                          aCCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGCGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGtt                                                                                                                         <  4:80903/169‑1 (MQ=255)
                                                                                                                                                          aCCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGCGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGtt                                                                                                                         >  3:80903/1‑169 (MQ=255)
                                                                                                                                                                      gCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGAt                                                                                                     >  4:213168/1‑177 (MQ=255)
                                                                                                                                                                      gCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGAt                                                                                                     <  3:213168/177‑1 (MQ=255)
                                                                                                                                                                         aTCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAgg                                           >  3:284970/1‑232 (MQ=255)
                                                                                                                                                                         aTCCAAGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAgg                                           <  4:284970/232‑1 (MQ=255)
                                                                                                                                                                                                 cAGATTGACTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAAtt            >  3:431380/1‑239 (MQ=255)
                                                                                                                                                                                                        aCTTCAGTCAGCATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCtt     <  4:431380/239‑1 (MQ=255)
                                                                                                                                                                                                           tcaggcagCAGTGCGGGCCAGCCCAGACGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCgg  <  4:450012/239‑1 (MQ=255)
                                                                                                                                                                                                                   cATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGt                                                                                                 <  4:199489/136‑1 (MQ=255)
                                                                                                                                                                                                                   cATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGt                                                                                                 >  3:199489/1‑136 (MQ=255)
                                                                                                                                                                                                                   cATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGa                                                                                <  4:301554/153‑1 (MQ=255)
                                                                                                                                                                                                                   cATGGCGTGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGa                                                                                >  3:301554/1‑153 (MQ=255)
                                                                                                                                                                                                                       gcgtgcCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAgg                                           >  3:456551/1‑186 (MQ=255)
                                                                                                                                                                                                                       gcgtgcCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGATTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAgg                                           <  4:456551/186‑1 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                               
CTTCGCAAGCCATACCAATTGGGCGTAATCCGTGGAGTACCGCTTCTGGTGCCATTTCCGCACTACAGGCACCAAATTCAATGAAACGGGGTAGCCAGCCCTGAATCGCCAGCGCGCTGCGGTGGAAATCTTCCAGCGCCATATCTTTGTGCGCACCGCAGTTAATGCGATCCACGAGGCCTGGTTTCGGTGACAGATTGACTTCAGTCAGCATGGCGCGCCAGCCCAGCAGGGCGTACTCATCGATTAATGACGTCGCAAGCTTTGTGGTTTTAGTTGACGTTGCAGGCATCGACATCGTTCAGCAGTGCCTCCATGCGGTTGAGTAAATCGGTCAGTTGATGGGTTTTTCCACGCGCGCAGACGGCTGCGCTTTGTTCGCACAACAGGCAGCGGCGAGGCGGCAGTGAATAGTCGCGGCGGGAGAGAATTTCGCCTTCGG  >  NC_000913/647212‑647653

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: