Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913718,6900TG43.8% 74.8 / 14.3 48T591P (ACC→CCC) speFornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base T (5/22);  new base G (21/0);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 3.22e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGG  >  NC_000913/718545‑718838
                                                                                                                                                   |                                                                                                                                                     
tCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAAAAATGCATTTTCTGACAC‑AAATGACGCAGGGGATaa                                                                                                                                                   >  2:827640/1‑150 (MQ=255)
     gcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGGATACAAAACATGGATTTTCTGACAC‑AACTGACGCAGGGGATAAACGGc                                                                                                                                              >  1:823015/1‑150 (MQ=255)
             aGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGctc                                                                                                                                     <  2:301247/151‑1 (MQ=255)
             aGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGct                                                                                                                                      <  1:679719/150‑1 (MQ=255)
                       ggCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTg                                                                                                                            <  1:192020/150‑1 (MQ=255)
                       ggCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTg                                                                                                                            <  1:901649/150‑1 (MQ=255)
                       ggCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGt                                                                                                                           <  1:718548/151‑1 (MQ=255)
                        gCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGtt                                                                                                                          <  2:582659/151‑1 (MQ=255)
                         cttcttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCCTACAAAAAATGCATTTTCTTACAC‑AACAGACGCAAGGGATAAACGGCGTAGCGGTCTTCCTTCTGttt                                                                                                                         >  1:299259/1‑151 (MQ=255)
                              tGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAAAAATTCATTTTCTGACAC‑CACCCGCGCCGGGGATAAAAGGCGGAGCGGCCCTCATGCTGgtt                                                                                                                         >  2:78666/1‑146 (MQ=255)
                                 ggATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAAACATGCATTTTCTTACAC‑AACTGACGCAGGGGATAAACCGCGGAGCGCTCTTCATGCTGttttt                                                                                                                       >  1:883548/1‑143 (MQ=255)
                                    ttCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAAAAATGCTTTTTCTCACAC‑AAACTACGCAGGGGATAAACGGCGTTGCGCTCTTTTTTCTGttttt                                                                                                                       >  2:759362/1‑140 (MQ=255)
                                     tCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAAACATGCATTTTCTGACAC‑AACTGACGCAGGGGATAAACCGCGGAGCGCTCTTCATGCTGTTTGTAGATgg                                                                                                                 >  1:270623/1‑147 (MQ=255)
                                      cATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGGATAAAAAAAAAGCATTTTTTGACAC‑AACCGACGCCGGGGATAAACAGCGTTGGGCCCTTCTTGCTGttttttggtggga                                                                                                               >  2:817549/1‑138 (MQ=255)
                                        tGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCACACAAAAAAAGCCTTTCCTGCCCC‑CAACGCCGCAAGGGGTAACAAGCGTAGGGCCCTTCATTCTGttttttgtggg                                                                                                                 >  2:899034/1‑136 (MQ=255)
                                           tGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGGATACAAAAAAATCATTTTCTGACAC‑AAAATACGCAGGGGATAAACGGCGGAGCGCTCTTCATGCTGTTTTTAGATGggagggagga                                                                                                        >  2:549291/1‑148 (MQ=255)
                                                 gTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCACAAAAATAAATCATTTTCTGAAAC‑AAACCACCCCGGGGATAAACCGCGTGGCGCCACTTACGCCGGTTTTTGATGGAAGGGAGCaactc                                                                                                    >  2:752030/1‑144 (MQ=255)
                                                   gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCt                                                                                                <  2:252141/150‑1 (MQ=255)
                                                      aaGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATg                                                                                            <  2:80033/151‑1 (MQ=255)
                                                       aGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATgg                                                                                           <  2:409024/151‑1 (MQ=255)
                                                        gTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGgc                                                                                          <  2:207805/151‑1 (MQ=255)
                                                          gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACCAAACATTCCTTTTCTTACAC‑AACATACGCAGGGGGTAACCGGCGTAGCGCCCCTCATGCTGTTTGTTGATGGAAGGAAGGACCTCTGCTTATTgc                                                                                          >  2:630136/1‑149 (MQ=255)
                                                          gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGgcg                                                                                         <  1:387672/150‑1 (MQ=255)
                                                           cTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAAATAATTATTTTCTTACACAAAATGCCGC‑GGGGGTAAACCGCGGAGCGCGCTTCTTTCTGTTTTTTGAAGGGAGGGAGGAcat                                                                                                     >  2:325814/1‑135 (MQ=255)
                                                               ttACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAAACATGCATTTTCTGACAC‑AACTGACGCAGGGGATAAACGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGGAGCACTTCTGCTAATAGCGcggcgg                                                                                   >  1:75863/1‑151 (MQ=255)
                                                                  cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACt                                                                                <  1:845164/151‑1 (MQ=255)
                                                                   ggAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCCTACAAAAAATGCAATTTCTTCCAC‑AACTGACGCAGGGGATAAACGGCGGAGGGCTCTTCATGCTGTTTTTTGATGGGAGGGAGCACTTCTGCTAATTGCGCGGCGGagt                                                                                >  1:275573/1‑148 (MQ=255)
                                                                   ggAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTc                                                                               <  2:447052/151‑1 (MQ=255)
                                                                       cATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAAACATGCATTTTCTTACAC‑AACTGACGCCGGGGATAAACGGCGTAGGGCTCTTCATGGTGTTTGTTGATGGAAAGTAGCACTTCTTCTAATGGCGCGTCGGGCCCaa                                                                             >  1:688902/1‑149 (MQ=255)
                                                                         tctcTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCCTACAAAACATTCCTTTTCTGACAC‑AAACGACGCAGGGGATAAACCGCGGAGCGCTCTTCCTGCTGTTTTTTGATGGGAGGGAGCACTTCTGCTAATGGCGCGTGGGGCTCaa                                                                             >  1:643007/1‑147 (MQ=255)
                                                                           tctTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATTCATTTTCTGGCAC‑AACTGACGCAGGGGATAAACGGCGTAGCGCTCTTCTTGCTGTTTTTAGATGGAAGGGAGGACTTCTGCTAATGGCGCGGCGGGCTCAAGCA‑Gtt                                                                      >  2:481609/1‑151 (MQ=255)
                                                                              ttttGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑Gtt                                                                      <  2:151657/148‑1 (MQ=255)
                                                                               tttGCAGTTGTTTCACGTTGTGGCGGGGCATA‑AAAAAATAACTTTCCTTTCAC‑CAATCACGCAGGGGATAAACGGCGGAGCGGCCTTCTTGCTTTTTTGTGGTTGAAGGGAGCACTTCTCCTCAATGCGCGTTGGaattaaaacggtttt                                                                    >  1:284069/1‑136 (MQ=255)
                                                                               tttGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑Gttttt                                                                   <  2:491772/150‑1 (MQ=255)
                                                                               tttGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑Gttttt                                                                   <  2:486081/150‑1 (MQ=255)
                                                                               tttGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAAACATTCATTTTCTTACAC‑AACTGACGCAGGGGATAAACCGCGGAGCGCTCTTCATGCTGTTTGTTGATGGAAGGGAGGACTTCTGCTTATGGGGCGGCGGGCTCAAGCA‑Ggtttt                                                                   >  1:383347/1‑150 (MQ=255)
                                                                                ttGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTc                                                                  <  2:737271/150‑1 (MQ=255)
                                                                                       tgttTCACGTTGTGGC‑GGGCCTACAAAAAAAGCATTTTCTGGCAA‑AACAGACGCCGGGGATAAACCGCGGAGCGCCCTTCATTCTGTTTTTAGATGGAAGGGAGCAATTTTGCTAATTGCGCGGCGGGGTCAAACAGGTTTTTCGagc                                                              >  2:42966/1‑146 (MQ=255)
                                                                                               gTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTCGAAGCGTACCAGCAgg                                                  <  2:819249/151‑1 (MQ=255)
                                                                                                      c‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTCGAAGCGTACCAGCAGGGcaacaa                                           <  1:868410/151‑1 (MQ=255)
                                                                                                          gCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTCGAAGCGTACCAGCAGGGCAACAAGtt                                        >  1:643378/1‑151 (MQ=255)
                                                                                                                 aaTCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAg                                 <  1:935509/151‑1 (MQ=255)
                                                                                                                  aTCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGt                                >  1:121308/1‑151 (MQ=255)
                                                                                                                      tGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTgg                            >  1:841522/1‑151 (MQ=255)
                                                                                                                                 acac‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAtcttct                 >  2:894447/1‑151 (MQ=255)
                                                                                                                                           aCGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTc        <  1:420894/151‑1 (MQ=255)
                                                                                                                                            cGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTca       <  2:153339/151‑1 (MQ=255)
                                                                                                                                                  gTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAgg  >  2:820065/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGC‑GGGCATACAAATCATGCATTTCCTGACAC‑AACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCA‑GTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGG  >  NC_000913/718545‑718838

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