Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,769,478 (G)7→9 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.8% 75.6 / ‑3.1 24coding (398/1113 nt)ydiKUPF0118 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (15/8);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01
*NC_0009131,769,4712.G91.7% 76.5 / 0.4 24coding (398/1113 nt)ydiKUPF0118 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (14/8);  total (15/9)
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.34e-01

TATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCT  >  NC_000913/1769337‑1769600
                                                                                                                                         ||                                                                                                                                  
tATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑gggg                                                                                                                                <  2:525801/139‑1 (MQ=255)
  tcgtcgACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑gggggg                                                                                                                              >  1:648859/1‑139 (MQ=255)
                                   tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTtat                                                                                              >  2:661817/1‑139 (MQ=255)
                                   tttCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCACTGCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTtat                                                                                              >  2:1169404/1‑139 (MQ=255)
                                          cGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGca                                                                                       <  2:648859/139‑1 (MQ=255)
                                          cGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGca                                                                                       >  2:513900/1‑139 (MQ=255)
                                                            ggATCTTGCCTGGCTTAATACCATTCCGGGTATTGGCGCGAAGCTGTAAGCAAGCTGGCACAACCTGCCTGGAATGG‑‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACGACCACCTGGTTCg                                                                     <  2:881492/138‑1 (MQ=255)
                                                                        gCTTAATACTATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGGCGGGCGTTCTTGGCGAAAGTCCGCCCTTATTTTGGCACCACCACAACCTGGTTCGTTGGACAggccc                                                         >  1:881492/1‑137 (MQ=255)
                                                                                  aTTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggg                                                                                                                                 <  2:697938/56‑1 (MQ=255)
                                                                                  aTTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑ggg                                                                                                                                 >  1:697938/1‑56 (MQ=255)
                                                                                   ttCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCggg                                              >  2:1217580/1‑139 (MQ=255)
                                                                                             tGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATgg                                    >  2:936346/1‑139 (MQ=255)
                                                                                                        tGTATGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTgcg                         <  1:1217580/139‑1 (MQ=255)
                                                                                                            tGCAGGCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGa                     <  1:431406/139‑1 (MQ=255)
                                                                                                               aGGCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                       >  1:956567/1‑54 (MQ=255)
                                                                                                               aGGCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTc                                                                                                       <  2:956567/54‑1 (MQ=255)
                                                                                                                 gCTTGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCtt                <  1:847374/139‑1 (MQ=255)
                                                                                                                 gCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                          >  1:725732/1‑113 (MQ=255)
                                                                                                                 gCTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCt                                          <  2:725732/113‑1 (MQ=255)
                                                                                                                  cTGGCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTtc               >  1:1213146/1‑139 (MQ=255)
                                                                                                                    ggCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTtctc             >  2:529216/1‑139 (MQ=255)
                                                                                                                    ggCACAACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTtctc             >  1:289879/1‑139 (MQ=255)
                                                                                                                      cactACTTGCTGGATATGGGGGGGGGCCGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCTTATCGGCCGCTTAATGGTGCATTGTGCGCTGATGCTTCTCtt           >  2:847374/1‑139 (MQ=255)
                                                                                                                      cacaACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccacc                                                                                  <  1:777333/68‑1 (MQ=255)
                                                                                                                      cacaACTTGCTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGcaccacc                                                                                  >  2:777333/1‑68 (MQ=255)
                                                                                                                                cTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                            <  2:446549/81‑1 (MQ=255)
                                                                                                                                cTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAgg                                                            >  1:446549/1‑81 (MQ=255)
                                                                                                                                cTGGATATGGGGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCct  >  1:169818/1‑139 (MQ=255)
                                                                                                                                         ||                                                                                                                                  
TATCGTCGACGGCAGCGGCCCGCTAATTAAAGCCATTTCCAGCGGTGACATGACGTTACCCGATCTGGCGTGGCTTAATACCATTCCGGTGATTGGCGCGAAGCTGTATGCAGGCTGGCACAACTTGCTGGATAT‑‑GGGGGGGACGGCGATCATGGCGAAAGTCCGCCCTTATATTGGCACCACCACCACCTGGTTCGTTGGGCAGGCGGCGCATATCGGGCGCTTTATGGTGCATTGTGCGCTGATGCTTCTCTTCAGTGCCCT  >  NC_000913/1769337‑1769600

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: