Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913718,6690CT43.8% 52.7 / 21.0 48E598K (GAA→AAA) speFornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base C (8/19);  new base T (21/0);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 3.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATC  >  NC_000913/718523‑718819
                                                                                                                                                  |                                                                                                                                                       
cGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTTGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGa                                                                                                                                                     <  2:542444/151‑1 (MQ=255)
  gCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGaca                                                                                                                                                   <  1:305974/151‑1 (MQ=255)
       cGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACt                                                                                                                                              <  1:256969/151‑1 (MQ=255)
         aaCCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGa                                                                                                                                            <  1:514222/151‑1 (MQ=255)
          aCCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGAc                                                                                                                                           <  2:713450/151‑1 (MQ=255)
           ccAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACg                                                                                                                                          <  2:225156/151‑1 (MQ=255)
            cAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACaaattac                                                                                                                                           >  1:449931/1‑144 (MQ=255)
                      tCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATaa                                                                                                                                <  1:478713/150‑1 (MQ=255)
                               aaaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGGTAAAAATCATTCATTTTCTTGCACCACCTACGCAGGGGATAACCgg                                                                                                                            >  2:284932/1‑145 (MQ=255)
                               aaaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAg                                                                                                                       <  1:751015/150‑1 (MQ=255)
                                       gTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAAAAAACATGCATTTTCTTGCACAACTGACGCAGGGGATAACCGGCGTAGCGctc                                                                                                                  >  2:701099/1‑147 (MQ=255)
                                         aGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCAt                                                                                                             <  1:351195/150‑1 (MQ=255)
                                         aGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATg                                                                                                            <  2:458581/151‑1 (MQ=255)
                                            tGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATTCATTTTCTGCCCCAACTGACGCAGGGGATAAACGGCGTAGCGCTCCTCATg                                                                                                            >  1:286521/1‑148 (MQ=255)
                                                   ttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTAAAAAACATGCATTTTCTTACACAACTGACGCAGGGGATAAACCGCGGAGCGCCCTTCATGCTGttttt                                                                                                   >  1:646759/1‑147 (MQ=255)
                                                    tGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATGCATTTCCTTACACAACTGACGCCGGGGATAAACGGCGTAGCGCTCTTCATTCTGGTTGt                                                                                                   >  1:423875/1‑148 (MQ=255)
                                                       ggATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACAAATGACGCAGGGGATAACCGGCGGAGCGCTCTTCATGCTGTTTTTTGATgg                                                                                             >  1:592021/1‑151 (MQ=255)
                                                        gATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGa                                                                                            <  2:754535/151‑1 (MQ=255)
                                                            cATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGt                                                                                        <  1:877407/151‑1 (MQ=255)
                                                            cATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGt                                                                                        <  2:825094/151‑1 (MQ=255)
                                                                cTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAATCATGCATTTCCTTCCACAACTGAAGCCGGGGGTAAACCGCCTAGCGCTCCTTCTGCTGTTTTTTGATGGAAGGGAg                                                                                      >  1:551825/1‑149 (MQ=255)
                                                                    cgcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCCAAAAAACAAGCGTTTCCTTCCACAAATTACGGCGGGGGTAAACGGCCGAGGGCGCCTTATGCTCTTTTGTGATAGAAGGTAGAAca                                                                                  >  2:813094/1‑148 (MQ=255)
                                                                     gcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACAACTGACGCAGGGGATAAACGGCGGAGCGCTCTTCCTGCTGTTTGTTGATGGGAGGGAGCAcct                                                                                 >  2:239345/1‑147 (MQ=255)
                                                                             aGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTTCTGACACAACAGACGCCGGGGATAAACCGCGTAGCGCTCTTCATGCTGGTTGTAGATGGGAGGGAGCACTTTTGCTTAttg                                                                       >  1:902549/1‑149 (MQ=255)
                                                                              gTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGgc                                                                      <  2:201326/151‑1 (MQ=255)
                                                                                gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTTCACAAATGACGCAGGGGATAAACGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGGAGCACTTCTGCTAATGgcg                                                                     >  1:87625/1‑150 (MQ=255)
                                                                                      tACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTTCTGACACAACTGACGCAGGGGATAAACGGCGGAGCGCTCTTTCTTCTGTTTTTAGATGGAAAGGAGGACTTCTTCTAATGGCGCGTCgg                                                               >  1:862701/1‑150 (MQ=255)
                                                                                      tACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCgg                                                               <  1:653299/150‑1 (MQ=255)
                                                                                       acggtgcATCTATTTTTGCATTTGAGTCACGTTGTGGAGAGCATACAAGTCATGCAGTTCCTGTCACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTc                                                                 <  2:633982/141‑1 (MQ=255)
                                                                                       aCGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTAAAAAATATTCATTTTCTTACACCACTGACGCAGGGGATAAACGGCGGAGCGCTCTTCTTGCTGTTTTTAGATGGGAGGGAGGACCTCTGCTTATGGGGcggcggg                                                              >  2:767964/1‑149 (MQ=255)
                                                                                        cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAACATTCATTTTCTTTCCCCACTTACGCAGGGGGTAAACGGCGTAGCGCCCTTCCTGCTGTTTTTAGATGGAAGGGAGCACTTTTGCTAATGGGGcggcgggc                                                             >  2:695120/1‑148 (MQ=255)
                                                                                         ggAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATTCATTTCCTTACACAACTGACGCAGGGGATAAACCGCGTAGCGCTCTTCATGCTGTTTGTAGATGGGAGGTAGGACTTCTGCTAATGGCGCGTCGGGCTc                                                           >  1:201988/1‑151 (MQ=255)
                                                                                             cATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATTATTCTTTTTCTGACACAAATGACGCAGGGGGTAAACGGCGGAGCGCTCTTCCTGCTGTTTGTAGATGGAAGGGAGCACTTCTTCTAATGGCGCGGCGGACTCaa                                                         >  2:826543/1‑149 (MQ=255)
                                                                                             cATCTCTTTTAGCAGTTGTTTCACGTTGTGGCGGGGCAACAAAACATGCATTTTCTTAAACAAATGAAGCAGGGGTTAAACCGCGTAGCGGCCTTCTTGATGTTTTTTGAAGGGAAGGAGGACTTCTTCTTATGGCGGGGCGGCGTCaa                                                         >  1:828752/1‑149 (MQ=255)
                                                                                                ctctTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATGCTTTTTCTTACCCAACCCACGCCGGGGATAAACAGCGGAGCGCTCTTCCTGCTCTTTGTAGATGGAAGGGAGCACCTCTGCCCATGGGGCGTCgg                                                               >  2:464286/1‑140 (MQ=255)
                                                                                                  ctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGt                                                    <  2:547818/149‑1 (MQ=255)
                                                                                                     tttGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATTGATTTTCTTACACCACTGACGCAGGGGATAACCCGCGGAGCGCTCTTTATTCTGTTTGTAGGTGGAAGGGAGGACTTCTGCTTATTGCGCGCCGGGCTTCAGCAGgttttt                                               >  2:760422/1‑150 (MQ=255)
                                                                                                       tGCAGTTGTTTCACGTTGTGGCGGGGATAAAAAACATGCATTTTCTGAAACAACAGCCGGAAGGGAAAAACCGCGTCGGGCTCTTTATGCTGGGTTGTGGAGGAAGGGAGGACTTCTGCTAAAGGCGCGTCGGGCTCAAGCAGgtttttg                                              >  1:152227/1‑148 (MQ=255)
                                                                                                           gttgttTCACGTTGTGGCGGGCCTACAAAACATGCATTTTCTTACACAACTGACGCAGGGGATAAACGGCGGAGCGCTCTTCATTCTGTTTTTAGATGGGAGGGAGCACTTCTTCTAATGGCGCGTCGGGCTCAAGCAGGTTTTTGAAGCg                                         >  2:13036/1‑151 (MQ=255)
                                                                                                           gttgttTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTGAAAAAAATAACGGCGGGGGTAAACGGCGGAGCACTCTTCTTGCTGTTTTTAGATGGGAGGGAGCACCTTTTCTTATGGGGCGTCGGACTCACACAGTTTTTTGAAGc                                          >  1:228654/1‑150 (MQ=255)
                                                                                                              gttTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTa                                       <  2:288756/150‑1 (MQ=255)
                                                                                                                ttCACGTTGTGGCGGGCATACAAAAAAATCATTTTCTGACACAAATGACGCAGGGGATAAACGGCGCAGCGCCCTTTATGCTGTTTTTAGATGGAAGGGAGCACTTCTGCTAATGGCGCGTCGGGCTCAAGCAGGTTTTTGAAGGGTAcc                                     >  1:605497/1‑150 (MQ=255)
                                                                                                                   aCGTTGTGGCGGGCATACAAATAATGCATTTTCTTACACAACTTACGCCGGGGATAAACCGCGTAGCGCCCTTTCTTCTGTTTTTTGATTGAAGGGAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTTTAAGCGTAAcagc                                  >  2:728294/1‑150 (MQ=255)
                                                                                                                            cGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca                         >  2:639069/1‑150 (MQ=255)
                                                                                                                             gggCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca                         >  2:728218/1‑149 (MQ=255)
                                                                                                                                   aCAAATCATGCATTTCCTGACACAACTGACGCAGGGAATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg                 >  2:384425/1‑151 (MQ=255)
                                                                                                                                            gCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTgg         <  2:553449/150‑1 (MQ=255)
                                                                                                                                            gCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGc        >  1:708916/1‑151 (MQ=255)
                                                                                                                                                  ccTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAtc  >  2:586470/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                       
CGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATC  >  NC_000913/718523‑718819

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