Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,600,491 C→A A251S (GCA→TCA)  lsrR ← DNA‑binding transcriptional repressor LsrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,600,4910CA95.2% 60.0 / ‑3.1 21coding (751/954 nt)lsrRDNA‑binding transcriptional repressor LsrR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (10/10);  minor base . (0/1);  total (10/11)
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.58e-01

GCCGCTGCTGTGTCCTGATCGGTAACCAGTGCGTTGATATAACCGCCTTTCATTGCAGCGGCAATTGCTTCGGCTTTATTTTCTCCCCCTGCCACGCCAACCCGGACGGGTATGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGCATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCGACAATCGCCACATCCGCTGCTTGCGCGGCTAACAGAACATCTTTGACGCAATTTTCATTTTTTAGCGTACGGGCAATGTCAGCGGAGGAT  >  NC_000913/1600304‑1600708
                                                                                                                                                                                             |                                                                                                                                                                                                                          
gccgcTGCGTGGTCCTGATCGGTACCCATTGTGTTGATATAACCGCCTTTCATGGCACCGGCAATTGCTTCGCTGTTATTTTCTCCCCCTGCCACGCCTCCCCGGACGGGTATGGTCTTCAGCGCGCTTAAATGTAGGCCAATCAGTTCGGTAGGTATTTTGATTTTCGTGACAGCGTCACTTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAAGCATTAATTGTTCGcc                                                                                                                                                            <  1:164481/251‑1 (MQ=255)
 ccgcTGCTGTGTCCTGATCGGTAACCAGTGCGTTGATATAACCGCCTTTCATTGCAGCGGCAATTGCTTCGGCTTTATTTTCTCCCCCTGCCACGCCAACCCGGACGGGTATGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGccc                                                                                                                                                           >  2:50918/1‑251 (MQ=255)
          gtgcccTGATCGGTACCCAGTGTGTTGATATAACCGCCTTTTATTGCAGCGGCAATTGCTTCGGCTTTATTTTCTCCCCCTGCCACGCCAACCGGGACGGTTATGGTCTTCAGCGCGCTTAAAGTTAAGCCAATCATTTCGTGATGTATTTTGATATTCGTGACAACGCCCCCTTTTGATTAAATAAAGTATCCTAAAATGCCGCCAACCGCCCCTTTTCGGCCAATCATTTACTGTTCGCCCTGGCTGat                                                                                                                                                   <  1:32789/247‑1 (MQ=255)
           gtCCTGTTCGGTAACCAGTGCGTTGATATACCCGCCTTTCATTGCAGCGGCAATTGCTTCGGCTTTATTTTCTCCCCCTGCCACGCCAACCCGGACGGGTATGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGata                                                                                                                                                  <  1:50918/250‑1 (MQ=255)
              ccTGATCGGTAACCAGTGCGTTGATATAACCGCCTTTCATTGCAGCGGCAATTGCTTCGGCTTTATTTTCTCCCCCTGCCACGCCAACCCGGACGGGTATGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGAtata                                                                                                                                                >  1:88892/1‑250 (MQ=255)
                        aaCCGGTGCGTTGATATACCCGCCTTTCATTGCAGCGGCATTGGTTTCGGCTTTATTTTCTCCCCCTGCCACGCCAACCCGGACGGGTATGGTCTTCAGCGCGCTTAAAGGTAAGCAAATCAGTTCGTTAGGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTGCGGCCAATCATTAACTGTTCGCCCTGGCTTATATAACCGTAGCGaa                                                                                                                                     <  2:209358/251‑1 (MQ=255)
                                       ataACCGCCTTTCATTGCAGCGGCAATTGCTTGGGTTTTATTTTCTCCCCCGGCCACGCCAACCCGGACGGGTATGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTGTGATATGCGGGACAACGTCACCATTTGAATCAAAAAAGTAGCCAAAAATTTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATAGAACCGGAGCGAATGATTGTCTCATCGt                                                                                                                      <  2:108977/251‑1 (MQ=255)
                                       ataACCGCCTGTCATTGCAGCGGCAATTGCTTCGGCTTTATTTTCTCCCCCTGCCACGCCAACCCGGACGGGTATGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTGCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCGTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGTCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCg                                                                                                                       <  1:123789/250‑1 (MQ=255)
                                                          gCGGCAATTGCTTCGGCTTTATTTTCTCCCCCTGCCACGCCAACCCGGACGGGTATGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCAcc                                                                                                   >  2:30610/1‑251 (MQ=255)
                                                           gggacatggcgtCGGCTTTATTGGCTCCCCCTGCCCCGCCAACCCGGTCGGGTTTGGTCTTTCGCGCCCTGAAAGGGAAGCCAATCAGTTCTTTAATTATTTTGATTATCGGGACAACGTCACCTTTTTAATCAAAAAAGTAGCCTAAAAGGTCGCCATCCGCCCCTTTCCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATTATTGTCGCATAGTCCTGTTTACTCACATCACCa                                                                                                  <  2:23294/240‑1 (MQ=255)
                                                                  tGCTTCGGCTTTATTTTCTCCCCCTGCCACGCCAACCCGGACGGGTATGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCg                                                                                            >  1:89354/1‑250 (MQ=255)
                                                                                    ctcCCCCTGCCACGCCAACCCGGACGGGTATGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCGACAATCGCCACATCCgct                                                                          >  1:134260/1‑251 (MQ=255)
                                                                                                  ccAACCCGGACGGGTAGGGGCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCGACAATCGCCACATCCGCTGCTTGCGCGGCTaa                                                            <  1:196185/251‑1 (MQ=255)
                                                                                                             gggTATGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCGACAATCGCCACATCCGCTGCTTGCGCGGCTAACAGAACATCtt                                                 >  1:46582/1‑251 (MQ=255)
                                                                                                              ggTATGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCGACAATCGCCACATCCGCAGCTTGCGCGGCTAACAGAACATCttt                                                >  2:146903/1‑251 (MQ=255)
                                                                                                                 aTGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCGACAATCGCCACATCCGCTGCTTGCGCGGCTAACAGAACATCTTTGa                                              <  1:174543/250‑1 (MQ=255)
                                                                                                                               gcTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCGACAATCGCCACATCCGCTGCTTGCGCGGCTAACAGAACATCTTTGACGCAATTTTCAtttt                               >  1:112297/1‑251 (MQ=255)
                                                                                                                                      ggTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCGACAATCGCCACATCCGCTGCTTGCGCGGCTAACAGAACATCTTTGACGCAATTTTCATTTTTTAGCGt                        >  1:33662/1‑251 (MQ=255)
                                                                                                                                         aaGCCAATCAGTTCGTTATGTATTTTGATATTCTTGCAAACGTCCCCTTTTGAATCAAAAAAGTAGCCTAAAAGGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCGACAATCGCCACATCCGCTGCTTGCGCGGCTAACAGAACATCTTTGACGCAATTTTCATTTTTTAGCGTAc                      <  2:81116/250‑1 (MQ=255)
                                                                                                                                           cccAACCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGAAACAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCGACAATCGCCACATCCGCTGCTTGCGCGGCTAACAGAACATCTTTGACGCAATTTTCATTTTTTAGCGTACggg                   <  1:230822/250‑1 (MQ=255)
                                                                                                                                                            gTATTTTGATATTCGTGACAACGTCACCTTTTGAATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCGACAATCGCCACATCCGCTGCTTGCGCGGCTAACAGAACATCTTTGACGCCATTTTCATTTTTTAGCGTACGGGCAATGTCAGCGGAGGAt  >  1:119228/1‑251 (MQ=255)
                                                                                                                                                                                             |                                                                                                                                                                                                                          
GCCGCTGCTGTGTCCTGATCGGTAACCAGTGCGTTGATATAACCGCCTTTCATTGCAGCGGCAATTGCTTCGGCTTTATTTTCTCCCCCTGCCACGCCAACCCGGACGGGTATGGTCTTCAGCGCGCTTAAAGGTAAGCCAATCAGTTCGTTATGTATTTTGATATTCGTGACAACGTCACCTTTTGCATCAAAAAAGTAGCCTAAAATGTCGCCAACCGCCCCTTTTCGGCCAATCATTAACTGTTCGCCCTGGCTGATATAACCGGAGCGAATGATTGTCGCATCGTCCTGTTGACTCACAGCACCAATGCCGACAATCGCCACATCCGCTGCTTGCGCGGCTAACAGAACATCTTTGACGCAATTTTCATTTTTTAGCGTACGGGCAATGTCAGCGGAGGAT  >  NC_000913/1600304‑1600708

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