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.C91.9% 118.6 / 2.9 37coding (500/732 nt)yfiHUPF0124 family protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (20/14);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGC  >  NC_000913/2734405‑2734662
                                                                                                                                    |                                                                                                                                
aaTTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑cccccccc                                                                                                                         <  1:931029/139‑1 (MQ=255)
         aaCACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGaa                                                                                                                 <  2:316662/139‑1 (MQ=255)
                           aCGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               <  2:1449872/139‑1 (MQ=255)
                           aCGCTGCCGGGCAAGCAGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCGC‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               <  2:1554310/139‑1 (MQ=255)
                                gCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                          >  1:1587324/1‑139 (MQ=255)
                                              tAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgaga                                                                            <  2:809710/139‑1 (MQ=255)
                                                    tCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               >  2:489680/1‑114 (MQ=255)
                                                    tCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               <  1:489680/114‑1 (MQ=255)
                                                                     aCCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaa                                                      <  2:1190786/139‑1 (MQ=255)
                                                                       cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGaa                                                                                                                 <  1:504711/77‑1 (MQ=255)
                                                                       cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGaa                                                                                                                 >  2:504711/1‑77 (MQ=255)
                                                                       cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                          <  1:1525939/100‑1 (MQ=255)
                                                                       cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                          >  2:1525939/1‑100 (MQ=255)
                                                                       cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAAc                                                   >  1:980822/1‑139 (MQ=255)
                                                                          gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAATACAgg                                                >  1:744146/1‑139 (MQ=255)
                                                                          gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATAAGCAAAACAgg                                                >  1:604977/1‑139 (MQ=255)
                                                                             gaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa                                             >  2:359784/1‑139 (MQ=255)
                                                                                     gCTGCACTTGCTTTCGCGTCTACTGCCATAAACGCCTCGCCAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTctct                                     >  2:979605/1‑139 (MQ=255)
                                                                                        gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGGGAATATTTTCCGGATTATCAGCAAAGCAGGAAACCGTCTCTTcc                                  >  2:1082079/1‑139 (MQ=255)
                                                                                          aCTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               >  1:501872/1‑76 (MQ=255)
                                                                                          aCTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               <  2:501872/76‑1 (MQ=255)
                                                                                              gCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACg                            >  1:1286900/1‑139 (MQ=255)
                                                                                                   aGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                              >  2:129168/1‑132 (MQ=255)
                                                                                                   aGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                              <  1:129168/132‑1 (MQ=255)
                                                                                                      gTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                    <  2:1527991/139‑1 (MQ=255)
                                                                                                         tACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAg                 >  2:687309/1‑139 (MQ=255)
                                                                                                         tACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAg                 >  2:1022199/1‑139 (MQ=255)
                                                                                                         tACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                         <  2:836992/98‑1 (MQ=255)
                                                                                                         tACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                         >  1:836992/1‑98 (MQ=255)
                                                                                                          aCTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                          <  2:152623/65‑1 (MQ=255)
                                                                                                          aCTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                          >  1:152623/1‑65 (MQ=255)
                                                                                                                   aaCGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACgccag       <  2:980822/139‑1 (MQ=255)
                                                                                                                        cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               <  2:1187140/46‑1 (MQ=255)
                                                                                                                        cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               >  1:1187140/1‑46 (MQ=255)
                                                                                                                        cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCcagccagc  >  1:43531/1‑139 (MQ=255)
                                                                                                                        cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGACCACCCcagcccgc  >  1:408648/1‑139 (MQ=255)
                                                                                                                          cgcgAACCTCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTg                                                                                                >  1:988858/1‑44 (MQ=255)
                                                                                                                          cgcgAACCTCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTg                                                                                                <  2:988858/44‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                
AATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGC  >  NC_000913/2734405‑2734662

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