Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 366,519 G→A A331V (GCG→GTG)  lacI ← DNA‑binding transcriptional repressor LacI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,5190GA93.8% 37.4 / ‑7.2 16A331V (GCG→GTG) lacIDNA‑binding transcriptional repressor LacI
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (2/13);  total (3/13)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCA  >  NC_000913/366412‑366759
                                                                                                           |                                                                                                                                                                                                                                                    
caatcgcgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGGCCAGGGGGG‑GTTTTTTTTTTACCAGAGAGACGGGCAACCGATGATTGCCCCTC‑ACCCCCTCGCCCTGCGGGAGGTGTGGCAAGGCGG‑CCCCGCTGGGTTGGCCCCGCCAGG‑GAAAATCCAGTTTGGTTGTGGTTAACgggggg                                                                                                   >  2:325992/8‑251 (MQ=255)
    cgcgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑Tttttcttttc                                                                                                                                                                                                                                 <  1:200506/121‑1 (MQ=255)
             cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGc                                                                                      <  1:157271/251‑1 (MQ=255)
                     gCCCGCCTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGACCAGGGGGGGTTTTTTTTTTC‑CCAGTGAGACGGGCAAAAGCTGATTGCCCTTCCACC‑CCCGGCCCCTAGGGAGGTGCAGCAAG‑CGGTCCACGCTGG‑GTTGCCCCAGCAGGCGGAAATCCTGGTTTATGGTGGTTTACGGCggggtgttacataatgtgtattct                                                                              >  1:120365/1‑230 (MQ=255)
                      cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGTTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                                                                                                        <  1:167617/224‑1 (MQ=255)
                       ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAAGCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGt                                                                            <  1:357201/251‑1 (MQ=255)
                             tacaggcGGGAAACCTGTCGTGCCAGCTGCATTAATGAACCGGCCAACGCGGGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGCCGGGCAACAGCTGATTGCCCTTC‑TTCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCTAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTAtcgt                                                                       <  1:361901/244‑1 (MQ=255)
                              ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGGGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCGCCAGCAGGCGAAAATCCTGTTTTATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTAtcgtc                                                                      <  1:381993/250‑1 (MQ=255)
                                    gggAAACCGTTCGGGCCCGCTGCATTAATGAATCGGCCAACGCGCGGGGAGGGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATccc                                                               <  2:225312/251‑1 (MQ=255)
                                     ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGTCTTGAGATAACATGAGCTGTCTTCGGTATCGTCGTATCCCa                                                              <  1:52818/251‑1 (MQ=255)
                                                     aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCca                                                                                                                                         <  1:247804/160‑1 (MQ=255)
                                                                   aTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGcc                                <  1:161686/251‑1 (MQ=255)
                                                                             gcgcgGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGt                      <  1:103237/251‑1 (MQ=255)
                                                                                                 cGTATTGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCa  <  1:213100/251‑1 (MQ=255)
                                                                                                      tGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                                                                                                                 >  1:324221/1‑135 (MQ=255)
                                                                                                      tGGGCACCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                                                                                                                 <  2:324221/135‑1 (MQ=255)
                                                                                                           |                                                                                                                                                                                                                                                    
ATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCGCCAGGGTGG‑TTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTC‑ACCGCCTGGCCCTGAGAGAGTTGCAGCAAG‑CGGTCCACGCTGG‑TTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCA  >  NC_000913/366412‑366759

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑