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.C96.7% 208.7 / ‑0.5 60coding (500/732 nt)yfiHUPF0124 family protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (27/31);  total (28/32)
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.95e-01

GCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACG  >  NC_000913/2734399‑2734653
                                                                                                                                         |                                                                                                                       
gccgAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑cc                                                                                                                      <  2:1525801/139‑1 (MQ=255)
          tGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑cccccccccact                                                                                                            <  2:1257620/139‑5 (MQ=255)
               aaCACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGaa                                                                                                        <  1:220363/139‑1 (MQ=255)
                      aCGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑cccccc                                                                                                                  >  1:1669344/1‑121 (MQ=255)
                      aCGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑cccccc                                                                                                                  <  2:1669344/121‑1 (MQ=255)
                          tCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGAc                                                                                             <  2:124017/139‑1 (MQ=255)
                           cGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGAcc                                                                                            <  2:83422/139‑1 (MQ=255)
                           cGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGAcc                                                                                            <  2:1388387/139‑1 (MQ=255)
                                               gCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                                        <  2:794940/139‑1 (MQ=255)
                                                cTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGAc                                                                                             <  1:491912/117‑1 (MQ=255)
                                                cTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGAc                                                                                             >  2:491912/1‑117 (MQ=255)
                                                 tGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCCACt                                                                                                             <  2:1664729/99‑1 (MQ=255)
                                                 tGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCCACt                                                                                                             >  1:1664729/1‑99 (MQ=255)
                                                         atCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      <  2:1068110/115‑1 (MQ=255)
                                                         atCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      >  1:1068110/1‑115 (MQ=255)
                                                           ccGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                 <  2:153934/118‑1 (MQ=255)
                                                           ccGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                                                 >  1:153934/1‑118 (MQ=255)
                                                             gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      >  1:1174948/1‑111 (MQ=255)
                                                             gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      <  2:1174948/111‑1 (MQ=255)
                                                                  aTACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGAtt                                                     <  2:1143675/139‑1 (MQ=255)
                                                                        aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      >  2:726478/1‑100 (MQ=255)
                                                                        aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      <  1:726478/100‑1 (MQ=255)
                                                                        aTCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGc                                               >  1:284348/1‑139 (MQ=255)
                                                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      >  1:108049/1‑93 (MQ=255)
                                                                               tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      <  2:108049/93‑1 (MQ=255)
                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGaa                                                                                                        >  2:709459/1‑74 (MQ=255)
                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGaa                                                                                                        <  1:709459/74‑1 (MQ=255)
                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      >  2:1671364/1‑92 (MQ=255)
                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      >  2:723241/1‑92 (MQ=255)
                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      <  2:1128487/92‑1 (MQ=255)
                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      <  1:723241/92‑1 (MQ=255)
                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      <  1:1671364/92‑1 (MQ=255)
                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      >  1:1128487/1‑92 (MQ=255)
                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                       >  1:1386900/1‑139 (MQ=255)
                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCCACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                                        >  1:1500747/1‑106 (MQ=255)
                                                                                gCTGAATGAAAGCTGCACTTGCTTTAGCGTCCACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                                        <  2:1500747/106‑1 (MQ=255)
                                                                                gCTGAATGAAACCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                                       <  2:284348/139‑1 (MQ=255)
                                                                                    aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAAc                                   >  2:1281956/1‑139 (MQ=255)
                                                                                              gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTc                          >  1:1479476/1‑139 (MQ=255)
                                                                                              gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGc                                                                                  >  1:1441832/1‑82 (MQ=255)
                                                                                              gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGc                                                                                  <  2:1441832/82‑1 (MQ=255)
                                                                                              gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTcc                         >  1:1620606/1‑139 (MQ=255)
                                                                                                  ttGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa                                    <  1:907010/124‑1 (MQ=255)
                                                                                                  ttGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa                                    >  2:907010/1‑124 (MQ=255)
                                                                                                            gTCTACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca           <  1:1281956/139‑1 (MQ=255)
                                                                                                             tCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca           <  2:1479476/139‑1 (MQ=255)
                                                                                                               tACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      >  1:1151341/1‑61 (MQ=255)
                                                                                                               tACTGCCATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      <  2:1151341/61‑1 (MQ=255)
                                                                                                                    ccATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      <  1:1646975/56‑1 (MQ=255)
                                                                                                                    ccATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      >  2:1646975/1‑56 (MQ=255)
                                                                                                                    ccATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTc                                                           >  2:144818/1‑83 (MQ=255)
                                                                                                                    ccATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTc                                                           <  1:144818/83‑1 (MQ=255)
                                                                                                                     cATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca           >  1:559848/1‑130 (MQ=255)
                                                                                                                     cATAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGca           <  2:559848/130‑1 (MQ=255)
                                                                                                                      aTAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAat                                                                 >  1:378704/1‑75 (MQ=255)
                                                                                                                      aTAAACGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAat                                                                 <  2:378704/75‑1 (MQ=255)
                                                                                                                        aaaCGCCTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTAt                                                   <  1:1540507/87‑1 (MQ=255)
                                                                                                                        aaaCGCCTAGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTCTCCGGATTAt                                                   >  2:1540507/1‑87 (MQ=255)
                                                                                                                              cTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      >  2:886813/1‑46 (MQ=255)
                                                                                                                              cTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      <  1:886813/46‑1 (MQ=255)
                                                                                                                              cTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACg  <  1:1322054/130‑1 (MQ=255)
                                                                                                                              cTCGCGAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACg  >  2:1322054/1‑130 (MQ=255)
                                                                                                                                  cgAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      >  1:93200/1‑42 (MQ=255)
                                                                                                                                  cgAACCTC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGc                                                                                      <  2:93200/42‑1 (MQ=255)
                                                                                                                                         |                                                                                                                       
GCCGAAAATTTGCTCAACACCCACGTTCGCCAGACGCTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCAGCACGCCTGCGCACAGTCCACG  >  NC_000913/2734399‑2734653

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