Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 236,911 T→C intergenic (+113/‑20) dnaQ → / → aspV DNA polymerase III subunit epsilon/tRNA‑Asp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913236,9110TC98.1% 185.2 / ‑4.4 53intergenic (+113/‑20)dnaQ/aspVDNA polymerase III subunit epsilon/tRNA‑Asp
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (24/28);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GGAAGTTGCCTCTGGCGAGCATAAATACCTGTGAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGTTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTCACTCATGAAAATGAGTTCAGAGAGCCGCAAGAT  >  NC_000913/236775‑237046
                                                                                                                                        |                                                                                                                                       
ggAAGTTGCCTCTGGCGAGCATAAATACCTGTGAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTcc                                                                                                                                      <  2:235642/140‑1 (MQ=255)
ggAAGTTGCCTCTGGCGAGCATAAATACCTGTGAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTcc                                                                                                                                      <  1:510484/140‑1 (MQ=255)
 gAAGTTGCCTCTGGCGAGCATAAATACCTGTGAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTccc                                                                                                                                     >  2:116676/1‑140 (MQ=255)
 gAAGTTGCCTCTGGCGAGCATAAATACCTGTGAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTccc                                                                                                                                     >  1:866493/1‑140 (MQ=255)
 gAAGTTGCCTCTGGCGAGCATAAATACCTGTGAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTccc                                                                                                                                     >  2:311195/1‑140 (MQ=255)
  aaGTTGCCTCTGGCGAGCATAAATACCTGTGAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCa                                                                                                                                    <  1:885256/140‑1 (MQ=255)
        ccTCTGGCGAGCATAAATACCTGTGAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGaa                                                                                                                              >  2:668211/1‑140 (MQ=255)
                   cATAAATACCTGTGAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGa                                                                                                                   >  1:439971/1‑140 (MQ=255)
                     tAAATACCTGTGAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGc                                                                                                                 <  2:828179/140‑1 (MQ=255)
                             tgtgAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTc                                                                                                         >  1:600411/1‑140 (MQ=255)
                             tgtgAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTc                                                                                                         >  1:996540/1‑140 (MQ=255)
                              gtgAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCa                                                                                                        >  2:788881/1‑140 (MQ=255)
                               tgAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAg                                                                                                       >  2:579807/1‑140 (MQ=255)
                               tgAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAg                                                                                                       <  1:997417/140‑1 (MQ=255)
                                 aaaGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGTTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTc                                                                                                     <  1:200448/140‑1 (MQ=255)
                                 aaaGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTc                                                                                                     <  2:242134/140‑1 (MQ=255)
                                 aaaGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTc                                                                                                     <  1:527411/140‑1 (MQ=255)
                                    ggCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGt                                                                                                  <  1:766537/140‑1 (MQ=255)
                                    ggCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGt                                                                                                  >  2:489616/1‑140 (MQ=255)
                                      cgcTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAGCCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTa                                                                                                >  2:729685/1‑140 (MQ=255)
                                      cgcTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTa                                                                                                <  2:742125/140‑1 (MQ=255)
                                      cgcTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTa                                                                                                >  1:976665/1‑140 (MQ=255)
                                                 gCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATAcctgcc                                                                                     <  2:199735/140‑1 (MQ=255)
                                                  cGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATAcctgcct                                                                                    <  1:520599/140‑1 (MQ=255)
                                                    aCTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGt                                                                                  <  2:187787/140‑1 (MQ=255)
                                                          gCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCa                                                                            <  1:788881/140‑1 (MQ=255)
                                                           cGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAg                                                                           <  2:429806/140‑1 (MQ=255)
                                                             aTTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAggg                                                                         >  1:363369/1‑140 (MQ=255)
                                                              tttttGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAgggg                                                                        <  2:996540/140‑1 (MQ=255)
                                                                  tGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTcg                                                                    >  2:150695/1‑140 (MQ=255)
                                                                  tGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTcg                                                                    <  1:311195/140‑1 (MQ=255)
                                                                    cagcaAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTcgcg                                                                  <  2:866493/140‑1 (MQ=255)
                                                                        aaaCGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGtt                                                              <  2:439971/140‑1 (MQ=255)
                                                                        aaaCGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCCCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGtt                                                              <  1:401415/140‑1 (MQ=255)
                                                                                   aaGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCccgtc                                                   >  1:23637/1‑140 (MQ=255)
                                                                                   aaGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCccgtc                                                   <  2:829877/140‑1 (MQ=255)
                                                                                        gagaAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTc                                              <  2:190447/140‑1 (MQ=255)
                                                                                        gagaAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTc                                              >  2:28083/1‑140 (MQ=255)
                                                                                             aaaCCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCa                                         <  1:668211/140‑1 (MQ=255)
                                                                                                  gTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTAtt                                    <  1:245906/140‑1 (MQ=255)
                                                                                                   ttGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTc                                   >  1:855246/1‑140 (MQ=255)
                                                                                                   ttGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGGCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTc                                   >  2:977595/1‑140 (MQ=255)
                                                                                                    ggACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACATCGCGGAGCGGTAGTTCAGTCGGTTAGAATCCCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTCa                                  <  1:116676/139‑1 (MQ=255)
                                                                                                     gACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTCAc                                 <  2:240869/140‑1 (MQ=255)
                                                                                                              cGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGTGTTCGAGTCCCGTCCGTTCCGCCACTATTCACTCATGaaaa                        <  1:721552/140‑1 (MQ=255)
                                                                                                              cGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTCACTCATGaaaa                        <  1:873616/140‑1 (MQ=255)
                                                                                                              cGAGGCAATCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTCACTCATGaaaa                        <  1:866930/140‑1 (MQ=255)
                                                                                                                      tCCGTAATATTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTCACTCATGAAAATGAGTTCa                >  2:947492/1‑140 (MQ=255)
                                                                                                                            atatTCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTCACTCATGAAAATGAGTTCAGAGAGc          >  1:749088/1‑140 (MQ=255)
                                                                                                                                tCGCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTCACTCATGAAAATGAGTTCAGAGAGCCGCa      >  1:41499/1‑140 (MQ=255)
                                                                                                                                  gCCTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTCACTCATGAAAATGAGTTCAGAGAGCCGCAAg    >  2:557090/1‑140 (MQ=255)
                                                                                                                                    cTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTCACTCATGAAAATGAGTTCAGAGAGCCGCAAGAt  >  1:379981/1‑140 (MQ=255)
                                                                                                                                    cTCGCTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTCACTCATGAAAATGAGTTCAGAGAGCCGCAAGAt  >  1:138030/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GGAAGTTGCCTCTGGCGAGCATAAATACCTGTGAAAGGCGCTAAAAATAGCGACTTGGGCGATTTTTGCAGCAAACGATTCAAAAGATGAGAAAAACCGTTGACGAAGGTCGAGGCAATCCGTAATATTCGCCTCGTTCCCAACGGAACACAACGCGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACTATTCACTCATGAAAATGAGTTCAGAGAGCCGCAAGAT  >  NC_000913/236775‑237046

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