Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,497,308 A→G E151E (GAA→GAG opgD → OPG biosynthetic periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,497,3080AG98.2% 188.7 / ‑7.1 55E151E (GAA→GAGopgDOPG biosynthetic periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (31/23);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCC  >  NC_000913/1497176‑1497443
                                                                                                                                    |                                                                                                                                       
tatGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCaa                                                                                                                                   >  2:645526/1‑139 (MQ=255)
tatGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGACaa                                                                                                                                   >  2:385709/1‑139 (MQ=255)
          ggaTTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCg                                                                                                                         >  1:87447/1‑139 (MQ=255)
                 gTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGc                                                                                                                  <  1:779656/139‑1 (MQ=255)
                            cGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTcgcg                                                                                                       <  2:697796/139‑1 (MQ=255)
                                tGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGtt                                                                                                   >  2:603721/1‑139 (MQ=255)
                                     tCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAg                                                                                              >  1:574062/1‑139 (MQ=255)
                                            aTCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTcgc                                                                                                        >  2:906643/1‑122 (MQ=255)
                                            aTCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTcgc                                                                                                        <  1:906643/122‑1 (MQ=255)
                                                 gCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTgg                                                                                  >  1:351041/1‑139 (MQ=255)
                                                 gCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTgg                                                                                  >  2:808549/1‑139 (MQ=255)
                                                   aaCACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                                >  1:463425/1‑139 (MQ=255)
                                                   aaCACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                                <  1:808549/139‑1 (MQ=255)
                                                   aaCACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                                >  1:773707/1‑139 (MQ=255)
                                                       caTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGcc                                                                            <  2:351041/139‑1 (MQ=255)
                                                        aTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCg                                                                           >  2:692722/1‑139 (MQ=255)
                                                        aTCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCg                                                                           >  2:528567/1‑139 (MQ=255)
                                                         tCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGACAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGt                                                                          <  1:385709/139‑1 (MQ=255)
                                                            ggCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCg                                                                           >  1:860955/1‑135 (MQ=255)
                                                            ggCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCg                                                                           <  2:860955/135‑1 (MQ=255)
                                                              cgcgTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATgt                                                                     <  2:773707/139‑1 (MQ=255)
                                                               gcgTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCaaa                                                                                                                                  <  1:828584/77‑1 (MQ=255)
                                                               gcgTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCaaa                                                                                                                                  >  2:828584/1‑77 (MQ=255)
                                                                  tGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGGGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATgtagta                                                                 <  1:645526/139‑1 (MQ=255)
                                                                  tGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATgtagta                                                                 >  2:414561/1‑139 (MQ=255)
                                                                  tGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATgtag                                                                   <  1:947151/137‑1 (MQ=255)
                                                                  tGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATgtag                                                                   >  2:947151/1‑137 (MQ=255)
                                                                               cGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTTATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGgcg                                                    <  2:634223/139‑1 (MQ=255)
                                                                               cGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGgcg                                                    <  2:574062/139‑1 (MQ=255)
                                                                                         tgttCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTAttt                                          <  1:603721/139‑1 (MQ=255)
                                                                                                aTACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGcgcg                                                  >  1:119174/1‑124 (MQ=255)
                                                                                                aTACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGcgcg                                                  <  2:119174/124‑1 (MQ=255)
                                                                                                aTACAACGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGcc                                   >  1:246331/1‑139 (MQ=255)
                                                                                                     aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTtga                              >  2:204355/1‑139 (MQ=255)
                                                                                                     aCGATGCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTtga                              <  1:414561/139‑1 (MQ=255)
                                                                                                          gCAGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATa                         <  2:928491/139‑1 (MQ=255)
                                                                                                            aGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGATTAAAGCCCCCGAACTGGCGCGCCGTGATgt                                                                     >  2:996912/1‑93 (MQ=255)
                                                                                                            aGGTGTTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGATTAAAGCCCCCGAACTGGCGCGCCGTGATgt                                                                     <  1:996912/93‑1 (MQ=255)
                                                                                                              gtgtTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTcgc                                                                                                        <  1:57879/56‑1 (MQ=255)
                                                                                                              gtgtTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTcgc                                                                                                        >  2:959429/1‑56 (MQ=255)
                                                                                                              gtgtTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTcgc                                                                                                        <  1:959429/56‑1 (MQ=255)
                                                                                                              gtgtTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTcgc                                                                                                        >  2:57879/1‑56 (MQ=255)
                                                                                                              gtgtTGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACata                     >  2:62105/1‑139 (MQ=255)
                                                                                                                 ttGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCg                                                                           >  2:46075/1‑82 (MQ=255)
                                                                                                                 ttGATACAAAACAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCg                                                                           <  1:46075/82‑1 (MQ=255)
                                                                                                                        aaaaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCcgcg                                     <  1:551172/113‑1 (MQ=255)
                                                                                                                        aaaaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCcgcg                                     >  2:551172/1‑113 (MQ=255)
                                                                                                                         aaaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGtt          >  2:928155/1‑139 (MQ=255)
                                                                                                                         aaaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTGTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGtt          >  1:346529/1‑139 (MQ=255)
                                                                                                                          aaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                                >  1:261066/1‑68 (MQ=255)
                                                                                                                          aaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGgcg                                                                                <  2:261066/68‑1 (MQ=255)
                                                                                                                          aaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGttt         >  2:769224/1‑139 (MQ=255)
                                                                                                                          aaCAATTAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGttt         >  2:181647/1‑139 (MQ=255)
                                                                                                                               ttAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTtga                              <  1:827580/113‑1 (MQ=255)
                                                                                                                               ttAGAGGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTtga                              >  2:827580/1‑113 (MQ=255)
                                                                                                                                 agagGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGcc  <  1:769224/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
TATGGGAATGGGATTCCGTCGCCGCGTTCGTATGTTTTCTGTAGATCCAGCAACACATCTGGCGCGTGAAATTCACTTTCGCCCGGAGTTGTTCAAATACAACGATGCAGGTGTTGATACAAAACAATTAGAAGGGCAAAGCGATCTCGGCTTTGCCGGTTTTCGCGTGTTTAAAGCCCCCGAACTGGCGCGCCGTGATGTAGTATCATTTCTCGGCGCGAGTTATTTCCGCGCCGTTGATGATACATATCAATACGGTTTGTCGGCC  >  NC_000913/1497176‑1497443

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