Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 335,166 T→A 8.9% intergenic (‑144/‑114) yahC ← / → yahD putative inner membrane protein/ankyrin repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913335,1660TA 8.9% 152.6 / 5.5 45intergenic (‑144/‑114)yahC/yahDputative inner membrane protein/ankyrin repeat protein
Reads supporting (aligned to +/- strand):  ref base T (19/22);  new base A (2/2);  total (21/24)
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.83e-01

CCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATT  >  NC_000913/335030‑335310
                                                                                                                                        |                                                                                                                                                
cccTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACaaa                                                                                                                                       <  1:695704/148‑1 (MQ=255)
   tctGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAg                                                                                                                                    <  1:444556/148‑1 (MQ=255)
                   cGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGc                                                                                                                                  >  4:674894/1‑134 (MQ=255)
                   cGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGc                                                                                                                                  <  3:674894/134‑1 (MQ=255)
                      aaGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTaa                                                                                                                 <  4:841757/148‑1 (MQ=255)
                        gCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTca                                                                                                                                            >  4:1068900/1‑119 (MQ=255)
                        gCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTca                                                                                                                                            <  3:1068900/119‑1 (MQ=255)
                          tACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAAtatt                                                                                                                                                 <  4:946561/112‑1 (MQ=255)
                          tACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAAtatt                                                                                                                                                 >  3:946561/1‑112 (MQ=255)
                                tACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAAtattat                                                                                                                                               <  2:304232/108‑1 (MQ=255)
                                tACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAAtattat                                                                                                                                               >  1:304232/1‑108 (MQ=255)
                                    tgtgCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTc                                                                                                   >  1:341090/1‑148 (MQ=255)
                                               taGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGAt                                                                                        <  2:341090/148‑1 (MQ=255)
                                                   gaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTaa                                                                                                           >  4:478881/1‑125 (MQ=255)
                                                   gaAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTaa                                                                                                           <  3:478881/125‑1 (MQ=255)
                                                        aaaaCCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGc                                                                               >  4:130944/1‑148 (MQ=255)
                                                           aCCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTca                                                                                                                                            <  1:216117/84‑1 (MQ=255)
                                                           aCCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATAATATTca                                                                                                                                            >  2:216117/1‑84 (MQ=255)
                                                                      tGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAatat                                                                                                                                                  <  1:790586/67‑1 (MQ=255)
                                                                      tGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAatat                                                                                                                                                  >  2:790586/1‑67 (MQ=255)
                                                                      tGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGt                                                                                  >  1:1064266/1‑131 (MQ=255)
                                                                      tGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGt                                                                                  <  2:1064266/131‑1 (MQ=255)
                                                                      tGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCAt                                                                 >  4:421486/1‑148 (MQ=255)
                                                                        aaaGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCtttt                                                                                               <  1:605667/116‑1 (MQ=255)
                                                                        aaaGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCtttt                                                                                               >  2:605667/1‑116 (MQ=255)
                                                                           gggTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGaaa                                                            <  3:130944/148‑1 (MQ=255)
                                                                               aTCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAg                                                        <  3:421486/148‑1 (MQ=255)
                                                                                     tCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACaaaa                                                                                                                                      >  4:898890/1‑64 (MQ=255)
                                                                                     tCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACaaaa                                                                                                                                      <  3:898890/64‑1 (MQ=255)
                                                                                             aCAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTaa                                                                                                                 <  4:202696/77‑1 (MQ=255)
                                                                                             aCAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTaa                                                                                                                 >  3:202696/1‑77 (MQ=255)
                                                                                              cAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAg                                                        >  4:1140010/1‑133 (MQ=255)
                                                                                              cAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAg                                                        <  3:1140010/133‑1 (MQ=255)
                                                                                                aaTATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCAt                                                                 <  2:14617/122‑1 (MQ=255)
                                                                                                aaTATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCAt                                                                 >  1:14617/1‑122 (MQ=255)
                                                                                                    ttGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGAtt                                                                                       >  4:1119615/1‑96 (MQ=255)
                                                                                                    ttGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGAtt                                                                                       <  3:1119615/96‑1 (MQ=255)
                                                                                                           aCCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAg                                                                                                          <  3:846052/70‑1 (MQ=255)
                                                                                                           aCCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAg                                                                                                          >  4:846052/1‑70 (MQ=255)
                                                                                                           aCCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGa                                 >  1:124778/1‑143 (MQ=255)
                                                                                                           aCCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGa                                 <  2:124778/143‑1 (MQ=255)
                                                                                                                        aaaaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCt                                                                                                  <  1:1074522/65‑1 (MQ=255)
                                                                                                                        aaaaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCt                                                                                                  >  2:1074522/1‑65 (MQ=255)
                                                                                                                         aaaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAAc                                              <  3:555080/116‑1 (MQ=255)
                                                                                                                         aaaGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAAc                                              >  4:555080/1‑116 (MQ=255)
                                                                                                                             aTTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGACTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAg                                           <  1:48653/115‑1 (MQ=255)
                                                                                                                             aTTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGACTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAg                                           >  2:48653/1‑115 (MQ=255)
                                                                                                                                  tgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCtttt                                                                                               <  4:1029234/58‑1 (MQ=255)
                                                                                                                                  tgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCtttt                                                                                               >  3:1029234/1‑58 (MQ=255)
                                                                                                                                  tgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCa                                            <  1:212204/109‑1 (MQ=255)
                                                                                                                                  tgAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCa                                            >  2:212204/1‑109 (MQ=255)
                                                                                                                                     atatTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAAcc                                             >  1:752584/1‑105 (MQ=255)
                                                                                                                                     atatTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAAcc                                             <  2:752584/105‑1 (MQ=255)
                                                                                                                                     atatTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACtatt  >  4:381279/1‑148 (MQ=255)
                                                                                                                                        |                                                                                                                                                
CCCTCTGGTTAGACATGGACGAAAGCTACTCTTACATGTGCAGACTATAGAGAAATAAAACCTGATTGAGTGAAAGGGTATCGGGTCAAAGAAACAAATATTGCACTACCGCACACTGCGAAAAGATTGTTGAATATTATTCAACAAAAAGGCGAGATACGGCGAATTAAATTTAAGCGGGTTCTTTTATCAGATTATCGTTGCTATTCTCAATTCATCGAAAACAGCCACAAAAACCAGGGGAATGTGAATGAGTATAAAAAATCTACCTGCCGACTATT  >  NC_000913/335030‑335310

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