Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 314,816 (G)5→6 pseudogene (460/872 nt) eaeH → pseudogene, attaching and effacing protein homology;factor; Not classified

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA314,8111.G98.2% 190.1 / ‑3.2 54pseudogene (455/872 nt)eaeHpseudogene, attaching and effacing protein homology;factor; Not classified
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (26/27);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATCCCGTAGTCATACCCGCATTGGTGTT  >  NC_000913_3_pae_tpiA/314663‑314953
                                                                                                                                                     |                                                                                                                                              
aaagctaaCCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAg                                                                                                                                               >  2:429683/1‑150 (MQ=255)
 aagctaaCCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAGTGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAgg                                                                                                                                              <  2:554830/150‑1 (MQ=255)
      aaCCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAggggggc                                                                                                                                         >  2:40159/1‑148 (MQ=255)
          aGGAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATa                                                                                                                                      <  1:559282/150‑1 (MQ=255)
           ggAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATAc                                                                                                                                     <  3:506998/150‑1 (MQ=255)
               aTACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCg                                                                                                                                 <  3:184553/150‑1 (MQ=255)
                   aGGATTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATTTTGTTCACTCAGGGGGGCAATACATCGTACa                                                                                                                             <  4:514304/150‑1 (MQ=255)
                     gAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGa                                                                                                                           <  3:302529/150‑1 (MQ=255)
                           cTCGGGAAATATGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCg                                                                                                                     <  4:420428/150‑1 (MQ=255)
                                     atGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCa                                                                                                               <  4:98832/146‑1 (MQ=255)
                                     atGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCa                                                                                                               >  3:98832/1‑146 (MQ=255)
                                       ggTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCaaa                                                                                                         <  1:362216/150‑1 (MQ=255)
                                       ggTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCaaa                                                                                                         <  1:779126/150‑1 (MQ=255)
                                              cgcgcgTCAAACTGGATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCGGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTGAGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGtt                                                                                                 <  2:544557/150‑1 (MQ=255)
                                              cgcgcgTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGt                                                                                                  <  4:275115/150‑1 (MQ=255)
                                              cgcgcgTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAAGGCTTTATCCGATTTATGATACGACGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACGGACGATGGTACTCAGTTAAATATTGGt                                                                                                  >  2:415230/1‑150 (MQ=255)
                                               gcgcgTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGtt                                                                                                 <  1:41085/150‑1 (MQ=255)
                                               gcgcgTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGtt                                                                                                 >  2:564478/1‑150 (MQ=255)
                                                  cgTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATc                                                                                                                      <  4:249194/126‑1 (MQ=255)
                                                  cgTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATc                                                                                                                      >  3:249194/1‑126 (MQ=255)
                                                  cgTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTtg                                                                                              >  1:108138/1‑150 (MQ=255)
                                                  cgTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTtg                                                                                              <  3:504437/150‑1 (MQ=255)
                                                    tCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTtggc                                                                                            >  3:509216/1‑150 (MQ=255)
                                                     cAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTtggct                                                                                           >  1:291948/1‑150 (MQ=255)
                                                      aaaCTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTtggctg                                                                                          <  1:691244/150‑1 (MQ=255)
                                                        aCTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTtggctggc                                                                                        >  1:341419/1‑150 (MQ=255)
                                                               gTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAggg                                                                                                                                             >  2:744412/1‑89 (MQ=255)
                                                               gTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAggg                                                                                                                                             <  1:744412/89‑1 (MQ=255)
                                                                   aTAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCa                                                                             <  1:427006/150‑1 (MQ=255)
                                                                   aTAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCa                                                                             >  2:361003/1‑150 (MQ=255)
                                                                       aGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGaa                                                                         >  1:769975/1‑150 (MQ=255)
                                                                               cGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTgg                                                                 >  2:191323/1‑150 (MQ=255)
                                                                                  tGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATg                                                              >  2:166199/1‑150 (MQ=255)
                                                                                    aaGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGc                                                            >  3:315444/1‑150 (MQ=255)
                                                                                     aGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCg                                                           <  3:105206/150‑1 (MQ=255)
                                                                                      ggATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTAAGGAAATGACTGGATGGCgg                                                          >  3:373872/1‑150 (MQ=255)
                                                                                       gATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCggg                                                         <  2:275935/150‑1 (MQ=255)
                                                                                         ttcttcGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCggggg                                                       <  2:414072/150‑1 (MQ=255)
                                                                                          tcttcGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGt                                                      <  2:499509/150‑1 (MQ=255)
                                                                                          tcttcGCTGGAAATGCATTATCCGATTTATGATACGCCGACAAATCTGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGt                                                      <  1:463056/150‑1 (MQ=255)
                                                                                           cttcGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTg                                                     >  2:792508/1‑150 (MQ=255)
                                                                                              cGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAAt                                                  <  1:40159/150‑1 (MQ=255)
                                                                                               gCTGGAAATGCTTCATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTCTTCAGGAAATGACTGGATGGCTGGGGTGAATa                                                 >  2:661330/1‑150 (MQ=255)
                                                                                                     aaTGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTa                                           <  3:508001/150‑1 (MQ=255)
                                                                                                          tttATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGAt                                      >  2:83398/1‑150 (MQ=255)
                                                                                                           ttATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACATATTTTCGTTTTCTTTTAAATATTGGTTTTGGCTGGCGTCATTTTTTAGTAAATGACTGGATGGCGGGGGTGAATTCTTTTATCGATc                                     >  4:452206/1‑150 (MQ=255)
                                                                                                              tCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATg                                  <  4:271954/150‑1 (MQ=255)
                                                                                                                agATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGAt                                <  2:137128/149‑1 (MQ=255)
                                                                                                                      aTGATACGCCGACAAATATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATcc                          >  2:675016/1‑150 (MQ=255)
                                                                                                                         aTACGCCGACAAATATGTTGTTCACTCAGGGGGGCACTACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATCCCGt                       >  4:415746/1‑150 (MQ=255)
                                                                                                                                   aaaTATGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATCCCGTAGTCATAccc             <  4:235061/150‑1 (MQ=255)
                                                                                                                                      tatGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATCCCGTAGTCATACCCGCa          <  1:792508/150‑1 (MQ=255)
                                                                                                                                       atGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATCCCGTAGTCATACCCGCAt         >  4:435687/1‑150 (MQ=255)
                                                                                                                                       atGTTGTTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATCCCGTAGTCATACCCGCAt         >  4:497908/1‑150 (MQ=255)
                                                                                                                                        tgttgtTCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATCCCGTAGTCATACCCGCAtt        >  3:341291/1‑150 (MQ=255)
                                                                                                                                         gttgttCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATCCCGTAGTCATACCCGCAttg       >  1:806762/1‑150 (MQ=255)
                                                                                                                                            gttCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATCCCGTAGTCATACCCGCAttggtg    >  1:603501/1‑150 (MQ=255)
                                                                                                                                             ttCACTCAGGGGGGCAATACATCGTCCAGACGATCGTACTCAGTCAAATATTGGTCATGGCAGGCGTCATTTTTAAGGAAATGTCTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATCCCGTAGTCATACCCGCAttggtgt   <  1:415230/150‑1 (MQ=255)
                                                                                                                                              tCACTCAGGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATCCCGTAGTCATACCCGCAttggtgtt  >  2:127908/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
AAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAATGTCGATAAAGATTTCTCGCTGAAGGATTCTTCGCTGGAAATGCTTTATCCGATTTATGATACGCCGACAAATATGTTGTTCACTCAGGGGGCAATACATCGTACAGACGATCGTACTCAGTCAAATATTGGTTTTGGCTGGCGTCATTTTTCAGGAAATGACTGGATGGCGGGGGTGAATACTTTTATCGATCATGATTTATCCCGTAGTCATACCCGCATTGGTGTT  >  NC_000913_3_pae_tpiA/314663‑314953

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