Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,535,441 C→T pseudogene (439/1029 nt) yjhR → pseudogene, helicase family, putative frameshift suppressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,535,4410CT98.3% 187.9 / ‑3.7 60pseudogene (439/1029 nt)yjhRpseudogene, helicase family, putative frameshift suppressor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/31);  minor base . (1/0);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

TCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGCCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCA  >  NC_000913/4535319‑4535565
                                                                                                                          |                                                                                                                             
tCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCcgcg                                                                                                               <  2:287679/139‑1 (MQ=255)
  tACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAg                                                                                                             >  2:124297/1‑139 (MQ=255)
       aaaaCATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTa                                                                                                        <  2:735102/139‑1 (MQ=255)
         accATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCTTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGTGATATGGACCTTATCGAGATGCAGTCCGCATTTTTGCCGCGAGGGTTAtt                                                                                                      >  1:1415054/3‑139 (MQ=255)
         aaCATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTAtt                                                                                                      >  1:62173/1‑139 (MQ=255)
          aCATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTg                                                                                                     >  1:1567519/1‑139 (MQ=255)
             tGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCa                                                                                                  <  2:1415054/139‑1 (MQ=255)
              gAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTcgc                                                                                                                   <  2:1499981/121‑1 (MQ=255)
              gAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTcgc                                                                                                                   >  1:1499981/1‑121 (MQ=255)
                      cgggTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCaa                                                                                                 <  2:674906/129‑1 (MQ=255)
                      cgggTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCaa                                                                                                 >  1:674906/3‑131 (MQ=255)
                      cAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATAtctc                                                                                        >  1:943945/1‑139 (MQ=255)
                       aGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATctct                                                                                       <  2:227975/139‑1 (MQ=255)
                          tttCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtttt                                                                                    <  2:62173/139‑1 (MQ=255)
                                   aGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAatat                                                                                            >  2:406601/1‑123 (MQ=255)
                                   aGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAatat                                                                                            <  1:406601/123‑1 (MQ=255)
                                   aGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGata                                                                           >  2:1111412/1‑139 (MQ=255)
                                    gCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGataa                                                                          <  2:1567519/139‑1 (MQ=255)
                                        cAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCaaaa                                                                                               >  2:1343423/1‑115 (MQ=255)
                                        cAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCaaaa                                                                                               <  1:1343423/115‑1 (MQ=255)
                                                  cTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAAt                                                            <  1:124297/139‑1 (MQ=255)
                                                     aaTGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCaaaa                                                                                               >  2:1032341/1‑102 (MQ=255)
                                                     aaTGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCaaaa                                                                                               <  1:1032341/102‑1 (MQ=255)
                                                     aaTGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATctct                                                                                       <  2:174518/109‑1 (MQ=255)
                                                     aaTGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATGTTCGCCGCGAGGGTTATTGGCAAAATATctct                                                                                       >  1:174518/1‑109 (MQ=255)
                                                          tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttcttc                                                                                <  1:1550152/111‑1 (MQ=255)
                                                          tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttcttc                                                                                >  2:1550152/1‑111 (MQ=255)
                                                                 tcACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCGCAATTTTGCCGGCGGGGGTATTTGCCAAATATCCTCTTTTTTCAGGTTATTATTGGTTTCAATTTGaactgttcccttt                                            >  2:1536358/1‑127 (MQ=255)
                                                                    cgcgCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAatat                                                                                            <  1:1501631/90‑1 (MQ=255)
                                                                    cgcgCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAatat                                                                                            >  2:1501631/1‑90 (MQ=255)
                                                                    cgcgCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGAc                                          >  1:743517/1‑139 (MQ=255)
                                                                    cgcgCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGAc                                          <  2:943945/139‑1 (MQ=255)
                                                                           aaGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTgg                                                                                                    >  1:1486772/1‑75 (MQ=255)
                                                                           aaGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTgg                                                                                                    <  2:1486772/75‑1 (MQ=255)
                                                                                   tttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCg                           <  1:1111412/139‑1 (MQ=255)
                                                                                     tCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTcgcc                                                                                                                  >  1:1662071/1‑51 (MQ=255)
                                                                                     tCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTcgcc                                                                                                                  <  2:1662071/51‑1 (MQ=255)
                                                                                             ttGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATAcaca                 <  2:657347/139‑1 (MQ=255)
                                                                                              tGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACaa                >  1:905378/1‑139 (MQ=255)
                                                                                               ggCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGAc                                          <  1:24950/112‑1 (MQ=255)
                                                                                               ggCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGAc                                          >  2:24950/1‑112 (MQ=255)
                                                                                                  gATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAAttt            <  2:905378/139‑1 (MQ=255)
                                                                                                    tatGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCa          >  1:582760/1‑139 (MQ=255)
                                                                                                            ttATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa  <  2:582760/139‑1 (MQ=255)
                                                                                                             tATCGAGATGCAGTCCACATTTTCGCCGCGAGGGTTa                                                                                                        >  1:1218575/1‑37 (MQ=255)
                                                                                                             tATCGAGATGCAGTCCACATTTTCGCCGCGAGGGTTa                                                                                                        <  2:1218575/37‑1 (MQ=255)
                                                                                                              aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGtt                                                     <  1:1490780/86‑1 (MQ=255)
                                                                                                              aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGtt                                                     >  2:1490780/1‑86 (MQ=255)
                                                                                                              aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAg                                        <  1:905111/99‑1 (MQ=255)
                                                                                                              aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAg                                        >  2:905111/1‑99 (MQ=255)
                                                                                                              aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACAt                      >  1:834855/1‑117 (MQ=255)
                                                                                                              aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACAt                      <  2:834855/117‑1 (MQ=255)
                                                                                                              aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa  <  1:1279966/137‑1 (MQ=255)
                                                                                                              aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa  >  2:1279966/1‑137 (MQ=255)
                                                                                                                    aTGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttcttc                                                                                <  1:1338497/53‑1 (MQ=255)
                                                                                                                    aTGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttcttc                                                                                >  2:1338497/1‑53 (MQ=255)
                                                                                                                     tGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGAc                                          >  1:465836/1‑90 (MQ=255)
                                                                                                                     tGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGAc                                          <  2:465836/90‑1 (MQ=255)
                                                                                                                     tGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa  >  1:557584/1‑130 (MQ=255)
                                                                                                                     tGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa  <  2:557584/130‑1 (MQ=255)
                                                                                                                          |                                                                                                                             
TCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGCCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCA  >  NC_000913/4535319‑4535565

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