Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,647,949 C→T intergenic (‑161/‑114) dinQ ← / → agrA membrane toxin DinQ/small RNA AgrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,647,9490CT95.5% 57.2 / ‑4.3 23intergenic (‑161/‑114)dinQ/agrAmembrane toxin DinQ/small RNA AgrA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (9/12);  total (9/14)
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.33e-01

CTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAG  >  NC_000913/3647702‑3648191
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
cTATCAGTTCAGCAGCTGACGCCGAAATCGGATCAGTTCCCGCAGCGCAATTAACTACCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACCGCCTGTCACGTCGCCTGTCTCATTGTTAGATAAGAACTCTCTCACTCCTGCCAGAGCATCCGTTCAAGGCATCACACGTACCTCTGCCCAGGCCTTTGATAGCGTTTTTgctct                                                                                                                                                                                                                                                 >  2:3710/1‑248 (MQ=255)
   ccAGTTCAGAAGCGGAAGCAGAAAGCGGCTCAGTCCCATCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCAAGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCGGTCACGGCGCATGTCTCATTGTTAGAGAAGAACTCTCTCACTCTGGCCAGAGCATCAGTTAACGGCACCACCCGTCCTTCTGACCAGGGCTTTGAAAGCGTTTATGCGCATTg                                                                                                                                                                                                                                              <  2:151278/250‑1 (MQ=255)
                            cGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTa                                                                                                                                                                                                                     <  1:336124/251‑1 (MQ=255)
                              gATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATCAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTAAGAGCTAAc                                                                                                                                                                                                                   >  1:93187/1‑251 (MQ=255)
                                      ccAGCAGCGCAATTAACGCCCCTAGAACGATGATTGTTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATaa                                                                                                                                                                                                            <  1:245249/250‑1 (MQ=255)
                                         gcagcGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCAt                                                                                                                                                                                                        <  1:325376/251‑1 (MQ=255)
                                            gcgcAATTAAAGCCCCTATAACAATGATTGCTTTATCAATCACAAGATTTCTCCATGCGATGGAGAGAGAATGCATCCGCTTACTCATCCCATGCCTGTAACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGACCCACCCTCACTTCTGACCCGGACTTTGAAAGAGTTTATGTGCATTGTCAGATTGACAATGCCACAGATCTAACTCAATAATCGTAc                                                                                                                                                                                                      >  2:70614/1‑250 (MQ=255)
                                                aaTTTACGACCCTAGAACGATGATTTCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATTCGCTTACTCATCCTCTGCCTGTCACTGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCAACACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTAc                                                                                                                                                                                                 <  2:208924/251‑1 (MQ=255)
                                                         cccTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGTATGCACCCGCTTACGAATCCACTGCTTGTCACGGCGCATGTCTCATGGTTAGATAAGAACTCTCTCTCTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATGGCCAGAGTGAAAATGCCTCAGAGCTAACTGGAGAATCATACATTACATGCAgggt                                                                                                                                                                                        <  2:152097/251‑3 (MQ=255)
                                                                   atgatTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAAcc                                                                                                                                                                                >  1:254270/1‑249 (MQ=255)
                                                                                                ccATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATg                                                                                                                                                    <  1:301452/248‑1 (MQ=255)
                                                                                                ccATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCGGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATg                                                                                                                                                    <  1:141224/248‑1 (MQ=255)
                                                                                                        tGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCAtt                                                                                                                                         >  1:204684/1‑251 (MQ=255)
                                                                                                                      atccggatccTCATCCACTGCCGGTTACGGCGTATGTCTCATTGTTAGATAAGAACTCTCTCAATCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGCTAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCGTGCAATAGTTTCAGTATGGTATTTGCATTGATGCGTTAgatg                                                                                                                            <  2:277497/241‑1 (MQ=255)
                                                                                                                             tacgcacccACTGACTGTCACGGCGCATGTCTGAATGTTATAGAAGAACGCTCGCACTCCGGCCAGAGCATCAGTTAACGGCACCACACGTACTTCGGACCAGGACTTTGAAAGCGTGTATGTTTATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCATCACGTCCTTGCAAGAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATTATGGCTa                                                                                                                     <  2:322621/243‑1 (MQ=255)
                                                                                                                                       ccgtcTGTAACGTCGCATGTCTCAATATTAGAAAAGATCTTTCTCAATCCGGCCAGAGCATCAGTTAACGGCACCACCCTGACTTCTGACCAGGACTTTGAAAGCGTTTACGTGCATTGCCAGAGTGAAAATGACTAAGAGCTAACTGGATAATCAGACAGTACATGCAGGTTTTAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGAAGATGGCTATCTCTCTCCat                                                                                                          <  2:110647/247‑2 (MQ=255)
                                                                                                                                           ccgTCACCGCGCATGTCTCATTGTTCGCTAAGAACTCTATCACTCCGGCCAGCGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGATTATGTGCATTGCCAGATTAAAAATGCCTCAGAGCTAACTGGCTAATCATACAGTACATTCCTGTTATAACACCAGCTCGACCTTGCAATAGTTTCAGAATGGAATTAGCATTGATGAGTTAGATGATTGCAATCTAACTCCAgct                                                                                                        >  2:248939/3‑247 (MQ=255)
                                                                                                                                                                   tAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAAc                                                                                                                                                                                 >  1:94868/1‑152 (MQ=255)
                                                                                                                                                                   tAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAAc                                                                                                                                                                                 <  2:94868/152‑1 (MQ=255)
                                                                                                                                                                        aataaCTCTCGCTCTCCTGCCAGAGCATCAGTTCACGGCACCACCCGTACTTCTTACCAGGACGTGGAAAGCGGTTATGTGCATTGCCAGAGTGAAAATTCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTTCTATAGTTTCATTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGATCTCCCAACTCAGGGCTGGAAAGTaaaa                                                                         <  2:326312/248‑1 (MQ=255)
                                                                                                                                                                                                                   cccGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCggcca                                >  2:382125/1‑246 (MQ=255)
                                                                                                                                                                                                                                ccAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                 >  1:92081/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg      >  1:23866/1‑250 (MQ=255)
                                                                                                                                                                                                                                               cGTTTATGGTCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATTCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCTACCAACTCATGGCTGGAAAGGAAAAAACCGACGCAAAGACGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTTAAAg  <  2:263568/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
CTATCAGTTCAGAAGCTGAAGCAGAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAG  >  NC_000913/3647702‑3648191

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