Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 364,170 2 bp→AA coding (2135‑2136/3075 nt) lacZ ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913364,1700GA92.9% 66.8 / ‑1.0 29A712V (GCC→GTT) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (10/16);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-02
*NC_000913364,1710GA93.1% 75.4 / ‑1.0 29A712V (GCC→GTT) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (11/16);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-01

ACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCGGCTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGAT  >  NC_000913/363945‑364410
                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                
aCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCg                                                                                                                                                                                                                          <  1:145702/251‑1 (MQ=255)
        aCTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACt                                                                                                                                                                                                                  <  1:49896/251‑1 (MQ=255)
           gATCGCGCAGAGGCGTCAGCAGTTGTGTTTTATCGCCAATACACATCTGTGAAAGAAAGCCGGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTa                                                                                                                                                                                                                 <  2:109311/249‑1 (MQ=255)
                 gcAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACt                                                                                                                                                                                                         <  1:42129/251‑1 (MQ=255)
                                            cGCCAATCCACATCCGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGATGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTTGCCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCtctc                                                                                                                                                                               <  1:262656/250‑1 (MQ=255)
                                                   ccACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCg                                                                                                                                                                       <  1:399174/251‑1 (MQ=255)
                                                           gggaaagaaaGCCTGCTTGGCGGTTAAATGGACAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTG‑GTGGTCAGATGCGGGATGGCATGGGACGCGGCGGGGAGCGTCTCACTGAGGTTTTCCGCCAGACTCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGGGGGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGtt                                                                                                                                                                 <  2:310309/247‑1 (MQ=255)
                                                                   aagCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTG‑GGGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCa                                                                                                                                                           <  1:298889/247‑1 (MQ=255)
                                                                          actggcgtttaactggcctacgactattaCCCAGCTCGATGCAAAAATCCATTTCTCTGCGCGG‑CAGGTGCGGGATGGCGGTGGACTCGGCGGGGAGCGTCAAACTGAGCTGTTCCGCCAGACGCCCCTGCTGCCAGGCTCTGCTGTGCCCAACTTCTGACCATGCGGTCGCTGTCGGTTGCTCTTCGAGTACTGTGCGCCCGAATTGCCCGTCGCTCTCCGGTTGCGGTAGTTCAGTCAGTTCAATCAAc                                                                                                                                                <  2:363775/227‑1 (MQ=255)
                                                                            gggCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCGGCTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTg                                                                                                                                              <  2:168357/250‑1 (MQ=255)
                                                                             ggcgttgaaaTTACCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTGCGCTG‑GTGGTCAGATGCGGGATGTCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTTCGGTAGTTCAGGCAGTGCAATCAACTGt                                                                                                                                             <  2:223561/244‑1 (MQ=255)
                                                                             ggCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGGTTTCCGCCAGACGCCCCTGCTGCCAGGCGCTGATGTGCCAAACTTCTGGCCCTGCGGGCGCGTTCGGTTGCACTACGCGTACTGTGAGCCCGAGTTGTCCGGCGCTCTCCGGCCGCGGTAGGTCAGGCCGGTCAATCAAATGt                                                                                                                                             >  1:44674/1‑251 (MQ=255)
                                                                                ggTTAAATTGCCAACGCTTATTTCCCAGCTCGATGCAAAAAGCCATTTCGCGG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACGCTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCGGCTTCTGACCATGCGGTCGCGTTCGGTTGTACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCTGTAGTTCAGGCAGTTCAATCAACTGTTTa                                                                                                                                          <  2:77816/251‑1 (MQ=255)
                                                                                                                         tCCATTTCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCTGCGGGGTGCGTCACCCTGAGGTTTTTCGCCAGACGCCACTGCTGGCAAGCGCTGATGTTCCCCAATTCTGGCCATGCGGTCGCGTTCGGTTGCCCTACGCGTACTGTGTGCCCGAGTTTCCCTGCGCTCTCCGGCTGCGGGAGTTTATTCAGTTCAATCATCCGTTTAACTTTTGGAGCGCCTTCCAGAGGGCATTTCACC‑CTTGcccg                                                                                                >  2:242022/1‑249 (MQ=255)
                                                                                                                           cATTTCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGATGAAGTGCCCAACTTCTGACCATGCGGGCGCGTTCGGTTGCACTACGCGGACTGTGAGACAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTTGGGAGCGACATCCAGGGG‑CCCTTCACCGATTGCCAGCg                                                                                              >  1:27007/1‑251 (MQ=255)
                                                                                                                             tttCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAAACGCCACTGCTGCCCGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGGTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCCGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGGGGCATCCGGAGG‑CACTTCCCCGCTTGCCAGccgc                                                                                            >  1:216395/1‑248 (MQ=255)
                                                                                                                                                     aTGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCCGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCCTcc                                                                      >  1:422267/1‑250 (MQ=255)
                                                                                                                                                               aCGCGGCGGGGAGCGTCACACTGAGGTTTTCCGACAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGc                                                           <  1:314805/251‑1 (MQ=255)
                                                                                                                                                                              tCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGCGCGACATCCAGCGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCCGCGCCACCATCCAGCGCAGGAGCTCGTTATCGCTATg                                             >  2:305720/1‑250 (MQ=255)
                                                                                                                                                                                    tGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGGTGTGCCCAACTTCTGCCCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTCCCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCTGCGGCTGACCATCCCGCGCCCCCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAAc                                      >  2:178466/1‑251 (MQ=255)
                                                                                                                                                                                     gAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAAc                                      <  1:439500/250‑1 (MQ=255)
                                                                                                                                                                                     gAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACa                                     <  1:311629/251‑1 (MQ=255)
                                                                                                                                                                                          tttCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTAt                                >  1:136050/1‑251 (MQ=255)
                                                                                                                                                                                          ttgccGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTa                                 <  1:205639/247‑1 (MQ=255)
                                                                                                                                                                                                   gACGCCACTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTc                       <  1:143542/251‑1 (MQ=255)
                                                                                                                                                                                                         aCTGCTGCCAGGCGCTGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCg                 <  1:180880/251‑1 (MQ=255)
                                                                                                                                                                                                            gctgcCAGGCGCTGATGAGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTCCGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATAGCTATGACGGAACCGGTATTCGCTGGTCACTTCgttg              >  2:21136/1‑248 (MQ=255)
                                                                                                                                                                                                               gcCAGGTGCAGATGAGCCCAACTTCTGACCATACGGTCCCGTCCGGTTGCACTACGCGTACTGTGAGCCAGAGATGCCCGGAGCTCTCCGCCAGCGGTAGTTCAGGCAGCTCCATCAACTGTTTCCCTTGGGGCGCGACATCCAGAGG‑CACTTAACCGCTGGCCAGCGGCTTACCATCCAGCGCCCCCACACAGGACAGGAACTAGTTAGCCCTATGACGGAACAGGttctagctggacccttgtctggtt           >  2:417689/1‑228 (MQ=255)
                                                                                                                                                                                                                        tGATGTGCCCAACTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGCACAGGTATTCGTTGGTCACTTCGATGGTTTGCCCGGAt  >  1:116694/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                
ACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTG‑GTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCACACTGAGGTTTTCCGCCAGACGCCACTGCTGCCAGGCGCTGATGTGCCCGGCTTCTGACCATGCGGTCGCGTTCGGTTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGG‑CACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGAT  >  NC_000913/363945‑364410

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