Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,601,098 C→A S48S (TCG→TCT)  lsrR ← DNA‑binding transcriptional repressor LsrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,601,0980CA96.2% 81.4 / ‑6.5 26S48S (TCG→TCT) lsrRDNA‑binding transcriptional repressor LsrR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (10/15);  minor base G (0/1);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATGGTTGCCTCGCCAAAACCAATCGCCAGCATCTGTTGTGGTTGAAGTAAACTCATCAACATATGCGCCGCGCCTATCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAATTCATTCTTCACTTTGAACATATTTAAATCTTTAATGCAATTGTTCAGTTCTTGCTCATTTATATCTGTGATGGCAACCACAGTTTGACTCT  >  NC_000913/1600854‑1601335
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
tCATGGTTGCCTCGCCAAAACCAATCGCCAGCATCTGTTGTGGTTGAAGTACACTCATCAACATATGCGCCGCGCCTATCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAACCTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGacatt                                                                                                                                                                                                                                         >  1:437340/1‑249 (MQ=255)
 cATGGTTGCCTCGCCAAAACCAATCGCCAGCATCTGTTGTGGTTGAAGTAAACTCATCAACATATGCGCCGCGCCTATCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGacattt                                                                                                                                                                                                                                        >  1:379660/1‑248 (MQ=255)
              ccAAAACCAATCGCCAGCATCTGTTGTGGTTGAAGTAAACTCATCAACATATGCGCCGCGCCTATCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAg                                                                                                                                                                                                                           <  1:366436/251‑1 (MQ=255)
                                                    aCTCATCAACATATGCGCCGCGCCTATCCCCAGTCGCCCACCGACCTCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCcgt                                                                                                                                                                                     >  2:84228/1‑251 (MQ=255)
                                                                          ccTATCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGAt                                                                                                                                                               <  1:72579/251‑1 (MQ=255)
                                                                             aTCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATc                                                                                                                                                              <  2:379660/249‑1 (MQ=255)
                                                                                cccAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCcgcgc                                                                                                                                                         >  1:188130/1‑251 (MQ=255)
                                                                                         ccACCGACATCAGCCTCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGttc                                                                                                                                                >  2:181420/1‑251 (MQ=255)
                                                                                         ccACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGttc                                                                                                                                                <  2:437340/251‑1 (MQ=255)
                                                                                                       aTCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTc                                                                                                                                  >  2:151460/1‑251 (MQ=255)
                                                                                                                cccAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAACCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGGTGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGCCGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTg                                                                                                                         >  1:158214/1‑251 (MQ=255)
                                                                                                                   aGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGaaa                                                                                                                      <  1:151460/251‑1 (MQ=255)
                                                                                                                                      tGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTccc                                                                                                                                >  1:299669/1‑222 (MQ=255)
                                                                                                                                      tGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTccc                                                                                                                                <  2:299669/222‑1 (MQ=255)
                                                                                                                                        cAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCa                                                                                                 <  2:188130/251‑1 (MQ=255)
                                                                                                                                                   tacgacgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCTACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAattcattc                                                                                       <  1:389763/249‑1 (MQ=255)
                                                                                                                                                   gacgacGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTGCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAattcatt                                                                                        <  2:65916/249‑1 (MQ=255)
                                                                                                                                                       acgTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATACTGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAATTCATTCTTCAc                                                                                  >  2:229398/1‑251 (MQ=255)
                                                                                                                                                                  gttCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTTTTAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAATTCACTCTTCACCTTGAACatat                                                                       <  2:162997/250‑1 (MQ=255)
                                                                                                                                                                          tCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAatt                                                                                            >  2:20854/1‑222 (MQ=255)
                                                                                                                                                                          tCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAatt                                                                                            <  1:20854/222‑1 (MQ=255)
                                                                                                                                                                           ccAGACAGCCCTCAAAGCGAGAATTAATCTGGACGAGAATAAGGCCGGACTGATGCCCTTTCTCCAGCCAGCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTGCTGAAATTGCCGAATCGTTGATTGTCATAATTCATTCTTCACTTTGAACATATTTAAATCt                                                               <  2:421705/250‑1 (MQ=255)
                                                                                                                                                                                cAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCGTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGAGTGTCATAATTCATGCTTCACTTTGAACATATTTAAATCTTTAATg                                                         <  1:181420/251‑1 (MQ=255)
                                                                                                                                                                                  gCCTTCAAATCGAGAATTAATCTGGACGCGAATAATGCCTGACTGATGACCTTTCTCCAGCAATAGAGACAGTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCGCGTTCTGTGTCAGCCCGTCGTTATAGTAAAACCACTCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTTTGAAATTGCCGAATTTTTGATTGTCAGAATTCATTTTTCACTTTTAACATATTTAAATCCTTAATGCa                                                       <  2:248278/251‑1 (MQ=255)
                                                                                                                                                                                                 aTTAAACTGTAGGCAAATAATGCCTGACTGATGCCCTTTCCCCAGCAATCTAGCCATTTTCAACCGTTTCAGGCCAAGACGATCGCTGATCTCTCTCTGGGTTAGCCCGTCGTGATAGTAAAACCACGCGATCCACGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAAGCGTTGATTGTCATAATTCATTTTTCACTTGTATCATATTTAAATCTTTAATGCAATTTTTCATTTCTTg                                        <  2:54134/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ccAGCAATCGGGACATTTTCAAACGTGGCAGGCCGAGACGATCGCTGATCTCGCTCTGGGGCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAATTCATTCTTCACTTTGAACATAGTTAAATCTTTAATGCAATTGTTCAGTTTTTGCTCATTTATATCGGTTAGGGCAACCGCAGTGTGGctct  <  2:137626/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
TCATGGTTGCCTCGCCAAAACCAATCGCCAGCATCTGTTGTGGTTGAAGTAAACTCATCAACATATGCGCCGCGCCTATCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAATTCATTCTTCACTTTGAACATATTTAAATCTTTAATGCAATTGTTCAGTTCTTGCTCATTTATATCTGTGATGGCAACCACAGTTTGACTCT  >  NC_000913/1600854‑1601335

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