Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,632,877 A→G 37.5% Y294C (TAC→TGC)  nadR → bifunctional DNA‑binding transcriptional repressor and NMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,632,8770AG37.5% 15.9 / 22.5 24Y294C (TAC→TGC) nadRbifunctional DNA‑binding transcriptional repressor and NMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (8/7);  new base G (4/5);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACT  >  W3110S.gb/4632809‑4632941
                                                                    |                                                                
ccTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATAc                                                                 <  1:1104823/70‑1 (MQ=255)
           atgaGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCg                                                     >  1:361957/1‑71 (MQ=255)
           atgaGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACAt                                                               >  1:775913/1‑61 (MQ=255)
           atgaGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGc                                                      <  1:1369230/70‑1 (MQ=255)
                    cATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAatat                                            <  1:1213922/71‑1 (MQ=255)
                        gCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGcc                                         >  1:625207/1‑70 (MQ=255)
                         cAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCa                                        <  1:667072/70‑1 (MQ=255)
                         cAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCa                                        <  1:1667893/70‑1 (MQ=255)
                            tATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATaa                                    <  1:558867/71‑1 (MQ=255)
                              ttCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGaaa                                               >  1:557982/1‑58 (MQ=255)
                                cTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGAtt                                                          >  1:1434476/1‑45 (MQ=255)
                                 tGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGa                                                 <  1:89570/53‑1 (MQ=255)
                                 tGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGaaa                                               >  1:1631581/1‑55 (MQ=255)
                                  gACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTgg                               >  1:579079/1‑70 (MQ=255)
                                   aCTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGc                              <  1:382298/70‑1 (MQ=255)
                                        gATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATa                                     <  1:966344/58‑1 (MQ=255)
                                          tAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATaaa                                   >  1:1633356/1‑58 (MQ=255)
                                          tAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATaaa                                   >  1:956/1‑58 (MQ=255)
                                           aaaaTCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTgg                               >  1:496074/1‑61 (MQ=255)
                                                 gcgcTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCACTTTATTGGCAtat                          <  1:1083208/60‑3 (MQ=25)
                                                 gcgcTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGc                              >  1:1458178/1‑56 (MQ=255)
                                                        gCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGt         >  1:143175/1‑70 (MQ=255)
                                                        gCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTc        <  1:716406/71‑1 (MQ=255)
                                                              cTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACt  <  1:895410/71‑1 (MQ=255)
                                                                    |                                                                
CCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACT  >  W3110S.gb/4632809‑4632941

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