Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,9170AC50.0% 51.0 / 17.9 50V43G (GTG→GGG) ssuBaliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (19/6);  new base C (0/25);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 1.16e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA  >  NC_000913/993772‑994060
                                                                                                                                                        |                                                                                                                                               
cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCaccacc                                                                                                                                            <  2:19710/151‑1 (MQ=255)
           ccTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                 >  1:716137/1‑151 (MQ=255)
                       ttCGCGTTTCTCCTTAATTTCCCGCCGGCGGTGTGGTCCCCGCAAACACATCCCGGGCGGTTGGCTTTCCCAACCCTGCCAGCAGGCGCAGCAGGGTTTTTTTCCCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                         <  2:828300/147‑1 (MQ=255)
                            gtttctcCCTAATTTTCAGCCAGGGGTGGGGTGCCCGCAAACACATCGCCTGCGTTGGGCGTTCCCAGCCCTGCCAGCAGGCGCAGCGGGTTACTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                 <  2:101775/147‑1 (MQ=255)
                            gtgtCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                 >  2:400890/1‑150 (MQ=255)
                              tgtttCCCGAATTTTCGCCCGGGGTTGGGGTGCCCGCAAACACCCCGCCGGGGTTGGGCTTTCCCAGCCCGGCCGGCAGGCGCAGCAGGTTCCTTTTCCCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                   <  1:373524/144‑1 (MQ=255)
                                cTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                            >  1:774042/1‑151 (MQ=255)
                                   ccTGATTTCCAGCCGGGGGGGGGGGCCCCGCAACCACCCCGCCTGGGTTGGGGTTTCCCAGCCCGCCCAGAGGGCGCGGGGGGGTCTTTTGCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                   <  1:499874/141‑1 (MQ=255)
                                   ccTGATTTCCACCCAGGGGTGTGGTCCCCGCTAACACACCCCCGGGGGTTGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCAGGTTCCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAg                                                                                                         <  2:306828/151‑1 (MQ=255)
                                   ccTGATTTCCACCCAGCGGTGGGGTCCCCGCAACCACCGCCCGGGGGTTTGGCTTTTCCAGCCCTGCCAGAGGGCGCAGGAGGGTTCTTTTCCCCCCACCGCTGCGCCCCACCTCCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                            <  2:204026/148‑1 (MQ=255)
                                     taattttCAGCGAGCGGGGTGGTCCCCGAAAACACATCGCGGGGGGTGGGCTTTCCCAGCCCTCCCAGAAGGCGCAGGGGGTTTCTTTTGCCACCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                          <  1:547616/144‑1 (MQ=255)
                                       attttCACCCAGGGGTGTGGTCCCCGCAACCACACCGCCGGGGGTTGGCTTTCCCAGCCCTGCCAGCGGGCGCAGAGGGTTCCTTTGCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                           <  2:458013/143‑1 (MQ=255)
                                         tttCGGCCGGGGGTGGGTCGCCCACAAACACACCGCTTGGGGTGGGCTTTCCAGGCCCTCCGGCCAGGCGAGGTGGGTTCCTTTTGCCACACCGGCGGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                       <  1:309657/147‑1 (MQ=255)
                                         tttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                   >  1:384035/1‑151 (MQ=255)
                                          ttccagccgggGGTGTGGTCCCCGAAAACACAGCGCCTGGGGTGGGCTTTCCCAGCCCTCCCGGAAGGCGCAGCAGGTTTCTTTTGCCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                      <  1:351808/145‑1 (MQ=255)
                                                ccgggggTGGGGTCCCCGCAAACACACCGCCTGCGGTTGGCTTTCCAAGCCCTGCCAGCAGGCGCAGCAGGTTACTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGa                                                                                              <  2:87522/145‑1 (MQ=255)
                                                 caggGGTGTGGTGCCCGTAAACACACCGCCGGGGGTTGGCCTTTCCACCCTGGCCGGCAGGCCCGGCGGGTTCCTTTCGCCCCCCCCGCCGCCGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                            <  2:382676/151‑1 (MQ=255)
                                                          gggtcCCCGCAAACACACCTCCTGGGTTGGGCTTTTCCAACCCTCCCAGCAGGCCCAGCGGGTTTCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                            <  2:703920/142‑1 (MQ=255)
                                                           ggtcCCCCCTAACCCGCCGCCGGTGGTTGGCTTTCCCAGCCCGGCCAGCGGGCGCAGCGGGGTTTTTTTGCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                   <  1:484043/151‑1 (MQ=255)
                                                              gCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                                >  1:613852/1‑151 (MQ=255)
                                                               cccGCTAACACATCGCTTGCGTTGGGGTTTTCGGGCCCTGCGAGCTGGCGGAGCGGGTTGTTTTCCCCCCCCCCGTGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                                 <  2:670113/149‑1 (MQ=255)
                                                               cccGCAACCACTCCGCCTGCGTTTGCCTTTCCCAGCCCTCCCAGCGGGCGCGGCGGGGTATTTTTCCCACCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                                <  1:175622/150‑1 (MQ=255)
                                                                ccGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                               >  1:727573/1‑150 (MQ=255)
                                                                 cGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa                                                                              >  1:8976/1‑150 (MQ=255)
                                                                   cTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                           >  2:456380/1‑151 (MQ=255)
                                                                       cacaTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                       >  1:260445/1‑151 (MQ=255)
                                                                        cccggcccTGGGGTTTGGCTTTCCAGGCCCTGCCAGCAGGCGCAGAGGGTTTTTTTCGCCACCCCCGCCGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                          <  1:313636/142‑1 (MQ=255)
                                                                                   ggtgggtgTTCCAAGCCCTGCCAGCGGGCGCAGCAGGGTACTTTGGCCACCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                           <  2:603078/151‑1 (MQ=255)
                                                                                    gTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                           >  2:758978/1‑150 (MQ=255)
                                                                                       ggCTTTCCCAGCCCTGCCAGCAGGCGCGGGGGGTTTTTTTGGCCCCCCCCGCTGCGCCCCACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                           <  1:587880/147‑1 (MQ=255)
                                                                                       ggCTTTCCCAGCCCTGCCAGCAGGCGCAGCGGGGTTCTTTCGCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                         <  2:599983/149‑1 (MQ=255)
                                                                                        gCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                       >  1:402133/1‑150 (MQ=255)
                                                                                         ggctttcccacccccgCCAGCGGGCGCAGCAGGTTCCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                        <  2:326081/135‑1 (MQ=255)
                                                                                          gTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                     >  1:489921/1‑150 (MQ=255)
                                                                                           tttCCGGCCCGGCCAGCAGGCGCGGCGGGGTACTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                     <  2:702106/149‑1 (MQ=255)
                                                                                               ccAGCCCGGCCAGCAGGCGCAGCGGGTTACTTTTCCCCCCACCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                <  2:644085/151‑1 (MQ=255)
                                                                                                 aGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                              >  2:744994/1‑151 (MQ=255)
                                                                                                    aCCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                            >  2:647227/1‑151 (MQ=255)
                                                                                                    aCCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                            >  1:96581/1‑151 (MQ=255)
                                                                                                               agGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                   >  1:377923/1‑150 (MQ=255)
                                                                                                                 gcgcAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                   >  2:763913/1‑148 (MQ=255)
                                                                                                                          ggTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                        >  1:845559/1‑150 (MQ=255)
                                                                                                                           gtgttttttgtccacccCGGCCGCCGCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                          <  2:446871/138‑1 (MQ=255)
                                                                                                                           gtgtttttCCCACCCCCTCTGGCACCCACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAg                         <  1:340917/144‑1 (MQ=255)
                                                                                                                           ggttattttccccaccgCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                          <  1:188855/138‑1 (MQ=255)
                                                                                                                             ttcttttcccccccaccGTTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                             <  1:766132/132‑1 (MQ=255)
                                                                                                                                      ccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATa                 <  1:327836/141‑1 (MQ=255)
                                                                                                                                       ccccaccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg            <  2:796320/149‑1 (MQ=255)
                                                                                                                                       caccacCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg            >  2:464224/1‑151 (MQ=255)
                                                                                                                                         cccccGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg            <  2:535049/146‑1 (MQ=255)
                                                                                                                                         ccaccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg           <  2:300310/150‑1 (MQ=255)
                                                                                                                                                     ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct     <  2:232022/144‑1 (MQ=255)
                                                                                                                                                     ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa  <  2:855553/147‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                               
CAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA  >  NC_000913/993772‑994060

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