Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,667,460 G→A Q369* (CAG→TAG)  mlc ← DNA‑binding transcriptional repressor Mlc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,667,4600GA95.0% 50.6 / NA 20Q369* (CAG→TAG) mlcDNA‑binding transcriptional repressor Mlc
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/9);  minor base C (1/0);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAA  >  NC_000913/1667223‑1667682
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                              
tAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTAATCTCAACGCTGATGTGCTAACTATACGCAgg                                                                                                                                                                                                                    >  1:61980/1‑250 (MQ=255)
                 ttCACAAATATAAACCAGGAAAATAATTACCCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGAt                                                                                                                                                                                                  >  2:134798/1‑251 (MQ=255)
                      aaaTATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGt                                                                                                                                                                                             <  1:108681/251‑1 (MQ=255)
                               ccaagAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAGCCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGCCGGATGCTGTCGGAGATGa                                                                                                                                                                                    <  2:260275/247‑1 (MQ=255)
                                     aaaTAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGa                                                                                                                                                                              <  2:188352/251‑1 (MQ=255)
                                                                              cAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGTTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATc                                                                                                                                        >  2:91912/1‑248 (MQ=255)
                                                                              cAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAg                                                                                                                                      <  1:91912/250‑1 (MQ=255)
                                                                                   tGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACACAGAACCGTTATACATCGCGTCTTTTAACAGTGAAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTATCAACGCTGATGTGCTCACTATACGCAGGAAGGGCCTGCTGACGGATGCTGACTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACTGTGAGCCACTCAGTAtttt                                                                                                                                >  2:119616/1‑251 (MQ=255)
                                                                                                       tAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCCGCCCTCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGCCCTGCCGACGGATGCTGTCTGAGATGACCGGGCAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATACATTc                                                                                                             >  1:138199/1‑250 (MQ=255)
                                                                                                                           aaCCCTGCAACATACGAAGCAACAAAGAACCTTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCGCATCGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAATAATgcg                                                                                         <  2:115718/250‑1 (MQ=255)
                                                                                                                                            aTCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGACTGCGCCCGACATGCGCGCCCAc                                                                       >  2:223838/1‑251 (MQ=255)
                                                                                                                                                     gAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGCCGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAat                                                              >  1:241185/1‑251 (MQ=255)
                                                                                                                                                                      gTCTTTTACCAGTGCAGCGCCTGCCATCGAGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGACGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGCGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCAGTTAATGCTGTCTTTTGCCAGt                                              >  1:132676/1‑250 (MQ=255)
                                                                                                                                                                                            gCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCt                        <  1:221254/250‑1 (MQ=255)
                                                                                                                                                                                            gCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTg                       <  1:168854/251‑1 (MQ=255)
                                                                                                                                                                                                    gCCCTGGTTAGAAAACTGAGTACTCTCCACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATgccg                           >  2:183333/1‑239 (MQ=255)
                                                                                                                                                                                                    gCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATgccg                           <  1:183333/239‑1 (MQ=255)
                                                                                                                                                                                                     cccTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATg               >  1:106213/1‑250 (MQ=255)
                                                                                                                                                                                                               aaaaCTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACCCAATGAGTCCACGGtt    >  1:228777/1‑251 (MQ=255)
                                                                                                                                                                                                                 aaCTGAGTACTCTCAACGCTTATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGATGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCTCAATGCCGCCTGACACAATGAGTCCACGGTTaa  <  2:258538/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                              
TAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAA  >  NC_000913/1667223‑1667682

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