Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,210,374 C→T 8.1% pseudogene (1275/2001 nt) yehQ → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,210,3740CT 8.1% 146.9 / 3.6 49pseudogene (1275/2001 nt)yehQpseudogene
Reads supporting (aligned to +/- strand):  ref base C (28/17);  new base T (2/2);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

AGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGG  >  NC_000913/2210235‑2210522
                                                                                                                                           |                                                                                                                                                    
aGGAAAACTCACCGGCAGCTACGCGTCGGATGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAACCACCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCg                                                                                                                                             >  8:81742/1‑149 (MQ=255)
  gAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGTGTTGTGCCGCt                                                                                                                                           >  3:177135/1‑149 (MQ=255)
  gAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCt                                                                                                                                           >  7:124491/1‑149 (MQ=255)
  gAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGATACCCGCAACCGCTTAAGCAGTGTTGTGCCGCt                                                                                                                                           >  3:174064/1‑149 (MQ=255)
        tCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGCGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGcc                                                                                                                                     >  2:72308/1‑149 (MQ=255)
        tCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGcc                                                                                                                                     >  7:166002/1‑149 (MQ=255)
        tCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGcc                                                                                                                                     >  7:56785/1‑149 (MQ=255)
                        gTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGAt                                                                                                                                 <  7:52689/137‑1 (MQ=255)
                        gTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGAt                                                                                                                                 >  8:52689/1‑137 (MQ=255)
                        gTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATg                                                                                                                     >  2:67124/1‑149 (MQ=255)
                         tCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGt                                                                                                                    <  5:118836/149‑1 (MQ=255)
                           ggCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTg                                                                                                                  >  3:61847/1‑149 (MQ=255)
                                        tCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTa                                                                                                      <  7:128041/148‑1 (MQ=255)
                                                gcgcACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCg                                                                                             >  6:180409/1‑149 (MQ=255)
                                                 ctcACTTGCGGGCGGGAAAATTGTTTCACAAGCATCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGTGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCgg                                                                                            <  4:174064/147‑1 (MQ=255)
                                                 cgcACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGTGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCgg                                                                                            <  4:177135/149‑1 (MQ=255)
                                                         cgggcgggCAAATTGTTTCACAAGCAGCAAAACGCAGGGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATtgtg                                                                                    <  6:96863/149‑1 (MQ=255)
                                                                gcAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGc                                                                             <  8:56785/149‑1 (MQ=255)
                                                                 cAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCGTAAGCAGCGTTGTGCGGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCATCCGGGCATTGTGGCTTTGCg                                                                            <  2:134643/149‑1 (MQ=255)
                                                                    aTTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGa                                                                         >  1:130239/1‑149 (MQ=255)
                                                                    aTTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGa                                                                         <  5:180409/149‑1 (MQ=255)
                                                                       gTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCAGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAata                                                                      >  6:88537/1‑149 (MQ=255)
                                                                          tCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATAtt                                                                    <  3:230549/148‑1 (MQ=255)
                                                                                gcagcaAAACGCAGTGACGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGc                                                               >  3:68027/1‑147 (MQ=255)
                                                                                  agcaAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGc                                                           >  2:33404/1‑149 (MQ=255)
                                                                                    caAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGt                                                         >  6:155088/1‑149 (MQ=255)
                                                                                       aaCGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTccc                                                      >  6:143202/1‑149 (MQ=255)
                                                                                        aCGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTcccc                                                     <  4:61847/149‑1 (MQ=255)
                                                                                              tGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCt                                               <  7:81742/149‑1 (MQ=255)
                                                                                                 cGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCa                                            >  7:219979/1‑149 (MQ=255)
                                                                                                 cGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCa                                            >  4:80803/1‑149 (MQ=255)
                                                                                                         aGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTcc                                                       >  5:85431/1‑130 (MQ=255)
                                                                                                         aGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTcc                                                       <  6:85431/130‑1 (MQ=255)
                                                                                                           ctgctgTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCctct                                  >  3:65979/1‑149 (MQ=255)
                                                                                                              ctgTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTaa                               >  1:187818/1‑149 (MQ=255)
                                                                                                                 ttAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATc                             <  8:219979/148‑1 (MQ=255)
                                                                                                                    gCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGt                         <  1:67124/149‑1 (MQ=255)
                                                                                                                     cTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTc                        <  8:124491/149‑1 (MQ=255)
                                                                                                                         ccGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGAt                     >  5:69082/1‑148 (MQ=255)
                                                                                                                           gCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAAc                  >  4:176356/1‑149 (MQ=255)
                                                                                                                           gCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAAc                  >  7:181164/1‑149 (MQ=255)
                                                                                                                             aaCCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACt                 >  1:129516/1‑148 (MQ=255)
                                                                                                                               ccGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACttat              <  3:80803/149‑1 (MQ=255)
                                                                                                                               ccGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACttat              <  4:68027/149‑1 (MQ=255)
                                                                                                                                cGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACttatt             >  4:113999/1‑149 (MQ=255)
                                                                                                                                cGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACttatt             <  5:155088/149‑1 (MQ=255)
                                                                                                                                 gCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGTGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACttatt             >  5:172225/1‑148 (MQ=255)
                                                                                                                                   ttAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACttatttat          >  2:122376/1‑149 (MQ=255)
                                                                                                                                   ttAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACttatttat          >  1:25450/1‑149 (MQ=255)
                                                                                                                                          gcGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGACTGCATACGGCCTCTTAATCAGGTCGATAACTTATTTATTCTGCCg   <  6:69082/149‑1 (MQ=255)
                                                                                                                                           cGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACTTATTTATTCTGCCgg  <  6:36149/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
AGGAAAACTCACCGGCAGCTACGCGTCGGCTGTTTAGTTTTCAGGCTGGCGCACTGGCGGGCGGGCAAATTGTTTCACAAGCAGCAAAACGCAGTGCCGATGGCGAGCTGCTGTTAGCTACCCGCAACCGCTTAAGCAGCGTTGTGCCGCTGTCGCCTGATGCCTGGCAAATGTTGAGCGCGCCGTTACGCCAGCCGGGCATTGTGGCTTTGCGGGAATATTTACGCCAGCGTCCCCCCGCCTGCATACGGCCTCTTAATCAGGTCGATAACTTATTTATTCTGCCGG  >  NC_000913/2210235‑2210522

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