Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,632,997 (T)5→6 pseudogene (1462/1607 nt) yhiL ← putative uncharacterized protein YhiL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,632,9921.T89.8% 148.0 / 9.3 49pseudogene (1467/1607 nt)yhiLputative uncharacterized protein YhiL
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base T (22/22);  total (25/24)
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

GTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATA‑TT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAAATAATTAGCATATACAGGTGTGGGTTTC  >  NC_000913/3632860‑3633119
                                                                                                                                     |                                                                                                                                  
gTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCata                                                                                                                                     >  2:81349/1‑133 (MQ=255)
        caacTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATA‑tt‑t‑tttgt                                                                                                                          <  2:353144/133‑4 (MQ=255)
          acTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATA‑tt‑t‑tttgtct                                                                                                                        <  1:101140/133‑6 (MQ=255)
               gCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGAAGCCCATATTT‑T‑TT‑GTCTAATGt                                                                                                                   >  1:233129/1‑133 (MQ=255)
                cccTGAGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGTCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTa                                                                                                                  >  2:501569/1‑133 (MQ=255)
                cccTGAGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTa                                                                                                                  >  1:286330/1‑133 (MQ=255)
                   tGAGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGAGCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTAttt                                                                                                               >  2:350285/1‑133 (MQ=255)
                     aGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACACCGGGCCCATGACGCCAATA‑TT‑T‑TTCGTCTAATGTATTTGc                                                                                                             >  2:220872/1‑133 (MQ=255)
                     aGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACACCGGGCCAATGAGGCCTATA‑TT‑T‑TTAGTCTAATGTATTTGc                                                                                                             >  2:183749/1‑133 (MQ=255)
                          tAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTgggg                                                                                                        >  2:373170/1‑133 (MQ=255)
                          tAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTgggg                                                                                                        >  1:445786/1‑133 (MQ=255)
                             attattttataatttcTATCACGTACTTTCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCt                                                                                                     <  2:210421/133‑1 (MQ=255)
                              ttattttattattttTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATAATCTCCCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCtt                                                                                                    <  2:7128/133‑1 (MQ=255)
                                 ttttattattttTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATA‑TTCT‑TT‑GTCTAATGAATGTGCTGGGGTCTTCAc                                                                                                 >  1:210421/1‑133 (MQ=255)
                                        attttTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAg                                                                                          <  1:367140/133‑1 (MQ=255)
                                           tttATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAgtgt                                                                                       <  2:286330/133‑1 (MQ=255)
                                               tCACGTGCTTCCTTCCGGAAATGTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTCCCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTcc                                                                                   <  1:232533/133‑1 (MQ=255)
                                                cACGTGCTACCTTCATGGCATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTccc                                                                                  <  1:144536/133‑1 (MQ=255)
                                                  cGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCca                                                                                <  1:466852/133‑1 (MQ=255)
                                                      cttccttcAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTGCACCTTTGGGTGTATCCCCACata                                                                            >  2:134210/1‑133 (MQ=255)
                                                             cAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGAt                                                                     >  1:103697/1‑133 (MQ=255)
                                                                       aGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGt                                                           >  2:399963/1‑133 (MQ=255)
                                                                        gCCCCATACTACCATTATCAAACTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTg                                                          <  1:220872/133‑1 (MQ=255)
                                                                        gCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTg                                                          <  2:103697/133‑1 (MQ=255)
                                                                               aataCCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCt                                                   <  1:496107/133‑1 (MQ=255)
                                                                                      ttATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCa                                            >  1:345744/1‑133 (MQ=255)
                                                                                      ttATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCa                                            >  1:62790/1‑133 (MQ=255)
                                                                                       tATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCag                                           >  1:45075/1‑133 (MQ=255)
                                                                                       tATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCag                                           >  2:403998/1‑133 (MQ=255)
                                                                                       tATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATAGGAAAAAATGTTCACGTCCTGGCTTGAAATCag                                           >  2:232629/1‑133 (MQ=255)
                                                                                       tATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGGCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCag                                           >  2:7647/1‑133 (MQ=255)
                                                                                        aTCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCGCCACATATCAAGATACGTTCACGTGATGGTCTGAAATCagt                                          >  2:191867/1‑133 (MQ=255)
                                                                                        aTCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGaa                                                <  2:401671/127‑1 (MQ=255)
                                                                                        aTCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGaa                                                >  1:401671/1‑127 (MQ=255)
                                                                                               tgttgATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGGTCACGTGGTGGTCTGAAATCAATAGTGTTa                                   >  1:142666/1‑133 (MQ=255)
                                                                                                 ttgATTATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTGACGTGCTGGTCAGAAATCAGTAGTGTTATc                                 >  1:205732/1‑133 (MQ=255)
                                                                                                      tATGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAaat                            <  2:62790/133‑1 (MQ=255)
                                                                                                        tGATCTACCCCGGGTCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAaataa                          <  1:501569/133‑1 (MQ=255)
                                                                                                        tGATCTACCCCGGGCTCATGATGCACATC‑TT‑TGTT‑GTCTAAGGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAaataa                          <  2:223378/133‑1 (MQ=255)
                                                                                                        tGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAaataa                          <  1:481709/133‑1 (MQ=255)
                                                                                                        tGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAaataa                          <  1:435288/133‑1 (MQ=255)
                                                                                                        tGATCTACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAaataa                          <  1:399963/133‑1 (MQ=255)
                                                                                                           tCTTCCCCGGGCCCATGATTCCCATATTT‑T‑TT‑GTCTAATGTATTTTCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAAATAATTa                       <  2:44637/133‑1 (MQ=255)
                                                                                                             tACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAAATAATTAGc                     <  1:81349/133‑1 (MQ=255)
                                                                                                             tACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAAATAATTAGc                     <  2:345744/133‑1 (MQ=255)
                                                                                                             tACCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAAATAATTAGc                     <  1:318711/133‑1 (MQ=255)
                                                                                                              aCCCCGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAAATAATTAGCa                    >  1:380434/1‑133 (MQ=255)
                                                                                                               ccccGGGCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAAATAATTAGCat                   >  2:353875/1‑133 (MQ=255)
                                                                                                                    ggCCCATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAAATAATTAGCATATACa              >  2:91977/1‑133 (MQ=255)
                                                                                                                       ccATGATGCCCATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAAATAATTAGCATATACAGgt           <  1:363975/133‑1 (MQ=255)
                                                                                                                              gtccATA‑TTCT‑TT‑GTCTAATGTATGTCCGGGGGTCTTCACCCTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAAATAATTAGCATATACAGGTGTGGGtt    <  2:233129/131‑1 (MQ=255)
                                                                                                                                ccATATTT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAAATAATTAGCATATACAGGTGTGGGTTTc  <  1:350285/133‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                  
GTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCAGGGAATTTAGCCCCATAATACCATTATCAAATTGTTGATTATGATCTACCCCGGGCCCATGATGCCCATA‑TT‑T‑TT‑GTCTAATGTATTTGCTGGGGTCTTCACCTTTGAGTGTCTCCCCACATATGAAGATACGTTCACGTGCTGGTCTGAAATCAGTAGTGTTATCAAAATAATTAGCATATACAGGTGTGGGTTTC  >  NC_000913/3632860‑3633119

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