Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,037,3330T.20.4% 136.6 / 31.8 49coding (820/1359 nt)yedVputative sensory kinase in two‑component regulatory system with YedW
Reads supporting (aligned to +/- strand):  ref base T (17/22);  new base . (4/6);  total (21/28)
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.59e-01

CAAAACATTGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCGAGATCGT  >  NC_000913/2037241‑2037431
                                                                                             |                                                                                                  
cAAAACATTGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTg                                                                                              <  2:543521/99‑1 (MQ=255)
     cATTGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTc                                                                                        >  2:355574/1‑100 (MQ=255)
       ttGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAg                                                                                      >  1:335912/1‑100 (MQ=255)
       ttGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAg                                                                                      >  1:280905/1‑100 (MQ=255)
        tGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGTCAATTGTTTTTGATATTCAGCg                                                                                    <  2:757258/100‑1 (MQ=255)
        tGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTTGATATTCAg                                                                                      <  1:389260/100‑1 (MQ=255)
        tGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGc                                                                                     <  2:656943/100‑1 (MQ=255)
         gTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCg                                                                                    >  1:11871/1‑100 (MQ=255)
                 cTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCt                                                                            >  1:410359/1‑100 (MQ=255)
                 cTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCt                                                                            >  2:1083868/1‑100 (MQ=255)
                    cccTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTGATATTCAGCGATACTTCTGGtt                                                                        >  1:1155671/1‑100 (MQ=255)
                    cccTGGCAAGAAACAGTATATTCTATGTTAAACGCAGAATAGTTTCCAGATCTTCAATGTTTCCGTCAACTGTTTCTTGATATTCAGCGATACttat                                                                            >  2:648855/1‑95 (MQ=255)
                     ccTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGtt                                                                        <  1:814021/100‑1 (MQ=255)
                     ccTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGt                                                                         <  2:1027956/99‑1 (MQ=255)
                       tGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTg                                                                      <  2:125141/100‑1 (MQ=255)
                        ggCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGa                                                                     >  2:237372/1‑100 (MQ=255)
                         gCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGAc                                                                    <  2:363598/100‑1 (MQ=255)
                            aGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTg                                                                  >  1:635721/1‑99 (MQ=255)
                            aGAAACAGTATGTTCTCTGTTAACCGCGAAATATCTCCGAGCTCGTCAATGTTTCCGGCAATTGTTTTTTTATATTCAGCGACACTTCTAGtcaggct                                                                   >  2:1293258/1‑92 (MQ=255)
                               aaCAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCg                                                              <  1:374997/100‑1 (MQ=255)
                               aaCAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCg                                                              <  2:141428/100‑1 (MQ=255)
                                  aGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAAc                                                          <  2:491875/100‑1 (MQ=255)
                                   gTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAAcc                                                         <  1:408769/100‑1 (MQ=255)
                                    tATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAAcc                                                         <  1:984819/100‑1 (MQ=255)
                                       gTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGa                                                      >  2:1032671/1‑100 (MQ=255)
                                         tctcTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGAtt                                                    >  1:670975/1‑100 (MQ=255)
                                          ctctGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTTTTTTTTTTCTATACTTCTGGTTTGACTGAGCGTAACCTGATTc                                                   >  2:874694/1‑100 (MQ=255)
                                          ctctGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTc                                                   <  2:997627/100‑1 (MQ=255)
                                          ctctGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTc                                                   <  1:90680/100‑1 (MQ=255)
                                              gTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGAc                                               <  2:1080673/100‑1 (MQ=255)
                                                 aaCCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAg                                           <  2:904846/100‑1 (MQ=255)
                                                  aCCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAg                                           <  2:145250/100‑1 (MQ=255)
                                                  aCCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAg                                           <  1:1092702/100‑1 (MQ=255)
                                                   ccGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTa                                         >  1:897081/1‑100 (MQ=255)
                                                            atTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTaa                                >  2:954719/1‑100 (MQ=255)
                                                              tttCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAAt                               <  1:982732/100‑1 (MQ=255)
                                                               ttCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAAtt                              <  1:354882/100‑1 (MQ=255)
                                                                tCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTgg                            <  1:33856/100‑1 (MQ=255)
                                                                  cAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGc                           <  1:286323/100‑1 (MQ=255)
                                                                   aGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCg                          <  2:635717/100‑1 (MQ=255)
                                                                             tGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTc                >  1:1064665/1‑100 (MQ=255)
                                                                                    ggCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCg         >  1:233166/1‑100 (MQ=255)
                                                                                    ggCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCg         <  1:1032680/100‑1 (MQ=255)
                                                                                    ggCAATTGTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCga        >  1:470102/1‑100 (MQ=255)
                                                                                     gCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCga        <  1:355577/100‑1 (MQ=255)
                                                                                     gCAATTGTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCgag       <  1:971919/100‑1 (MQ=255)
                                                                                      cAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCgag       <  2:747440/100‑1 (MQ=255)
                                                                                      cAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCgag       >  2:74987/1‑100 (MQ=255)
                                                                                           gTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCGAGATCGt  >  2:334318/1‑100 (MQ=255)
                                                                                             |                                                                                                  
CAAAACATTGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCGAGATCGT  >  NC_000913/2037241‑2037431

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