Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,769,478 (G)7→8 coding (405/1113 nt) ydiK → UPF0118 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,769,4711.G95.6% 147.5 / ‑0.2 45coding (398/1113 nt)ydiKUPF0118 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (23/20);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

TCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTG  >  NC_000913/1769339‑1769596
                                                                                                                                     |                                                                                                                             
tcgtcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATgggggg                                                                                                                         >  1:483708/1‑139 (MQ=255)
   tcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATgggggggga                                                                                                                      >  1:262835/1‑137 (MQ=255)
   tcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATgggggggga                                                                                                                      >  1:749584/1‑137 (MQ=255)
      aCGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCGGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGACACAACTTGCTGGATATggggggggacgg                                                                                                                   >  2:540578/1‑134 (MQ=255)
             cGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCAt                                                                                                             >  2:285115/1‑139 (MQ=255)
                         tAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCg                                                                                                 <  1:294355/139‑1 (MQ=255)
                                 tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTtata                                                                                         >  2:359884/1‑139 (MQ=255)
                                 tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTtata                                                                                         >  1:665979/1‑139 (MQ=255)
                                 tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTCGCTGGATATGGGGGGGGACGGCGATCATGGCGCAAGTCCGCCCTtata                                                                                         >  1:539120/1‑139 (MQ=255)
                                 tttCCAGCGGTGACAAGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACACCTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTtata                                                                                         >  1:533665/1‑139 (MQ=255)
                                        cGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   <  2:679596/122‑1 (MQ=255)
                                        cGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   >  1:679596/1‑122 (MQ=255)
                                           tGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   >  2:504924/1‑119 (MQ=255)
                                           tGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                   <  1:504924/119‑1 (MQ=255)
                                             aCATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccacca                                                                             <  1:640278/139‑1 (MQ=255)
                                                             tCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTgg                                                             >  2:801936/1‑139 (MQ=255)
                                                               tGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGc                                                           <  2:483708/139‑1 (MQ=255)
                                                                 gCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAg                                                         <  1:46972/139‑1 (MQ=255)
                                                                  cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                        <  1:359884/139‑1 (MQ=255)
                                                                  cGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGCCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                        >  1:533475/1‑139 (MQ=255)
                                                                     ggCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAggcgg                                                     <  1:285115/139‑1 (MQ=255)
                                                                               cATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCgg                                           <  1:290333/139‑1 (MQ=255)
                                                                                aTTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCggg                                          <  1:214173/139‑1 (MQ=255)
                                                                                aTTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCggg                                          >  2:741784/1‑139 (MQ=255)
                                                                                 ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGGACCACCACCACCTAGTTCGTGGGGCAGGCGGCGCATATCGGgc                                         >  1:779550/1‑139 (MQ=255)
                                                                                 ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTg                                                                      <  1:509252/110‑1 (MQ=255)
                                                                                 ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTg                                                                      >  2:509252/1‑110 (MQ=255)
                                                                                 ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCggg                                          <  2:756737/138‑1 (MQ=255)
                                                                                 ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCggg                                          >  1:756737/1‑138 (MQ=255)
                                                                                  tCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgcg                                        <  2:234919/139‑1 (MQ=255)
                                                                                  tCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGgcg                                        <  1:51809/139‑1 (MQ=255)
                                                                                    cGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                      <  1:766794/139‑1 (MQ=255)
                                                                                              cgcgAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATAtt                                                                                       >  1:527399/1‑80 (MQ=255)
                                                                                              cgcgAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATAtt                                                                                       <  2:527399/80‑1 (MQ=255)
                                                                                                 gAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATtg                         >  1:690498/1‑139 (MQ=255)
                                                                                                        tATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGTGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGa                 >  2:293221/1‑139 (MQ=255)
                                                                                                        tATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGa                 >  2:50813/1‑139 (MQ=255)
                                                                                                               gCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                      >  1:72545/1‑112 (MQ=255)
                                                                                                               gCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                      <  2:72545/112‑1 (MQ=255)
                                                                                                               gCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTtc           >  2:67759/1‑139 (MQ=255)
                                                                                                               gCTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGCCACCACCACCACCTGGTTCGTTGTGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTtc           <  2:419863/139‑1 (MQ=255)
                                                                                                                cTGGCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTtct          >  2:747742/1‑139 (MQ=255)
                                                                                                                  ggCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAGAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTctct        >  2:467191/1‑139 (MQ=255)
                                                                                                                  ggCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                      <  1:100282/109‑1 (MQ=255)
                                                                                                                  ggCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTctct        >  2:593553/1‑139 (MQ=255)
                                                                                                                  ggCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTctct        >  1:402710/1‑139 (MQ=255)
                                                                                                                  ggCACAACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGCTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                      >  2:100282/1‑109 (MQ=255)
                                                                                                                    cacaACTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTc      >  2:36156/1‑139 (MQ=255)
                                                                                                                        aCTTGCTGGATATGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTg  <  1:206943/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                             
TCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTG  >  NC_000913/1769339‑1769596

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