Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,130,088 T→G pseudogene (56/1101 nt) yghO ← pseudogene, IS‑interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,130,0880TG98.2% 194.9 / NA 55pseudogene (56/1101 nt)yghOpseudogene, IS‑interrupted
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (28/26);  minor base C (1/0);  total (29/26)
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.99e-01

GGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTC  >  NC_000913/3129943‑3130234
                                                                                                                                                 |                                                                                                                                                  
ggTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAgcgtcg                                                                                                                                               >  2:402196/1‑151 (MQ=255)
       aTTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATaa                                                                                                                                         <  1:635563/150‑1 (MQ=255)
       aTTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATaa                                                                                                                                         <  1:604149/150‑1 (MQ=255)
           gcgcGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGtg                                                                                                                                     <  2:589363/150‑1 (MQ=255)
           gcgcGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGCATCCAGTTTGGAGCGTCGGGATAAAGtg                                                                                                                                     >  2:2397/1‑150 (MQ=255)
               ggTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGAc                                                                                                                                <  2:666544/151‑1 (MQ=255)
                     aCGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGaaaa                                                                                                                          <  1:622780/151‑1 (MQ=255)
                     aCGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGaaaa                                                                                                                          <  1:581708/151‑1 (MQ=255)
                      cGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAg                                                                                                                         <  1:115724/151‑1 (MQ=255)
                           aaCTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATa                                                                                                                    <  1:92592/151‑1 (MQ=255)
                           aaCTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCAGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGAt                                                                                                                     >  2:373088/1‑150 (MQ=255)
                             cTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATaa                                                                                                                   >  2:62916/1‑150 (MQ=255)
                              tATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAAt                                                                                                                 >  2:442180/1‑151 (MQ=255)
                                     cccGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTa                                                                                                           >  2:229179/1‑150 (MQ=255)
                                         cTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTc                                                                                                      >  1:206057/1‑151 (MQ=255)
                                                  ccacccacGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTttatt                                                                                             <  1:844058/151‑1 (MQ=255)
                                                    acccacGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTttatta                                                                                            >  2:746681/1‑150 (MQ=255)
                                                    acccacGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTttatta                                                                                            >  1:698021/1‑150 (MQ=255)
                                                       cacGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGGCATTTTTATTaa                                                                                           <  1:526611/148‑1 (MQ=255)
                                                             tGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGc                                                                                  >  1:134704/1‑151 (MQ=255)
                                                              gCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTg                                                                                   <  1:316918/149‑1 (MQ=255)
                                                                cAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTTGTTGCGTTCAATGAATAAAGGGGGTATTCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCCCTTTTATTAAGGAAGGTGc                                                                                  <  1:793342/148‑1 (MQ=255)
                                                                   atgatATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACaa                                                                            <  2:113223/151‑1 (MQ=255)
                                                                       tatGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTcc                                                                        >  2:376731/1‑151 (MQ=255)
                                                                         tGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCt                                                                      >  2:338600/1‑151 (MQ=255)
                                                                              cgtccCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGatat                                                                 <  2:715567/151‑1 (MQ=255)
                                                                               gtccCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATg                                                                <  2:111301/151‑1 (MQ=255)
                                                                                tccCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATg                                                                <  2:530243/150‑1 (MQ=255)
                                                                                 ccccGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATg                                                                >  1:645914/1‑149 (MQ=255)
                                                                                 ccccGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATg                                                                >  2:742936/1‑149 (MQ=255)
                                                                                    cGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATg                                                                >  2:618008/1‑146 (MQ=255)
                                                                                    cGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGAATCCAGTTTGGACCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTACGGTCATTTTTATTAAGGAAAGTGCGAACAAGTCCCTGatat                                                                 >  1:432831/1‑145 (MQ=255)
                                                                                     ggATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATc                                                           >  1:248049/1‑150 (MQ=255)
                                                                                               cAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCAt                                                >  1:195137/1‑151 (MQ=255)
                                                                                                 gACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCt                                              >  2:127472/1‑151 (MQ=255)
                                                                                                  aCAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTg                                             >  1:439454/1‑151 (MQ=255)
                                                                                                   cAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTgg                                            <  2:340975/151‑1 (MQ=255)
                                                                                                    aaaTGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAGGTCCCTGATATGAGATCATGTTTGTCATCTgg                                            >  1:411984/1‑150 (MQ=255)
                                                                                                    aaaTGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTgg                                            >  2:62836/1‑150 (MQ=255)
                                                                                                     aaTGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAg                                          <  1:144482/151‑1 (MQ=255)
                                                                                                      aTGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGc                                         <  2:776495/151‑1 (MQ=255)
                                                                                                        gCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGcc                                        <  2:177052/150‑1 (MQ=255)
                                                                                                        gCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCa                                       <  2:248074/151‑1 (MQ=255)
                                                                                                               gCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACa                                <  1:499411/151‑1 (MQ=255)
                                                                                                                      aTGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtc                          >  2:279721/1‑150 (MQ=255)
                                                                                                                        gAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtcat                        <  2:63481/150‑1 (MQ=255)
                                                                                                                         aaTAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtcatc                       >  2:271736/1‑150 (MQ=255)
                                                                                                                             aaGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATg                    >  2:699928/1‑149 (MQ=255)
                                                                                                                                   gTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtcatca            <  2:432375/151‑1 (MQ=255)
                                                                                                                                   gTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtcatca            <  1:167568/151‑1 (MQ=255)
                                                                                                                                     aTCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAAc          >  1:38111/1‑151 (MQ=255)
                                                                                                                                       ccAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACtt        >  1:127683/1‑151 (MQ=255)
                                                                                                                                           tttGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCt    >  1:669212/1‑151 (MQ=255)
                                                                                                                                            ttGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCtt   >  2:129665/1‑151 (MQ=255)
                                                                                                                                             tGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCtt   <  1:727450/150‑1 (MQ=255)
                                                                                                                                             tGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTc  >  2:363294/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTC  >  NC_000913/3129943‑3130234

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