Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,736,157 G→A *393* (TAG→TAA xylR → xylose divergent operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,736,1570GA98.0% 155.0 / ‑5.3 49*393* (TAG→TAAxylRxylose divergent operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/27);  minor base C (1/0);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCACCTTGTCGATCAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAGCGATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAA  >  NC_000913/3736022‑3736287
                                                                                                                                        |                                                                                                                                   
aaCCACCTTGTCGATCAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACga                                                                                                                                 <  2:1082960/139‑1 (MQ=255)
aaCCACCTTGTCGATCAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACga                                                                                                                                 >  1:344843/1‑139 (MQ=255)
              tCAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGt                                                                                                                    <  2:866632/139‑1 (MQ=255)
                     gaTATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTc                                                                                                             <  1:395952/139‑1 (MQ=255)
                        atCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          <  2:708385/139‑1 (MQ=255)
                        atCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          <  1:515246/139‑1 (MQ=255)
                          cGCAAATTTGCGGTTATCCGTCGCTGCAAGATTTCTACTCAGTTTTTAAAAAAGCATATGACACAACGCCAAAAGAGTATCGCGAAGTACATAGCCATGCCATGTTGTAACCAACGATGATAATTCTCTGCGTcataac                                                                                                        >  2:405862/1‑135 (MQ=255)
                                       ttATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           <  2:344843/139‑1 (MQ=255)
                                            cATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTgg                                                                                        <  2:371116/137‑1 (MQ=255)
                                            cATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTgg                                                                                        >  1:371116/1‑137 (MQ=255)
                                                    cAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATg                                                                              >  2:736784/1‑139 (MQ=255)
                                                            cTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATgcg                                                                            <  2:764916/133‑1 (MQ=255)
                                                            cTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATgcg                                                                            >  1:764916/1‑133 (MQ=255)
                                                              aCTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAAc                                                                    >  2:375753/1‑139 (MQ=255)
                                                                tctGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGc                                                                  <  2:329400/139‑1 (MQ=255)
                                                                 tgttttttAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          >  1:921758/4‑98 (MQ=255)
                                                                 tgttttttAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          <  2:921758/95‑1 (MQ=255)
                                                                 ctGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGcc                                                                 >  2:499069/1‑139 (MQ=255)
                                                                            aaaGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCt                                                      <  2:294145/139‑1 (MQ=255)
                                                                              aGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACgaatga                                                                                                                             <  1:759468/65‑5 (MQ=255)
                                                                              aGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACgaatga                                                                                                                             >  2:759468/1‑61 (MQ=255)
                                                                                          gacgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATgcg                                                                            >  1:989832/1‑103 (MQ=255)
                                                                                          gacgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATgcg                                                                            <  2:989832/103‑1 (MQ=255)
                                                                                            cgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAAtgatga                                                                                                                           >  1:191614/1‑54 (MQ=255)
                                                                                            cgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAAtgatga                                                                                                                           <  2:191614/54‑1 (MQ=255)
                                                                                                    gagTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTc                                                                                                                   >  1:1022826/1‑54 (MQ=255)
                                                                                                    gagTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTc                                                                                                                   <  2:1022826/54‑1 (MQ=255)
                                                                                                      gTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           >  2:846242/1‑76 (MQ=255)
                                                                                                      gTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           <  1:846242/76‑1 (MQ=255)
                                                                                                       tATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGAt                                                                               >  2:282524/1‑87 (MQ=255)
                                                                                                       tATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGAt                                                                               <  1:282524/87‑1 (MQ=255)
                                                                                                         tCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCg                                                                                                               >  2:520553/1‑53 (MQ=255)
                                                                                                         tCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCg                                                                                                               <  1:520553/53‑1 (MQ=255)
                                                                                                         tCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACcaac                         <  1:499069/139‑1 (MQ=255)
                                                                                                           gcgATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACcaacaa                       <  1:318276/139‑1 (MQ=255)
                                                                                                           gcgATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACcaacaa                       <  1:736784/139‑1 (MQ=255)
                                                                                                                       gCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATaa                                    >  1:104404/1‑114 (MQ=255)
                                                                                                                       gCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATaa                                    <  2:104404/114‑1 (MQ=255)
                                                                                                                        cGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATgcg                                                                            <  1:828293/73‑1 (MQ=255)
                                                                                                                        cGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATgcg                                                                            >  2:828293/1‑73 (MQ=255)
                                                                                                                        cGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCg          >  2:226124/1‑139 (MQ=255)
                                                                                                                         gAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAgg                                            >  2:373955/1‑104 (MQ=255)
                                                                                                                         gAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAgg                                            <  1:373955/104‑1 (MQ=255)
                                                                                                                               aTGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCa   >  2:1080092/1‑139 (MQ=255)
                                                                                                                                tgttgtAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGc                                                                  >  2:29269/1‑75 (MQ=38)
                                                                                                                                tgttgtAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGc                                                                  <  1:29269/75‑1 (MQ=38)
                                                                                                                                tgttgtAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGcc                                               <  2:683388/94‑1 (MQ=255)
                                                                                                                                tgttgtAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGcc                                               >  1:683388/1‑94 (MQ=255)
                                                                                                                                tgttgtAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCaa  <  1:375753/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                   
AACCACCTTGTCGATCAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAGCGATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAA  >  NC_000913/3736022‑3736287

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