Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,221,752 T→C pseudogene (2152/2648 nt) ycgH → putative transporter component, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,221,7520TC96.4% 93.2 / NA 28pseudogene (2152/2648 nt)ycgHputative transporter component, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (14/13);  minor base A (1/0);  total (15/13)
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

AGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGA  >  NC_000913/1221629‑1221878
                                                                                                                           |                                                                                                                              
aGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAgcg                                                                                                                       >  2:65686/1‑133 (MQ=255)
     aGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCt                                                                                                                             <  1:370862/122‑1 (MQ=255)
     aGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCt                                                                                                                             >  2:370862/1‑122 (MQ=255)
         cTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTAAGAAACCTATAGCGCGAAAGGTg                                                                                                              >  1:183430/1‑133 (MQ=255)
           gTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCa                                                                                                            <  1:65686/133‑1 (MQ=255)
                ccACCTGGCTTGCCGTTGCCATCAGTAAAACAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGt                                                                                                       <  1:161396/133‑1 (MQ=255)
                          tGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGaa                                                                                             <  1:43051/133‑1 (MQ=255)
                           gCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAg                                                                                            >  1:55480/1‑133 (MQ=255)
                             cGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCg                                                                                          <  2:316912/133‑1 (MQ=255)
                                                   gcATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATa                                                                    >  1:509824/1‑133 (MQ=255)
                                                   gcATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATa                                                                    >  1:118655/1‑133 (MQ=255)
                                                             aCTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTgg                                                          >  2:507378/1‑133 (MQ=255)
                                                               tCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCt                                                              >  1:384766/1‑127 (MQ=255)
                                                               tCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCt                                                              <  2:384766/127‑1 (MQ=255)
                                                                   gggCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGa                                                    >  1:97070/1‑133 (MQ=255)
                                                                   gggCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTAGAAGCGGGTTACGTTTTTGCCGTGACTACATCCTTAGGTCTGGa                                                    >  1:191421/1‑133 (MQ=255)
                                                                    ggCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGaa                                                   >  2:341403/1‑133 (MQ=255)
                                                                                 tGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAAtata                                      <  2:118655/133‑1 (MQ=255)
                                                                                     ttAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTg                                  >  2:348898/1‑133 (MQ=255)
                                                                                                        ggATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCcacaggca               <  2:55480/133‑1 (MQ=255)
                                                                                                         gATGAATTAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCcacaggcac              <  2:35161/133‑1 (MQ=255)
                                                                                                          aTGAATTAGCCTATGCACCCTATAGCGCAAAAGGTGCAACCATCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCcacaggcaca             <  2:183430/133‑1 (MQ=255)
                                                                                                           tGAATTAGCCCATGACCCCTATAGAGCGCAAGTTGCAACCGCCTCGCTAGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCcacaggcacag            <  2:191421/133‑1 (MQ=255)
                                                                                                                tAGCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTAt       <  1:493470/133‑1 (MQ=255)
                                                                                                                  gCCTATGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCt     >  1:22682/1‑133 (MQ=255)
                                                                                                                       tGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATc      <  1:80336/127‑1 (MQ=255)
                                                                                                                       tGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATc      >  2:80336/1‑127 (MQ=255)
                                                                                                                       tGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGa  <  1:487825/131‑1 (MQ=255)
                                                                                                                       tGAACCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGa  >  2:487825/1‑131 (MQ=255)
                                                                                                                           |                                                                                                                              
AGCGCAGGGCTGTATGCCACCTGGTTTGCCGATGACATCAGTAAAAAAGGCGCATACATTGACTCCTGGGCGCAATATAGCTGGTTTAAAAACTCGGTGAAAGGGGATGAATTAGCCTATGAATCCTATAGCGCGAAAGGTGCAACCGTCTCGCTGGAAGCGGGTTACGGTTTTGCCCTGAATAAATCCTTTGGTCTGGAAGCGGCGAAATATACGTGGATCTTCCAGCCACAGGCACAGGCTATCTGGA  >  NC_000913/1221629‑1221878

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