Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,734,543 (C)8→9 coding (492/732 nt) yfiH ← UPF0124 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,734,5351.C93.8% 155.7 / 2.8 48coding (500/732 nt)yfiHUPF0124 family protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (24/21);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

AAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAG  >  NC_000913/2734403‑2734657
                                                                                                                                      |                                                                                                                          
aaaaTTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGGTAAATTTCCGCCAGATACTTAACACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcccccc                                                                                                                      <  2:1576608/139‑1 (MQ=255)
                             aCGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                          >  1:397151/1‑139 (MQ=255)
                              cGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGcc                                                                                         >  1:436197/1‑139 (MQ=255)
                                             tGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGaa                                                                                                            >  2:446212/1‑105 (MQ=255)
                                             tGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGaa                                                                                                            <  1:446212/105‑1 (MQ=255)
                                                      tCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAc                                                                                                           >  2:930541/1‑97 (MQ=255)
                                                      tCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGGCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAc                                                                                                           <  1:930541/97‑1 (MQ=255)
                                                         gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCt                                                                                                                             <  2:1083269/76‑1 (MQ=255)
                                                         gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCt                                                                                                                             >  1:1083269/1‑76 (MQ=255)
                                                         gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTcc                                                              >  2:863554/1‑139 (MQ=255)
                                                              aTACTTATCACCATGCTGAATTAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGAtt                                                         <  2:314880/139‑1 (MQ=255)
                                                              aTACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGAtt                                                         <  2:436197/139‑1 (MQ=255)
                                                              aTACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACACCTCGCGAACCTCCCCCCCCCACTTCGaa                                                                                                            >  2:703254/1‑88 (MQ=255)
                                                              aTACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACACCTCGCGAACCTCCCCCCCCCACTTCGaa                                                                                                            <  1:703254/88‑1 (MQ=255)
                                                                         cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAc                                                                                                           >  2:646398/1‑78 (MQ=255)
                                                                         cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAc                                                                                                           <  1:646398/78‑1 (MQ=255)
                                                                           tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                                            >  1:1462350/1‑139 (MQ=255)
                                                                           tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                                            >  2:946322/1‑139 (MQ=255)
                                                                           tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGt                                                                                                     <  1:1337173/81‑1 (MQ=255)
                                                                           tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGt                                                                                                     >  2:1337173/1‑81 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAc                                                                                                           >  2:1238678/1‑75 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAc                                                                                                           <  1:1238678/75‑1 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgaga                                                                       >  1:891023/1‑111 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgaga                                                                       <  2:891023/111‑1 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                           >  2:872308/1‑139 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                           >  2:1125350/1‑139 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                           >  2:415606/1‑139 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                           >  1:267354/1‑139 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcc                                                                                                                          >  1:170395/1‑59 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcc                                                                                                                          <  2:170395/59‑1 (MQ=255)
                                                                              tgaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa                                         >  2:1156271/1‑139 (MQ=255)
                                                                                aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                    >  1:255791/1‑126 (MQ=255)
                                                                                aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                    <  2:255791/126‑1 (MQ=255)
                                                                                  tgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcc                                                                                                                          <  1:1082760/53‑1 (MQ=255)
                                                                                  tgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcc                                                                                                                          >  2:1082760/1‑53 (MQ=255)
                                                                                   gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGt                                    <  2:378564/139‑1 (MQ=255)
                                                                                   gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGt                                    >  2:523587/1‑139 (MQ=255)
                                                                                          gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTcc                             >  2:1539653/1‑139 (MQ=255)
                                                                                              ttGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                         <  1:80980/139‑1 (MQ=255)
                                                                                              ttGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                         <  2:397151/139‑1 (MQ=255)
                                                                                              ttGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGAGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                         <  1:1156271/139‑1 (MQ=255)
                                                                                                     aGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                     <  1:568651/72‑1 (MQ=255)
                                                                                                     aGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                     >  2:568651/1‑72 (MQ=255)
                                                                                                        gTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca               <  2:22788/139‑1 (MQ=255)
                                                                                                        gTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca               <  1:523587/139‑1 (MQ=255)
                                                                                                           tACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                          >  1:1234629/1‑61 (MQ=255)
                                                                                                           tACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                          <  2:1234629/61‑1 (MQ=255)
                                                                                                                  aTAAACGCCTCGCGAACCCCCCCCCCCCACTTCGAACGCGCGTGGACCAAGTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTAGCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgc     <  1:844344/139‑1 (MQ=255)
                                                                                                                   tAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca               >  1:1058572/1‑128 (MQ=255)
                                                                                                                   tAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca               <  2:1058572/128‑1 (MQ=255)
                                                                                                                     aaCGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgccag  <  2:267354/139‑1 (MQ=255)
                                                                                                                          cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgc     >  2:584830/1‑131 (MQ=255)
                                                                                                                          cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgc     <  1:584830/131‑1 (MQ=255)
                                                                                                                            cgcgAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                    >  1:515414/1‑82 (MQ=255)
                                                                                                                            cgcgAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                    <  2:515414/82‑1 (MQ=255)
                                                                                                                                      |                                                                                                                          
AAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAG  >  NC_000913/2734403‑2734657

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: