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

CAAAACATTGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCGAGATCGT  >  NC_000913/2037241‑2037431
                                                                                             |                                                                                                  
cAAAACATTGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTg                                                                                              <  2:543521/99‑1 (MQ=255)
     cATTGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTc                                                                                        >  2:355574/1‑100 (MQ=255)
       ttGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAg                                                                                      >  1:335912/1‑100 (MQ=255)
       ttGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAg                                                                                      >  1:280905/1‑100 (MQ=255)
        tGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGTCAATTG‑‑TTTTTGATATTCAGCg                                                                                    <  2:757258/100‑1 (MQ=255)
        tGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTGTTTTTTTGATATTCAg                                                                                      <  1:389260/100‑1 (MQ=255)
        tGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGc                                                                                     <  2:656943/100‑1 (MQ=255)
         gTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCg                                                                                    >  1:11871/1‑100 (MQ=255)
                 cTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCt                                                                            >  1:410359/1‑100 (MQ=255)
                 cTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCt                                                                            >  2:1083868/1‑100 (MQ=255)
                    cccTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑‑TTTTTGATATTCAGCGATACTTCTGGtt                                                                        >  1:1155671/1‑100 (MQ=255)
                    cccTGGCAAGAAACAGTATATTCTATGTTAAACGCAGAATAGTTTCCAGATCTTCAATGTTTCCGTCAACTG‑TTTCTTGATATTCAGCGATACttat                                                                            >  2:648855/1‑95 (MQ=255)
                     ccTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGtt                                                                        <  1:814021/100‑1 (MQ=255)
                     ccTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGt                                                                         <  2:1027956/99‑1 (MQ=255)
                       tGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTg                                                                      <  2:125141/100‑1 (MQ=255)
                        ggCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGa                                                                     >  2:237372/1‑100 (MQ=255)
                         gCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGAc                                                                    <  2:363598/100‑1 (MQ=255)
                            aGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTg                                                                  >  1:635721/1‑99 (MQ=255)
                            aGAAACAGTATGTTCTCTGTTAACCGCGAAATATCTCCGAGCTCGTCAATGTTTCCGGCAATTG‑TTTTTTTATATTCAGCGACACTTCTAGtcaggct                                                                   >  2:1293258/1‑92 (MQ=255)
                               aaCAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCg                                                              <  1:374997/100‑1 (MQ=255)
                               aaCAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCg                                                              <  2:141428/100‑1 (MQ=255)
                                  aGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑‑TTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAAc                                                          <  2:491875/100‑1 (MQ=255)
                                   gTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑‑TTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAAcc                                                         <  1:408769/100‑1 (MQ=255)
                                    tATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAAcc                                                         <  1:984819/100‑1 (MQ=255)
                                       gTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGa                                                      >  2:1032671/1‑100 (MQ=255)
                                         tctcTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGAtt                                                    >  1:670975/1‑100 (MQ=255)
                                          ctctGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTTTTTTTTTTCTATACTTCTGGTTTGACTGAGCGTAACCTGATTc                                                   >  2:874694/1‑100 (MQ=255)
                                          ctctGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTc                                                   <  2:997627/100‑1 (MQ=255)
                                          ctctGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTc                                                   <  1:90680/100‑1 (MQ=255)
                                              gTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGAc                                               <  2:1080673/100‑1 (MQ=255)
                                                 aaCCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑‑TTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAg                                           <  2:904846/100‑1 (MQ=255)
                                                  aCCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAg                                           <  2:145250/100‑1 (MQ=255)
                                                  aCCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAg                                           <  1:1092702/100‑1 (MQ=255)
                                                   ccGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑‑TTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTa                                         >  1:897081/1‑100 (MQ=255)
                                                            atTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑‑TTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTaa                                >  2:954719/1‑100 (MQ=255)
                                                              tttCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAAt                               <  1:982732/100‑1 (MQ=255)
                                                               ttCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAAtt                              <  1:354882/100‑1 (MQ=255)
                                                                tCCAGCTCTTCAATGTTTCCGGCAATTG‑‑TTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTgg                            <  1:33856/100‑1 (MQ=255)
                                                                  cAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGc                           <  1:286323/100‑1 (MQ=255)
                                                                   aGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCg                          <  2:635717/100‑1 (MQ=255)
                                                                             tGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTc                >  1:1064665/1‑100 (MQ=255)
                                                                                    ggCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCg         >  1:233166/1‑100 (MQ=255)
                                                                                    ggCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCg         <  1:1032680/100‑1 (MQ=255)
                                                                                    ggCAATTG‑‑TTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCga        >  1:470102/1‑100 (MQ=255)
                                                                                     gCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCga        <  1:355577/100‑1 (MQ=255)
                                                                                     gCAATTG‑‑TTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCgag       <  1:971919/100‑1 (MQ=255)
                                                                                      cAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCgag       <  2:747440/100‑1 (MQ=255)
                                                                                      cAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCgag       >  2:74987/1‑100 (MQ=255)
                                                                                           g‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCGAGATCGt  >  2:334318/1‑100 (MQ=255)
                                                                                             |                                                                                                  
CAAAACATTGTTTTTATCTGCCCTGGCAAGAAACAGTATGTTCTCTGTTAACCGCGAAATATTTTCCAGCTCTTCAATGTTTCCGGCAATTG‑TTTTTTGATATTCAGCGATACTTCTGGTTTGACTGAGCGTAACCTGATTCTGACCCAGTAATGCATTAATTGGCGTTCTAAGTTCATGAGCGAGATCGT  >  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: ‑