Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,277,850 A→G 11.4% intergenic (+448/‑498) pyrH → / → uppS UMP kinase/di‑trans, poly‑cis‑decaprenylcistransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,277,8500AG11.4% 79.2 / 4.2 36intergenic (+448/‑498)pyrH/uppSUMP kinase/di‑trans, poly‑cis‑decaprenylcistransferase
Reads supporting (aligned to +/- strand):  ref base A (15/16);  new base G (2/2);  total (17/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 = 6.79e-01

CCTGCCTTAACAGAAGAACGTAGAAAAGAGCGCGTTAAAGATGTTAAGAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATAAAGAAAAAGATATTATGTCAGTATAAAACTAATATACAATG  >  CP000730/1277724‑1277990
                                                                                                                              |                                                                                                                                            
ccTGCCTTAACAGAAGAACGTAGAAAAGTGCGCGTTAAAGATGTTAAGAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGTGATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACt                                                                                                                                >  1:98304/1‑141 (MQ=255)
   gCCTTAACAGAAGAACGTAGAAAAGAGCGCGTTAAAGATGTTAAGAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACGGaa                                                                                                                             >  2:244600/1‑141 (MQ=255)
     cTTAACAGAAGAACGAAGAAAAGAGGGCGTTAAAGATGTTAAGAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGa                                                                                                                           <  2:103055/141‑1 (MQ=255)
        aaCAGCAGAACGTAGAAAAGAGCGCGTTAAAGATGTTAACAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGAACAGTTGAAAACAGATGAAAAAAATGGCGACATTACTGAAGATGa                                                                                                                        >  2:81603/1‑141 (MQ=255)
        aaCAGAAGACCGTAGCCACGAGCGCGTCACAGCTGTTAAGAAAATCGGTGAAGCAGCTAAAGTATCTGTTCGAAATCTTCGTCGTGATCTGACTGATCAGTTGAAACAAGCTGCAACCAATGGCCACATCACTGACGCTGa                                                                                                                        >  2:349622/1‑141 (MQ=255)
        aaCAGAAGAACGTAGAAAAGAGCGCGTTAAAGATGTTAAGAAAATTGGTGACGAAGCTAAAGTATCTGTTAGAAAAATTCGTCGTGATATGAATGCTAAGTTGAACAAAGATGAAAAAAATGGCGACATTACTGAAGATGa                                                                                                                        >  2:168444/1‑141 (MQ=255)
          cAGAAGAACGTAGAAAAGAGCGCGTTAAAGATGTTAAGAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGaggcgg                                                                                                                      >  1:263153/1‑136 (MQ=255)
              agaaCGTAGAAAAGAGCGCGTTAAAGATGTTAAGAAAATGGGTGAAGAAGCTAAAGTATCGGTTCGAAATATTCATCCTGATATGAATGATCAGTTGAAAAAAGATGAAAAAACTGGACACATTACTGAAGAGGAGATgac                                                                                                                  >  1:135857/1‑140 (MQ=255)
              agaaCGTAGAAAAGAGCGCGATAAAGATGTTAAGAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGAGCAGGTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTgag                                                                                                                  <  1:32804/141‑1 (MQ=255)
                  cGTAGAAAAGGGCGCTGTAAAGATGTTAAGAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGt                                                                                                              <  1:313664/141‑1 (MQ=255)
                  cGTAGAAAAGGGCGCGTTTCTAGTGGTAATGAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGGAGATGAGTTGAGAAGt                                                                                                              <  2:172130/141‑1 (MQ=255)
                  cGTAGAAAAGAGCGCGTTAAAGATGTTAAGAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGt                                                                                                              <  1:34594/141‑1 (MQ=255)
                     aGAAAAGAGCGCGTTAAAGATGTTAAGAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAGATGGCGACATTACTGAAGATGAGTTGAGAAGTGGt                                                                                                           >  1:305050/1‑141 (MQ=255)
                     aGAAAAGAGCGCGTTAAAGATGTTAAGAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAGATGGCGACATTACTGAAGATGAGTTGAGAAGTGGt                                                                                                           >  1:208520/1‑141 (MQ=255)
                           gagCGCGTTAAAGATGTTAAGAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGaa                                                                                                     <  2:137376/141‑1 (MQ=255)
                                            tAAGAAAGTTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAAc                                                                                    <  2:90280/141‑1 (MQ=255)
                                                            gaagCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATaaaa                                                                    >  1:149039/1‑141 (MQ=255)
                                                             aagCGAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGAGCAGTTGAAAAAAGATGAACAATCTGGCGACACTACCGAAGATTAGTTGAGATGTGGTACTGAAGATGTTCACAAAGCAACAGACAATTCAATAAAAg                                                                   <  2:217/141‑1 (MQ=255)
                                                                cTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGaaa                                                                >  1:77978/1‑141 (MQ=255)
                                                                cTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTGGCGAGGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGaaa                                                                >  1:106115/1‑141 (MQ=255)
                                                                    aGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGa                                                            <  1:122772/141‑1 (MQ=255)
                                                                          ggTTCGAAATATTCGTCGTGGTATGAATGATAGGTGGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAAt                                                      <  2:5538/140‑1 (MQ=255)
                                                                                  atatACGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTg                                              >  2:264965/1‑141 (MQ=255)
                                                                                         tcgtGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATAAAGa                                       <  1:346896/141‑1 (MQ=255)
                                                                                         cccggATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGACTGCTGATAAAGa                                       <  1:196253/137‑1 (MQ=255)
                                                                                           gtGATATGAATGATCAGTTGAAAAAAGATGAAAAAGATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATAAAGaaa                                     <  2:208520/141‑1 (MQ=255)
                                                                                           gtGATATGAATGATCAGTTGAAAAAAGATGAAAAAGATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATAAAGaaa                                     <  2:305050/141‑1 (MQ=255)
                                                                                                atgaatgaTCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGACTGCTGATAAAGAAAAAGat                                <  1:210092/141‑1 (MQ=255)
                                                                                                 ttaatgaTCAGTTGAAAAAAAATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATAAAGAAAAAGata                               <  2:51059/139‑1 (MQ=255)
                                                                                                 tgaatgcTCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCGATAAAAGCAATTGATCAACTGCTTGCTCATAAAGAAAAAGata                               >  2:323023/1‑141 (MQ=255)
                                                                                                  gaatgaTCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATAAAGAAAAAGatat                              <  1:293859/141‑1 (MQ=255)
                                                                                                        tCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATAAAGAAAAAGATATTATGTc                        <  1:18687/141‑1 (MQ=255)
                                                                                                          aGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATAAAGAAAAAGATATTATGTCAg                      <  1:109096/141‑1 (MQ=255)
                                                                                                                   aaGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATACAGAAAAAGATATTATGTCAGTATAAAACt             >  2:78301/1‑141 (MQ=255)
                                                                                                                    aGATTAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATAAAGAAAAAGATATTATTTCAGCATAAAACTa            >  2:65162/1‑141 (MQ=255)
                                                                                                                        gAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATAAAGAAAAAGATATTATGTCAGTATAAAACTAatat        >  1:237112/1‑141 (MQ=255)
                                                                                                                              aaTGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATAAAGAAAAAGATATTATGTCAGTATAAAACTAATATACAAta  >  2:377641/1‑140 (MQ=255)
                                                                                                                              |                                                                                                                                            
CCTGCCTTAACAGAAGAACGTAGAAAAGAGCGCGTTAAAGATGTTAAGAAAATTGGTGAAGAAGCTAAAGTATCTGTTCGAAATATTCGTCGTGATATGAATGATCAGTTGAAAAAAGATGAAAAAAATGGCGACATTACTGAAGATGAGTTGAGAAGTGGTACTGAAGATGTTCAGAAAGCAACAGACAATTCAATAAAAGAAATTGATCAAATGATTGCTGATAAAGAAAAAGATATTATGTCAGTATAAAACTAATATACAATG  >  CP000730/1277724‑1277990

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