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,4410CT97.7% 136.1 / ‑6.4 44pseudogene (439/1029 nt)yjhRpseudogene, helicase family, putative frameshift suppressor
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (21/22);  total (22/22)
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

TGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGCCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGGTGTGGA  >  NC_000913/4535317‑4535574
                                                                                                                            |                                                                                                                                     
tGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTcgccg                                                                                                                         <  2:850332/139‑1 (MQ=255)
  tCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCcgcg                                                                                                                       >  1:938827/1‑139 (MQ=255)
  tCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCcgcg                                                                                                                       >  1:561210/1‑139 (MQ=255)
           aaCATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTAtt                                                                                                              >  1:42916/1‑139 (MQ=255)
             cATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTgg                                                                                                            <  1:211210/139‑1 (MQ=255)
             cATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTgg                                                                                                            <  2:42916/139‑1 (MQ=255)
                gAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCaa                                                                                                         >  1:434278/1‑139 (MQ=255)
                       gCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATAtct                                                                                                  >  1:579105/1‑139 (MQ=255)
                            tttCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtttt                                                                                             <  2:938827/139‑1 (MQ=255)
                            tttCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtttt                                                                                             <  2:561210/139‑1 (MQ=255)
                              tCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttct                                                                                           <  2:1056852/139‑1 (MQ=255)
                                                     tCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAAtt                                                                    >  1:17194/1‑139 (MQ=255)
                                                            tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCa                                                                                                          <  2:94981/94‑1 (MQ=255)
                                                            tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCa                                                                                                          >  1:94981/1‑94 (MQ=255)
                                                            tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTGTTTTCTTCAGATAATAATGCGTTGCAATTTGAGttt                                                             >  1:699720/1‑139 (MQ=255)
                                                            tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttct                                                                                           >  1:840300/1‑109 (MQ=255)
                                                            tGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttct                                                                                           <  2:840300/109‑1 (MQ=255)
                                                             gCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTc                                                            >  1:750547/1‑139 (MQ=255)
                                                             gCTGTCTCACGCACTAAAGATAGTTTTCTGGCATTTGGCGATATGGACCTTATCGAGCTGCAGCCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTc                                                            >  1:221434/1‑139 (MQ=255)
                                                              cTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCa                                                           >  1:371459/1‑139 (MQ=255)
                                                                 tctcACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttct                                                                                           <  1:256510/104‑1 (MQ=255)
                                                                 tctcACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTttct                                                                                           >  2:256510/1‑104 (MQ=255)
                                                                     aCGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtt                                                                                               >  1:895845/1‑96 (MQ=255)
                                                                     aCGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtt                                                                                               <  2:895845/96‑1 (MQ=255)
                                                                     aCGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGa                                                    >  1:338791/1‑139 (MQ=255)
                                                                     aCGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGa                                                    <  1:562264/139‑1 (MQ=255)
                                                                     aCGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGa                                                    <  2:434278/139‑1 (MQ=255)
                                                                      cgcgCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGAc                                                   <  1:1103401/139‑1 (MQ=255)
                                                                             aaGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGAttt                                            >  2:267077/1‑139 (MQ=255)
                                                                                     tttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCg                                    <  2:699720/139‑1 (MQ=255)
                                                                                     tttCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCg                                    >  1:306478/1‑139 (MQ=255)
                                                                                       tCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGca                                  >  2:712780/1‑139 (MQ=255)
                                                                                                            ccTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAg                                                 <  2:78279/103‑1 (MQ=255)
                                                                                                            ccTTATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAg                                                 >  1:78279/1‑103 (MQ=255)
                                                                                                              ttATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa           <  2:446040/139‑1 (MQ=255)
                                                                                                              ttATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa           <  2:221434/139‑1 (MQ=255)
                                                                                                               tATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTgg                                                                                                            <  2:841109/41‑1 (MQ=255)
                                                                                                               tATCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTgg                                                                                                            >  1:841109/1‑41 (MQ=255)
                                                                                                                aTCGAGATGTAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa           <  2:535044/137‑1 (MQ=255)
                                                                                                                aTCGAGATGTAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCa           >  1:535044/1‑137 (MQ=255)
                                                                                                                aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtttt                                                                                             <  1:889306/55‑1 (MQ=255)
                                                                                                                aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCtttt                                                                                             >  2:889306/1‑55 (MQ=255)
                                                                                                                aTCGAGATGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATg         <  1:712780/139‑1 (MQ=255)
                                                                                                                       tGCAGTCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGGTGTGGa  <  2:579105/139‑1 (MQ=255)
                                                                                                                            |                                                                                                                                     
TGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTGGCGATATGGACCTTATCGAGATGCAGCCCGCATTTTCGCCGCGAGGGTTATTGGCAAAATATCTCTTTTCTTCAGATAATAATGCGTTGCAATTTGAGTTTCAGAAACGACAGGATTTAATTTCCGCACATACACAAATTTCAACGCTGCATGGTGTGGA  >  NC_000913/4535317‑4535574

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