Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,210,811 C→T pseudogene (1712/2001 nt) yehQ → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,210,8110CT94.8% 246.8 / 2.8 77pseudogene (1712/2001 nt)yehQpseudogene
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base T (34/39);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 5.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTG  >  NC_000913/2210714‑2210907
                                                                                                 |                                                                                                
gTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGc                                                                                                <  2:1717451/100‑1 (MQ=255)
 tCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCt                                                                                               <  2:486345/100‑1 (MQ=255)
     ttttGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCt                                                                                           <  1:1029741/100‑1 (MQ=255)
     ttttGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCtt                                                                                          <  2:1054746/100‑1 (MQ=255)
        tgttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCt                                                                                        >  1:33698/1‑100 (MQ=255)
          ttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGc                                                                                      <  2:656610/100‑1 (MQ=255)
          ttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGc                                                                                      >  1:1562295/1‑100 (MQ=255)
            gAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCa                                                                                    >  2:1363009/1‑100 (MQ=255)
               nggCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCg                                                                                 >  1:212447/2‑100 (MQ=255)
               gggCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCCTCTGCCAGCg                                                                                 >  1:61398/1‑100 (MQ=255)
                  cAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGttt                                                                              <  2:32895/100‑1 (MQ=255)
                    aTTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTg                                                                            <  1:1556835/100‑1 (MQ=255)
                     ttGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGc                                                                           <  2:212448/100‑1 (MQ=255)
                     ttGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGc                                                                           >  2:363624/1‑100 (MQ=255)
                      tGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGcc                                                                          <  2:641640/100‑1 (MQ=255)
                       gACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGccn                                                                         <  1:771056/100‑2 (MQ=255)
                       gACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCt                                                                         >  2:496453/1‑100 (MQ=255)
                         cacTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgt                                                                       >  2:1873421/1‑100 (MQ=255)
                         cacTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgt                                                                       >  1:1613258/1‑100 (MQ=255)
                         cacTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgt                                                                       >  1:1102210/1‑100 (MQ=255)
                             ggAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGcc                                                                    >  1:1681554/1‑99 (MQ=255)
                             ggAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCg                                                                   <  2:1537877/100‑1 (MQ=255)
                             ggAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCg                                                                   <  2:339673/100‑1 (MQ=255)
                                  cACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTAc                                                              <  2:1387695/100‑1 (MQ=255)
                                   aCGGGTGATGATGACAAAACCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCTTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACt                                                             <  2:889903/100‑1 (MQ=255)
                                   aCGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACt                                                             <  2:1156744/100‑1 (MQ=255)
                                      ggTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCt                                                          <  1:1703622/100‑1 (MQ=255)
                                       gTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtc                                                         <  1:580614/100‑1 (MQ=255)
                                           tgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtcatca                                                     >  2:278537/1‑100 (MQ=255)
                                            gatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTCCCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAg                                                    >  1:938289/1‑100 (MQ=255)
                                             atgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGt                                                   >  1:1458192/1‑100 (MQ=255)
                                               gaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTg                                                 <  1:1264989/100‑1 (MQ=255)
                                                aCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGc                                                >  2:1624452/1‑100 (MQ=255)
                                                aCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGc                                                >  1:1137322/1‑100 (MQ=255)
                                                 cAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                               <  2:216345/100‑1 (MQ=255)
                                                 cAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                               >  2:1382157/1‑100 (MQ=255)
                                                 cAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                               >  1:1522422/1‑100 (MQ=255)
                                                 cAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                               >  2:252067/1‑100 (MQ=255)
                                                 cAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                               >  1:1360824/1‑100 (MQ=255)
                                                  aaaaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTg                                              <  1:783633/100‑1 (MQ=255)
                                                  aaaaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTg                                              >  2:1453746/1‑100 (MQ=255)
                                                  aaaaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTg                                              >  2:1878108/1‑100 (MQ=255)
                                                   aaaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGa                                             >  1:1461154/1‑100 (MQ=255)
                                                   aaaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGa                                             <  2:1321798/100‑1 (MQ=255)
                                                     aaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATa                                           <  1:29687/100‑1 (MQ=255)
                                                     aaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCCGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATa                                           <  1:1438450/100‑1 (MQ=255)
                                                             ccTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCa                                   >  1:993649/1‑100 (MQ=255)
                                                             ccTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCa                                   >  1:411162/1‑100 (MQ=255)
                                                             ccTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCa                                   >  2:1462944/1‑100 (MQ=255)
                                                              cTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAg                                  >  1:1434915/1‑100 (MQ=255)
                                                              cTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAg                                  >  1:698866/1‑100 (MQ=255)
                                                              cTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAg                                  <  1:496451/100‑1 (MQ=255)
                                                               tGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAgg                                 >  1:1081444/1‑100 (MQ=255)
                                                                 ggCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCg                               <  2:336004/100‑1 (MQ=255)
                                                                  gcgcTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGt                              <  1:1788757/100‑1 (MQ=255)
                                                                   cgcTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGtt                             <  1:1556702/100‑1 (MQ=255)
                                                                         aCGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTAt                       <  1:1830293/100‑1 (MQ=255)
                                                                         aCGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTAt                       >  1:819808/1‑100 (MQ=255)
                                                                            cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                    >  2:1346217/1‑100 (MQ=255)
                                                                            cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                    >  2:1788184/1‑100 (MQ=255)
                                                                            cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                    >  2:254863/1‑100 (MQ=255)
                                                                            cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                    <  1:373298/100‑1 (MQ=255)
                                                                              gAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAAc                  <  2:1276236/100‑1 (MQ=255)
                                                                              gAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAAc                  <  1:116575/100‑1 (MQ=255)
                                                                               aaaCTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACt                 <  1:363622/100‑1 (MQ=255)
                                                                                aaCTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTg                >  2:1079665/1‑100 (MQ=255)
                                                                                aaCTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTg                >  1:1430926/1‑100 (MQ=255)
                                                                                  cTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCt              <  1:1211017/100‑1 (MQ=255)
                                                                                   tACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTc             >  1:251032/1‑100 (MQ=255)
                                                                                   tACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTc             <  2:1618227/100‑1 (MQ=255)
                                                                                    aCGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTcc            <  2:1102211/100‑1 (MQ=255)
                                                                                    aCGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTcc            <  2:1613268/100‑1 (MQ=255)
                                                                                      gCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCAt          <  2:1522430/100‑1 (MQ=255)
                                                                                       ccTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATa         <  2:1623385/100‑1 (MQ=255)
                                                                                       ccTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATa         <  2:1461162/100‑1 (MQ=255)
                                                                                       ccTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATa         >  1:511924/1‑100 (MQ=255)
                                                                                              cACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTg  >  1:1829745/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTG  >  NC_000913/2210714‑2210907

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