Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 468,510 A→G E43G (GAG→GGG)  ybaO → putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913468,5100AG97.8% 153.4 / ‑4.6 45E43G (GAG→GGG) ybaOputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (25/19);  total (26/19)
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.21e-01

GGCTATGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTGCCGCTTTGTCACG  >  NC_000913/468379‑468637
                                                                                                                                   |                                                                                                                                
ggCTATGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                               <  2:475467/135‑1 (MQ=255)
ggCTATGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                               >  1:475467/1‑135 (MQ=255)
  cTATGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACgg                                                                                                                         >  1:116479/1‑139 (MQ=255)
    aTGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTAt                                                                                                                      >  2:556887/1‑139 (MQ=255)
                    aCCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                               >  2:519498/1‑115 (MQ=255)
                    aCCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                               <  1:519498/115‑1 (MQ=255)
                    aCCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGt                                                                                                      >  2:1473677/1‑139 (MQ=255)
                                  gCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATc                                                                                        >  2:271756/1‑139 (MQ=255)
                                      tACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGa                                                                                    <  1:1473677/139‑1 (MQ=255)
                                            agcagGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATc                                                                                        >  1:1474793/1‑129 (MQ=255)
                                            agcagGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATc                                                                                        <  2:1474793/129‑1 (MQ=255)
                                                 gATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCg                                                                                                    <  2:682189/112‑1 (MQ=255)
                                                 gATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCg                                                                                                    >  1:682189/1‑112 (MQ=255)
                                                   ttGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTc                                                                       >  2:378951/1‑139 (MQ=255)
                                                      cACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTa                                                                                                                       >  2:803942/1‑89 (MQ=255)
                                                      cACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTa                                                                                                                       <  1:803942/89‑1 (MQ=255)
                                                        cccTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCt                                                                  <  2:116479/139‑1 (MQ=255)
                                                                     ggCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCg                                                                                                    <  1:221724/92‑1 (MQ=255)
                                                                     ggCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCg                                                                                                    >  2:221724/1‑92 (MQ=255)
                                                                     ggCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGcc                                                                                                  <  1:937465/94‑1 (MQ=255)
                                                                     ggCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGcc                                                                                                  >  2:937465/1‑94 (MQ=255)
                                                                            gCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATaa                                              >  2:162086/1‑139 (MQ=255)
                                                                                 aGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACg                                         <  1:162086/139‑1 (MQ=255)
                                                                                      ttAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACgg                                                                                                                         <  2:1084421/55‑1 (MQ=255)
                                                                                      ttAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACgg                                                                                                                         >  1:1084421/1‑55 (MQ=255)
                                                                                      ttAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCgg                                                                                                            >  2:1247499/1‑67 (MQ=255)
                                                                                      ttAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCgg                                                                                                            <  1:1247499/67‑1 (MQ=255)
                                                                                      ttAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCg                                                                                                    <  2:1042619/75‑1 (MQ=255)
                                                                                      ttAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCg                                                                                                    >  1:1042619/1‑75 (MQ=255)
                                                                                                     ccccTTGCTGGAAGCGCCTGAAACGGCTGGGGGacgacg                                                                                                                          <  1:616729/39‑1 (MQ=255)
                                                                                                     ccccTTGCTGGAAGCGCCTGAAACGGCTGGGGGacgacg                                                                                                                          >  2:616729/1‑39 (MQ=255)
                                                                                                     ccccTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATg                     >  1:726138/1‑139 (MQ=255)
                                                                                                     ccccTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATg                     >  1:384597/1‑139 (MQ=255)
                                                                                                     ccccTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATg                     >  1:451745/1‑139 (MQ=255)
                                                                                                      cccTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATgg                    >  2:1290003/1‑139 (MQ=255)
                                                                                                        cTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTa                  >  1:602967/1‑139 (MQ=255)
                                                                                                          tGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTAtt                >  1:590395/1‑139 (MQ=255)
                                                                                                           gCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTg               >  2:776622/1‑139 (MQ=255)
                                                                                                               gAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCtt                                                                                                                 >  1:590366/1‑37 (MQ=255)
                                                                                                               gAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCtt                                                                                                                 <  2:590366/37‑1 (MQ=255)
                                                                                                                 aGCGCCTGAAACGGCTGGGGGACGACGGTATCCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTGCCGCt          <  1:556887/139‑1 (MQ=255)
                                                                                                                     ccTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGAc                                                               <  1:1467739/81‑1 (MQ=255)
                                                                                                                     ccTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGAc                                                               >  2:1467739/1‑81 (MQ=255)
                                                                                                                     ccTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCaa                                      <  2:97317/106‑1 (MQ=255)
                                                                                                                     ccTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCaa                                      >  1:97317/1‑106 (MQ=255)
                                                                                                                      cTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCt                                                   >  2:18747/1‑92 (MQ=255)
                                                                                                                      cTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCt                                                   <  1:18747/92‑1 (MQ=255)
                                                                                                                        gAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTGCCGCTTTGTCACg  <  2:451745/139‑1 (MQ=255)
                                                                                                                        gAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTGCCGCTTTGTCACg  <  1:378951/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                
GGCTATGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTGCCGCTTTGTCACG  >  NC_000913/468379‑468637

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