Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,113,668 C→T G30G (GGC→GGT yceK → DUF1375 domain‑containing lipoprotein YceK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,113,6680CT96.5% 82.8 / ‑6.0 29G30G (GGC→GGTyceKDUF1375 domain‑containing lipoprotein YceK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/14);  minor base G (1/0);  total (15/14)
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.78e-01

GATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTTCATCCCACTCATCAGCTGCTGCGCGTCCTTCCTCGGTATCCAGCGGTGGCTCAAGCTGAAATTCCCCCTCGTCCCATTCATGTAATGTATTC  >  NC_000913/1113421‑1113906
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
gATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGCGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAAc                                                                                                                                                                                                                                             >  2:424766/1‑251 (MQ=255)
             tATCAGTCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATccc                                                                                                                                                                                                                                 <  1:266152/250‑1 (MQ=255)
                                     aTTTATGGATTTCTCAAGTCTGTGTGTGTCGTATAAGGCGTCCACGCCGCACCCGGCAAAAACAACTAGCCAATAAAAATACCGTCGTTATGCCGGTATTTTTTTACGAAAGATTTATCAAATGCGTTTAATTGTGGTGAGCATCATTGTGACCTTGCTGAGTGGCTTTGGCAGCATTATTAGTCGCACTAGACCGGGGCAGGTGTTTGGTAACCAATCTTATCCCGTTGTGCAATGGTTTGTGCGTGAca                                                                                                                                                                                                        <  2:34503/251‑2 (MQ=255)
                                                          ttGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGCGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGACGCACTATACCGGGGCAGGGGCATGGTAACAAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACg                                                                                                                                                                                    >  2:160913/1‑250 (MQ=255)
                                                                    aaTAAGGCTTACACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATTCCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCt                                                                                                                                                                         <  2:29948/250‑1 (MQ=255)
                                                                                 cgccgcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTg                                                                                                                                                            <  1:269417/251‑1 (MQ=255)
                                                                                           ggCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGAtt                                                                                                                                                  >  2:128530/1‑250 (MQ=255)
                                                                                                 aacaacGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTAGCACTATACAGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGACTGGCGTTATGTAACGATCCTTGATATGCCATTCTCTCTGGTTTTTGATCCTTTAc                                                                                                                                            >  2:320558/1‑251 (MQ=255)
                                                                                                       gAGCCACTCAAACTACCGGCGCTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTAGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCCGGGGCATGGTAACCAATATTATCACGGTGTGCAATGGGCTGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATtct                                                                                                                                                                 >  2:513989/1‑224 (MQ=255)
                                                                                                       gAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGTGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATtct                                                                                                                                                                 <  1:513989/224‑1 (MQ=255)
                                                                                                                aaaaTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAAGTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACa                                                                                                                             <  1:412751/251‑1 (MQ=255)
                                                                                                                                gCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGt                                                                                                              >  1:122028/1‑250 (MQ=255)
                                                                                                                                            ttACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGAtg                                                                                                                                                                                                                 <  2:144519/139‑1 (MQ=255)
                                                                                                                                            ttACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGAtg                                                                                                                                                                                                                 >  1:144519/1‑139 (MQ=255)
                                                                                                                                                         aTCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGGGACTCCGCCTGGCGTTATGTCACGTTCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCGGATCGACATTCATCATGGCCCGAATGAGGGATTAACGTTCATCCCACtc                                                                                    >  1:423414/1‑251 (MQ=255)
                                                                                                                                                                                     catGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTCTACCGGGGCAGGGGCATGGGAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGCTACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTTCATCCCACTCCTCAGCTGCTGCGCGTCCTTCCTCGGt                                                         >  2:358948/1‑250 (MQ=255)
                                                                                                                                                                                           gACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTTCATCCCACTCATCAGCTGCTGCGCGTCCTTCCTCGGTATCCAGc                                                  <  1:128530/251‑1 (MQ=255)
                                                                                                                                                                                             ccTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGGCACGATCCTTGATCTGACATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGCTCGACATTCATCATGGCCCGTATGAGTGATTAACGTTCATCCCACTCATCAGCTGCTGCGCGTCCTTCCTCGGTATCCAGCgg                                                >  2:459512/1‑251 (MQ=255)
                                                                                                                                                                                                    gAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTTCATCCCACTCATCCGCTGCTGCGCGTCCTTCCTCGGTATCCAGCGGTGGCTc                                          >  2:47941/1‑250 (MQ=255)
                                                                                                                                                                                                     aGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTAGCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCTTTATGTCACGATCCTTGATCTGc                                                                                                                                                                       <  2:223596/124‑1 (MQ=255)
                                                                                                                                                                                                     aGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGTAAACAATATTATCCCGGTGTGCAATGGGCTGTGCGTGACTACGCCTGGCGTTATGTCACGATCCTTGATCTGc                                                                                                                                                                       >  1:223596/1‑124 (MQ=255)
                                                                                                                                                                                                           tgtgGCAGCATTATTAGACGCACTATACCGGGGCAGGGGCCTGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTTCATCCCACTCATCAGCTGCTGCGCGTCCTTCCTCGGTATCCAGCGGTGACTCCAGCTgt                                   >  2:450154/1‑249 (MQ=255)
                                                                                                                                                                                                                       attaGTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGAGTAACGTTCATCCCACTCATCAGCTGCTGCGCGTCCTTCCTCGGTATCCAGCGGTGGCTCAAGCTGAAATTCCCCCTCGt                      <  2:122028/251‑1 (MQ=255)
                                                                                                                                                                                                                           gTCGCACTATACCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTTCATCCCACTCATCAGCTGCTGCGCGTCCTTCCTCGGTATCCAGCGGTGGCTCAAGCTGAAATTCCCCCTCGTccc                   >  1:69888/1‑250 (MQ=255)
                                                                                                                                                                                                                           gTCGCACCATTCCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTTCATCCCACTCATCAGCTGCTGCGCGTCCTTCCTCGGTATCCAGCGGTGGCTCAAGCTGAAATTCCCCCTCGTCCCa                  <  2:423414/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ataCCGGGGCAGGGGCATGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTTCATCCCACTCATCAGCTGCTGCGCGTCCGTCCTCGGTATCCAGCTGTGGCTCAAGCTGAAATTCCCCCTCGTCCCATTCatgta          <  1:320558/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ggggcaggggcaTGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTGTACTGCTGCCGATCGACATTCAGCATGGCCCGTATGAGTGATTAACGTTCATCCCCCTCATCAGCTGCTGCTCGTCCTTCCTCGGTAGCCAGCGGTGGCTCAAGCTTAAATTCCCCCTCGTCCCATTCatgtaatgta     <  2:69888/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gggcaggggcaTGGTAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTTCATCCCACTCATCAGCTGCTGCGCGTCCTTCCTCGGTATCCAGCGGTGGCTCAAGCTGAAATTCCCCCTCGTCCCATTCATGTAATGTAt    <  1:358948/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gcaggggcaAGGTCACAAATATTATCCAGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTTCATCCCACTCCTCAGCTGCTGCGCGTCCTTCCTCGGTATCCAGCGGTGGCTCAAGCTTAAATTCCCCCTCGTCCCATTCatgtactgttttt  >  2:306787/1‑247 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
GATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGGATGTGCGTGACTCCGCCTGGCGTTATGTCACGATCCTTGATCTGCCATTCTCTCTGGTTTTTGATACTTTACTGCTGCCGATCGACATTCATCATGGCCCGTATGAGTGATTAACGTTCATCCCACTCATCAGCTGCTGCGCGTCCTTCCTCGGTATCCAGCGGTGGCTCAAGCTGAAATTCCCCCTCGTCCCATTCATGTAATGTATTC  >  NC_000913/1113421‑1113906

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