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.C98.0% 176.2 / ‑3.4 51coding (500/732 nt)yfiHUPF0124 family protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (28/22);  total (29/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 = 9.81e-01

AAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGC  >  NC_000913/2734403‑2734662
                                                                                                                                      |                                                                                                                                
aaaaTTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑cccccc                                                                                                                           <  1:635463/139‑1 (MQ=255)
                                  gCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGaa                                                                                                                 >  2:431696/1‑116 (MQ=255)
                                  gCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGaa                                                                                                                 <  1:431696/116‑1 (MQ=255)
                                       gCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                                                     <  1:1341848/139‑1 (MQ=255)
                                                tAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                          >  2:67660/1‑125 (MQ=255)
                                                tAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                          <  1:67660/125‑1 (MQ=255)
                                                         gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCt                                                                                                                                   <  2:776174/76‑1 (MQ=255)
                                                         gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCt                                                                                                                                   >  1:776174/1‑76 (MQ=255)
                                                             gATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGa                                                                <  1:415044/139‑1 (MQ=255)
                                                                    aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCa                                                                                   <  2:509780/113‑1 (MQ=255)
                                                                    aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCa                                                                                   >  1:509780/1‑113 (MQ=255)
                                                                    aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCCACTTAGAACGCGCGTGGACCAATTGCCGGCCCAAACCAGGCGAGAATATTTCCCGGATTAACAg                                                         >  2:196094/1‑139 (MQ=255)
                                                                    aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑cccccc                                                                                                                           >  2:220079/1‑71 (MQ=255)
                                                                    aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑cccccc                                                                                                                           <  1:220079/71‑1 (MQ=255)
                                                                       aCCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaaa                                                     >  1:873603/1‑139 (MQ=255)
                                                                        ccATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTCCGACCCGGCGGGACCCAATGGCCGCCCTTACCAAGCCAAGATTATTTTCCGGATAGCCAGCaaaa                                                    >  2:843560/1‑139 (MQ=255)
                                                                         cATGCTGCATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAAc                                                   >  2:938079/1‑139 (MQ=255)
                                                                         cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAcgcg                                                                                                             <  2:1347728/81‑1 (MQ=255)
                                                                         cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAAcgcg                                                                                                             >  1:1347728/1‑81 (MQ=255)
                                                                         cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                                                 >  1:477193/1‑109 (MQ=255)
                                                                         cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                                                 <  2:477193/109‑1 (MQ=255)
                                                                           tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCgg                                                                                            <  1:1198096/96‑1 (MQ=255)
                                                                           tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCgg                                                                                            >  2:1198096/1‑96 (MQ=255)
                                                                           tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                                                 >  1:371812/1‑139 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCt                                                                                                                                   <  2:584287/57‑1 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCt                                                                                                                                   >  1:584287/1‑57 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCCCTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                                >  1:1075758/1‑139 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTaa                                                                                      <  1:16371/101‑1 (MQ=255)
                                                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTaa                                                                                      >  2:16371/1‑101 (MQ=255)
                                                                                aatgaaAGCTGCACTTGCTTTAGCGTCTAGTGCCATAAACGCCTCCCGAACCAC‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAAAAGGAAAc                                            >  2:827241/1‑139 (MQ=255)
                                                                                   gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               >  1:284762/1‑85 (MQ=255)
                                                                                   gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               <  2:284762/85‑1 (MQ=255)
                                                                                      aGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtctc                                      >  2:182952/1‑139 (MQ=255)
                                                                                          gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                                                     <  2:574275/88‑1 (MQ=255)
                                                                                          gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                                                     >  1:574275/1‑88 (MQ=255)
                                                                                          gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTcc                                  >  2:88328/1‑139 (MQ=255)
                                                                                                  tttAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGcc                          <  1:182952/139‑1 (MQ=255)
                                                                                                   ttAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCg                                          <  1:87342/122‑1 (MQ=255)
                                                                                                   ttAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCg                                          >  2:87342/1‑122 (MQ=255)
                                                                                                        gTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               >  2:999865/1‑64 (MQ=255)
                                                                                                        gTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               <  1:999865/64‑1 (MQ=255)
                                                                                                        gTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca                    <  1:938079/139‑1 (MQ=255)
                                                                                                          cTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                              <  2:395159/127‑1 (MQ=255)
                                                                                                          cTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCa                              >  1:395159/1‑127 (MQ=255)
                                                                                                          cTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGcaca                  >  2:1288624/1‑139 (MQ=255)
                                                                                                           tACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               <  2:500406/61‑1 (MQ=255)
                                                                                                           tACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               >  1:500406/1‑61 (MQ=255)
                                                                                                           tACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                         >  2:1379177/1‑99 (MQ=255)
                                                                                                           tACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                                                         <  1:1379177/99‑1 (MQ=255)
                                                                                                                 cATAAACGCCTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACg           <  1:88328/139‑1 (MQ=255)
                                                                                                                          cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               <  2:15319/46‑1 (MQ=255)
                                                                                                                          cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                               >  1:15319/1‑46 (MQ=255)
                                                                                                                          cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCga                                                                              <  1:840537/63‑1 (MQ=255)
                                                                                                                          cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCga                                                                              >  2:840537/1‑63 (MQ=255)
                                                                                                                          cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCa                                                       >  2:754372/1‑86 (MQ=255)
                                                                                                                          cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCa                                                       <  1:754372/86‑1 (MQ=255)
                                                                                                                          cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCcagccagc  >  1:852491/1‑139 (MQ=255)
                                                                                                                          cTCGCGAACCTCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCcagccagc  >  1:1041291/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                
AAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACGCCAGCCAGC  >  NC_000913/2734403‑2734662

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