Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 528,172 (G)6→5 pseudogene (224/711 nt) ylbH → putative RHS domain‑containing protein YlbH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913528,1670G.94.6% 145.0 / 2.7 37pseudogene (219/711 nt)ylbHputative RHS domain‑containing protein YlbH
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (15/20);  total (16/21)
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

ACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATA  >  NC_000913/528031‑528296
                                                                                                                                        |                                                                                                                                 
aCCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAgggg                                                                                                                                >  2:116435/1‑140 (MQ=255)
  ccGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAggggga                                                                                                                              <  1:116435/140‑2 (MQ=255)
    gCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAgggggatg                                                                                                                            <  2:262832/140‑4 (MQ=255)
     cATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAgggggatgg                                                                                                                           <  2:79746/140‑5 (MQ=255)
      aTCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGa                                                                                                                         <  2:173963/140‑1 (MQ=255)
        cACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGt                                                                                                                       <  1:291235/140‑1 (MQ=255)
        cACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGt                                                                                                                       <  1:217608/140‑1 (MQ=255)
        cACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGt                                                                                                                       <  1:387572/140‑1 (MQ=255)
                   gCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTGTGCGTa                                                                                                            <  2:51688/140‑1 (MQ=255)
                   gCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTa                                                                                                            >  1:299187/1‑140 (MQ=255)
                    ccGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTAt                                                                                                           <  1:249861/140‑1 (MQ=255)
                    ccGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTAt                                                                                                           >  1:330617/1‑140 (MQ=255)
                             cTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGaa                                                                                                   >  1:140009/1‑140 (MQ=255)
                                     gcagcaGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGaa                                                                                          <  2:391259/140‑1 (MQ=255)
                                     gcagcaGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCAGTATGCGTATCCGCTGAATCCGGTGaa                                                                                          <  2:135389/140‑1 (MQ=255)
                                          agTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTa                                                                                     <  2:136289/140‑1 (MQ=255)
                                              tGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGaa                                                                                          >  2:304885/1‑131 (MQ=255)
                                              tGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGaa                                                                                          <  1:304885/131‑1 (MQ=255)
                                                aTAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGCATTGATc                                                                               <  2:287265/140‑1 (MQ=255)
                                                                 tactacAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTccc                                                              >  1:25021/1‑140 (MQ=255)
                                                                  actacAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCg                                                             <  2:330617/140‑1 (MQ=255)
                                                                           ggaaccggTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAg                                                    >  1:98771/1‑140 (MQ=255)
                                                                                cggTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCta                                               <  2:302388/140‑1 (MQ=255)
                                                                                  gTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCtaat                                             >  1:386292/1‑140 (MQ=255)
                                                                                   tactacGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCtaata                                            >  2:124021/1‑140 (MQ=255)
                                                                                          aTCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAgaagaa                                     >  2:53573/1‑140 (MQ=255)
                                                                                             cGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAg                                  <  1:53573/140‑1 (MQ=255)
                                                                                             cGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAg                                  >  1:376961/1‑140 (MQ=255)
                                                                                                   aGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCaa                             <  2:140009/140‑1 (MQ=255)
                                                                                                               tCACTCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCaaa                >  1:62966/1‑140 (MQ=255)
                                                                                                                   tCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGc            <  2:299187/140‑1 (MQ=255)
                                                                                                                   tCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGc            <  1:52564/140‑1 (MQ=255)
                                                                                                                   tCAGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGc            >  2:365711/1‑140 (MQ=255)
                                                                                                                     aGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGc            <  2:109462/138‑1 (MQ=255)
                                                                                                                     aGGACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGc            >  1:109462/1‑138 (MQ=255)
                                                                                                                       gACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATgg        <  2:25021/140‑1 (MQ=255)
                                                                                                                       gACCCGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATgg        >  1:314003/1‑140 (MQ=255)
                                                                                                                          ccGATAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGat     >  2:188689/1‑140 (MQ=255)
                                                                                                                             aTAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCa           <  2:314766/131‑1 (MQ=255)
                                                                                                                             aTAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCa           >  1:314766/1‑131 (MQ=255)
                                                                                                                             aTAGGGCTGGAGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGatata  <  1:124021/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                 
ACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCGGAACCGGTACTACGATCCGTTGCAGGGGCGGTATATCACTCAGGACCCGATAGGGCTGGAGGGGGGATGGAGTCTGTATGCGTATCCGCTGAATCCGGTGAATGGTATTGATCCATTAGGGTTAAGTCCCGCAGATGTAGCGCTAATAAGAAGAAAAGATCAACTAAACCATCAAAGAGCATGGGATATA  >  NC_000913/528031‑528296

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