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

TCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGG  >  NC_000913/718545‑718838
                                                                                                                                                   |                                                                                                                                                     
tCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCATTTTCTGACACAAATGACGCAGGGGATaa                                                                                                                                                   >  2:827640/1‑150 (MQ=255)
     gcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATGGATTTTCTGACACAACTGACGCAGGGGATAAACGGc                                                                                                                                              >  1:823015/1‑150 (MQ=255)
             aGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGctc                                                                                                                                     <  2:301247/151‑1 (MQ=255)
             aGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGct                                                                                                                                      <  1:679719/150‑1 (MQ=255)
                       ggCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTg                                                                                                                            <  1:192020/150‑1 (MQ=255)
                       ggCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTg                                                                                                                            <  1:901649/150‑1 (MQ=255)
                       ggCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGt                                                                                                                           <  1:718548/151‑1 (MQ=255)
                        gCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGtt                                                                                                                          <  2:582659/151‑1 (MQ=255)
                         cttcttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAAAATGCATTTTCTTACACAACAGACGCAAGGGATAAACGGCGTAGCGGTCTTCCTTCTGttt                                                                                                                         >  1:299259/1‑151 (MQ=255)
                              tGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATTCATTTTCTGACACCACCCGCGCCGGGGATAAAAGGCGGAGCGGCCCTCATGCTGgtt                                                                                                                         >  2:78666/1‑146 (MQ=255)
                                 ggATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACAACTGACGCAGGGGATAAACCGCGGAGCGCTCTTCATGCTGttttt                                                                                                                       >  1:883548/1‑143 (MQ=255)
                                    ttCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAAAATGCTTTTTCTCACACAAACTACGCAGGGGATAAACGGCGTTGCGCTCTTTTTTCTGttttt                                                                                                                       >  2:759362/1‑140 (MQ=255)
                                     tCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTGACACAACTGACGCAGGGGATAAACCGCGGAGCGCTCTTCATGCTGTTTGTAGATgg                                                                                                                 >  1:270623/1‑147 (MQ=255)
                                      cATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAAAAAAAAAGCATTTTTTGACACAACCGACGCCGGGGATAAACAGCGTTGGGCCCTTCTTGCTGttttttggtggga                                                                                                               >  2:817549/1‑138 (MQ=255)
                                        tGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCACACAAAAAAAGCCTTTCCTGCCCCCAACGCCGCAAGGGGTAACAAGCGTAGGGCCCTTCATTCTGttttttgtggg                                                                                                                 >  2:899034/1‑136 (MQ=255)
                                           tGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAAAAATCATTTTCTGACACAAAATACGCAGGGGATAAACGGCGGAGCGCTCTTCATGCTGTTTTTAGATGggagggagga                                                                                                        >  2:549291/1‑148 (MQ=255)
                                                 gTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCACAAAAATAAATCATTTTCTGAAACAAACCACCCCGGGGATAAACCGCGTGGCGCCACTTACGCCGGTTTTTGATGGAAGGGAGCaactc                                                                                                    >  2:752030/1‑144 (MQ=255)
                                                   gggAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCt                                                                                                <  2:252141/150‑1 (MQ=255)
                                                      aaGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATg                                                                                            <  2:80033/151‑1 (MQ=255)
                                                       aGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATgg                                                                                           <  2:409024/151‑1 (MQ=255)
                                                        gTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGgc                                                                                          <  2:207805/151‑1 (MQ=255)
                                                          gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACCAAACATTCCTTTTCTTACACAACATACGCAGGGGGTAACCGGCGTAGCGCCCCTCATGCTGTTTGTTGATGGAAGGAAGGACCTCTGCTTATTgc                                                                                          >  2:630136/1‑149 (MQ=255)
                                                          gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGgcg                                                                                         <  1:387672/150‑1 (MQ=255)
                                                           cTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAATAATTATTTTCTTACACAAAATGCCGCGGGGGTAAACCGCGGAGCGCGCTTCTTTCTGTTTTTTGAAGGGAGGGAGGAcat                                                                                                     >  2:325814/1‑135 (MQ=255)
                                                               ttACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTGACACAACTGACGCAGGGGATAAACGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGGAGCACTTCTGCTAATAGCGcggcgg                                                                                   >  1:75863/1‑151 (MQ=255)
                                                                  cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACt                                                                                <  1:845164/151‑1 (MQ=255)
                                                                   ggAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAAAATGCAATTTCTTCCACAACTGACGCAGGGGATAAACGGCGGAGGGCTCTTCATGCTGTTTTTTGATGGGAGGGAGCACTTCTGCTAATTGCGCGGCGGagt                                                                                >  1:275573/1‑148 (MQ=255)
                                                                   ggAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTc                                                                               <  2:447052/151‑1 (MQ=255)
                                                                       cATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACAACTGACGCCGGGGATAAACGGCGTAGGGCTCTTCATGGTGTTTGTTGATGGAAAGTAGCACTTCTTCTAATGGCGCGTCGGGCCCaa                                                                             >  1:688902/1‑149 (MQ=255)
                                                                         tctcTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATTCCTTTTCTGACACAAACGACGCAGGGGATAAACCGCGGAGCGCTCTTCCTGCTGTTTTTTGATGGGAGGGAGCACTTCTGCTAATGGCGCGTGGGGCTCaa                                                                             >  1:643007/1‑147 (MQ=255)
                                                                           tctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATTCATTTTCTGGCACAACTGACGCAGGGGATAAACGGCGTAGCGCTCTTCTTGCTGTTTTTAGATGGAAGGGAGGACTTCTGCTAATGGCGCGGCGGGCTCAAGCAGtt                                                                      >  2:481609/1‑151 (MQ=255)
                                                                              ttttGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGtt                                                                      <  2:151657/148‑1 (MQ=255)
                                                                               tttGCAGTTGTTTCACGTTGTGGCGGGGCATAAAAAAATAACTTTCCTTTCACCAATCACGCAGGGGATAAACGGCGGAGCGGCCTTCTTGCTTTTTTGTGGTTGAAGGGAGCACTTCTCCTCAATGCGCGTTGGaattaaaacggtttt                                                                    >  1:284069/1‑136 (MQ=255)
                                                                               tttGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGttttt                                                                   <  2:491772/150‑1 (MQ=255)
                                                                               tttGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGttttt                                                                   <  2:486081/150‑1 (MQ=255)
                                                                               tttGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATTCATTTTCTTACACAACTGACGCAGGGGATAAACCGCGGAGCGCTCTTCATGCTGTTTGTTGATGGAAGGGAGGACTTCTGCTTATGGGGCGGCGGGCTCAAGCAGgtttt                                                                   >  1:383347/1‑150 (MQ=255)
                                                                                ttGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTc                                                                  <  2:737271/150‑1 (MQ=255)
                                                                                       tgttTCACGTTGTGGCGGGCCTACAAAAAAAGCATTTTCTGGCAAAACAGACGCCGGGGATAAACCGCGGAGCGCCCTTCATTCTGTTTTTAGATGGAAGGGAGCAATTTTGCTAATTGCGCGGCGGGGTCAAACAGGTTTTTCGagc                                                              >  2:42966/1‑146 (MQ=255)
                                                                                               gTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAgg                                                  <  2:819249/151‑1 (MQ=255)
                                                                                                      cGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaacaa                                           <  1:868410/151‑1 (MQ=255)
                                                                                                          gCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGtt                                        >  1:643378/1‑151 (MQ=255)
                                                                                                                 aaTCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAg                                 <  1:935509/151‑1 (MQ=255)
                                                                                                                  aTCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGt                                >  1:121308/1‑151 (MQ=255)
                                                                                                                      tGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTgg                            >  1:841522/1‑151 (MQ=255)
                                                                                                                                 acacAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAtcttct                 >  2:894447/1‑151 (MQ=255)
                                                                                                                                           aCGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTc        <  1:420894/151‑1 (MQ=255)
                                                                                                                                            cGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTca       <  2:153339/151‑1 (MQ=255)
                                                                                                                                                  gTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAgg  >  2:820065/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGG  >  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: