Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 713,690 G→A 10.0% L94L (CTG→CTA seqA → regulatory protein for replication initiation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb713,6900GA10.0% 183.4 / 9.3 60L94L (CTG→CTAseqAregulatory protein for replication initiation
Reads supporting (aligned to +/- strand):  ref base G (22/32);  new base A (3/3);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACAC  >  W3110S.gb/713624‑713760
                                                                  |                                                                      
gCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                    >  1:1691132/1‑71 (MQ=255)
gCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCtgttgt                                                                    >  1:482381/1‑71 (MQ=255)
 cGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                   >  1:1693871/1‑71 (MQ=255)
 cGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                   >  1:1627325/1‑71 (MQ=255)
 cGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                   >  1:1620705/1‑71 (MQ=255)
 cGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                   >  1:809060/1‑71 (MQ=255)
 cGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                   >  1:1717100/1‑71 (MQ=255)
 cGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                   >  1:900135/1‑71 (MQ=255)
 cGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                   >  1:195166/1‑71 (MQ=255)
 cGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTc                                                                   >  1:1381066/1‑71 (MQ=255)
   tGAACTTCTGCTTTCGGATGAATTCGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTa                                                                 <  1:1064849/71‑1 (MQ=255)
   tGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTa                                                                 <  1:1233456/71‑1 (MQ=255)
   tGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTa                                                                 <  1:736023/71‑1 (MQ=255)
   tGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTa                                                                 <  1:166744/71‑1 (MQ=255)
   tGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTa                                                                 <  1:779114/71‑1 (MQ=255)
     aaCTTCTGCTTTCGGATGAATACGCAGAGCAAAGGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTaca                                                               >  1:1812856/1‑71 (MQ=255)
     aaCTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTaca                                                               >  1:1532612/1‑71 (MQ=255)
     aaCTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTaca                                                               >  1:1539935/1‑71 (MQ=255)
     aaCTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTaca                                                               >  1:619017/1‑71 (MQ=255)
     aaCTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTaca                                                               >  1:1157192/1‑71 (MQ=255)
     aaCTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTaca                                                               >  1:682023/1‑71 (MQ=255)
     aaCTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTaca                                                               >  1:1397169/1‑71 (MQ=255)
          cTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACtat                                                           <  1:1616534/70‑1 (MQ=255)
             cTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTATTGTCTACACtata                                                          <  1:1571492/68‑1 (MQ=255)
                  ggATGTATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTATTGTCTACACTATATTCTCtt                                                   <  1:1751642/70‑1 (MQ=255)
                  ggATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTATTGTCTACACTATATTCTCtt                                                   <  1:1689173/70‑1 (MQ=255)
                      gAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGc                                              >  1:192164/1‑71 (MQ=255)
                      gAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGc                                              >  1:1708766/1‑71 (MQ=255)
                         tACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                   <  1:610356/63‑1 (MQ=255)
                         tACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                   <  1:1198491/63‑1 (MQ=255)
                         tACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                   <  1:1803716/63‑1 (MQ=255)
                         tACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGa                                                 >  1:36098/1‑65 (MQ=255)
                         tACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGa                                                 >  1:615601/1‑65 (MQ=255)
                           cGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAg                                          >  1:897176/1‑70 (MQ=255)
                            gCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                   <  1:1616871/60‑1 (MQ=255)
                             cAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCg                                       <  1:1537950/71‑1 (MQ=255)
                             cAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCg                                       <  1:52812/71‑1 (MQ=255)
                             cAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCg                                       <  1:808014/71‑1 (MQ=255)
                             cAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCAACACTATATTCTCTTGACGCCCAGGCg                                       <  1:1417733/71‑1 (MQ=255)
                              agagCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCg                                       <  1:1127087/70‑1 (MQ=255)
                                    aaaGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACta                                                            <  1:816325/43‑1 (MQ=255)
                                        cgagcgGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGc                            <  1:1531972/71‑1 (MQ=255)
                                        cgagcgGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGc                            <  1:1783647/71‑1 (MQ=255)
                                        cgagcgGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGc                            <  1:1099146/71‑1 (MQ=255)
                                         gagcgGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGc                            <  1:319443/70‑1 (MQ=255)
                                             ggTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACgg                       <  1:986574/71‑1 (MQ=255)
                                              gTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                   <  1:1594077/42‑1 (MQ=255)
                                              gTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGa                      <  1:43605/71‑1 (MQ=255)
                                              gTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGa                      <  1:490819/71‑1 (MQ=255)
                                              gTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGa                      <  1:578736/71‑1 (MQ=255)
                                              gTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGa                      <  1:218128/71‑1 (MQ=255)
                                               tCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                   <  1:402453/41‑1 (MQ=255)
                                               tCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGa                      <  1:341115/70‑1 (MQ=255)
                                                cAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAAt                    >  1:1835944/1‑71 (MQ=255)
                                                 aaTCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                   <  1:133914/71‑1 (MQ=255)
                                                 aaTCGCTTTATGCTGCTATTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCTAAGCAACGGAATc                   >  1:1799453/1‑71 (MQ=255)
                                                 aaTCGCTTTATGCTGCTATTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAAt                    >  1:1797354/1‑70 (MQ=255)
                                                 aaTCGCTTTATGCTGCTATTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATc                   >  1:1342141/1‑71 (MQ=255)
                                                  aTCGCTTTATGCTGCTGTTGTCTACACTATATTCTCtt                                                   <  1:640390/38‑1 (MQ=255)
                                                     gCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCg                  >  1:1329760/1‑68 (MQ=255)
                                                          aTGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACg           <  1:694807/70‑1 (MQ=255)
                                                          aTGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACg           <  1:1717045/70‑1 (MQ=255)
                                                             ctgctgTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCg       <  1:710834/71‑1 (MQ=255)
                                                             ctgctgTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCg       <  1:1312875/71‑1 (MQ=255)
                                                             ctgctgTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCg       <  1:1082623/71‑1 (MQ=255)
                                                                  gttgtCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTacac  >  1:510312/1‑71 (MQ=255)
                                                                  |                                                                      
GCGTGAACTTCTGCTTTCGGATGAATACGCAGAGCAAAAGCGAGCGGTCAATCGCTTTATGCTGCTGTTGTCTACACTATATTCTCTTGACGCCCAGGCGTTTGCCGAAGCAACGGAATCGTTGCACGGTCGTACAC  >  W3110S.gb/713624‑713760

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