Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913718,6900TG52.2% 48.5 / 74.5 47T591P (ACC→CCC) speFornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base T (6/16);  new base G (24/0);  total (31/16)
Fisher's exact test for biased strand distribution p-value = 7.53e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-09

TCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGA  >  NC_000913/718545‑718839
                                                                                                                                                 |                                                                                                                                                      
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)
                             ttGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTAAAAAACATGCATTTTCTTACACAACTGACGCAGGGGATAAACCGCGGAGCGCCCTTCATGCT‑Gttttt                                                                                                                       >  1:646759/1‑147 (MQ=255)
                              tGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATGCATTTCCTTACACAACTGACGCCGGGGATAAACGGCGTAGCGCTCTTCATTCT‑GGTTGt                                                                                                                       >  1:423875/1‑148 (MQ=255)
                                 ggATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACAAATGACGCAGGGGATAACCGGCGGAGCGCTCTTCATGCT‑GTTTTTTGATgg                                                                                                                 >  1:592021/1‑151 (MQ=255)
                                  gATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGa                                                                                                                <  2:754535/151‑1 (MQ=255)
                                      cATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGt                                                                                                            <  1:877407/151‑1 (MQ=255)
                                      cATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGt                                                                                                            <  2:825094/151‑1 (MQ=255)
                                          cTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAATCATGCATTTCCTTCCACAACTGAAGCCGGGGGTAAACCGCCTAGCGCTCCTTCTGCT‑GTTTTTTGATGGAAGGGAg                                                                                                          >  1:551825/1‑149 (MQ=255)
                                              cgcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCCAAAAAACAAGCGTTTCCTTCCACAAATTACGGCGGGGGTAAACGGCCGAGGGCGCCTTATGCT‑CTTTTGTGATAGAAGGTAGAAca                                                                                                      >  2:813094/1‑148 (MQ=255)
                                               gcgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTACACAACTGACGCAGGGGATAAACGGCGGAGCGCTCTTCCTGCT‑GTTTGTTGATGGGAGGGAGCAcct                                                                                                     >  2:239345/1‑147 (MQ=255)
                                                       aGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTTCTGACACAACAGACGCCGGGGATAAACCGCGTAGCGCTCTTCATGCT‑GGTTGTAGATGGGAGGGAGCACTTTTGCTTAttg                                                                                           >  1:902549/1‑149 (MQ=255)
                                                        gTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGgc                                                                                          <  2:201326/151‑1 (MQ=255)
                                                          gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTTTCACAAATGACGCAGGGGATAAACGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGGAGCACTTCTGCTAATGgcg                                                                                         >  1:87625/1‑150 (MQ=255)
                                                                tACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTTCTGACACAACTGACGCAGGGGATAAACGGCGGAGCGCTCTTTCTTCT‑GTTTTTAGATGGAAAGGAGGACTTCTTCTAATGGCGCGTCgg                                                                                   >  1:862701/1‑150 (MQ=255)
                                                                tACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCgg                                                                                   <  1:653299/150‑1 (MQ=255)
                                                                 acggtgcATCTATTTTTGCATTTGAGTCACGTTGTGGAGAGCATACAAGTCATGCAGTTCCTGTCACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTc                                                                                     <  2:633982/141‑1 (MQ=255)
                                                                 aCGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTAAAAAATATTCATTTTCTTACACCACTGACGCAGGGGATAAACGGCGGAGCGCTCTTCTTGCT‑GTTTTTAGATGGGAGGGAGGACCTCTGCTTATGGGGcggcggg                                                                                  >  2:767964/1‑149 (MQ=255)
                                                                  cGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAAAAAACATTCATTTTCTTTCCCCACTTACGCAGGGGGTAAACGGCGTAGCGCCCTTCCTGCT‑GTTTTTAGATGGAAGGGAGCACTTTTGCTAATGGGGcggcgggc                                                                                 >  2:695120/1‑148 (MQ=255)
                                                                   ggAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAAACATTCATTTCCTTACACAACTGACGCAGGGGATAAACCGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGGAGGTAGGACTTCTGCTAATGGCGCGTCGGGCTc                                                                               >  1:201988/1‑151 (MQ=255)
                                                                       cATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATTATTCTTTTTCTGACACAAATGACGCAGGGGGTAAACGGCGGAGCGCTCTTCCTGCT‑GTTTGTAGATGGAAGGGAGCACTTCTTCTAATGGCGCGGCGGACTCaa                                                                             >  2:826543/1‑149 (MQ=255)
                                                                       cATCTCTTTTAGCAGTTGTTTCACGTTGTGGCGGGGCAACAAAACATGCATTTTCTTAAACAAATGAAGCAGGGGTTAAACCGCGTAGCGGCCTTCTTGAT‑GTTTTTTGAAGGGAAGGAGGACTTCTTCTTATGGCGGGGCGGCGTCaa                                                                             >  1:828752/1‑149 (MQ=255)
                                                                          ctctTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTACAAAACATGCTTTTTCTTACCCAACCCACGCCGGGGATAAACAGCGGAGCGCTCTTCCTGCT‑CTTTGTAGATGGAAGGGAGCACCTCTGCCCATGGGGCGTCgg                                                                                   >  2:464286/1‑140 (MQ=255)
                                                                            ctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGt                                                                        <  2:547818/149‑1 (MQ=255)
                                                                               tttGCAGTTGTTTCACGTTGTGGCGGGGATACAAAACATTGATTTTCTTACACCACTGACGCAGGGGATAACCCGCGGAGCGCTCTTTATTCT‑GTTTGTAGGTGGAAGGGAGGACTTCTGCTTATTGCGCGCCGGGCTTCAGCAGgttttt                                                                   >  2:760422/1‑150 (MQ=255)
                                                                                 tGCAGTTGTTTCACGTTGTGGCGGGGATAAAAAACATGCATTTTCTGAAACAACAGCCGGAAGGGAAAAACCGCGTCGGGCTCTTTATGCTGGGTTGTGGA‑GGAAGGGAGGACTTCTGCTAAAGGCGCGTCGGGCTCAAGCAGgtttttg                                                                  >  1:152227/1‑148 (MQ=255)
                                                                                     gttgttTCACGTTGTGGCGGGCCTACAAAACATGCATTTTCTTACACAACTGACGCAGGGGATAAACGGCGGAGCGCTCTTCATTCT‑GTTTTTAGATGGGAGGGAGCACTTCTTCTAATGGCGCGTCGGGCTCAAGCAGGTTTTTGAAGCg                                                             >  2:13036/1‑151 (MQ=255)
                                                                                     gttgttTCACGTTGTGGCGGGCATACAAAACATGCATTTTCTGAAAAAAATAACGGCGGGGGTAAACGGCGGAGCACTCTTCTTGCT‑GTTTTTAGATGGGAGGGAGCACCTTTTCTTATGGGGCGTCGGACTCACACAGTTTTTTGAAGc                                                              >  1:228654/1‑150 (MQ=255)
                                                                                        gttTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTa                                                           <  2:288756/150‑1 (MQ=255)
                                                                                          ttCACGTTGTGGCGGGCATACAAAAAAATCATTTTCTGACACAAATGACGCAGGGGATAAACGGCGCAGCGCCCTTTATGCT‑GTTTTTAGATGGAAGGGAGCACTTCTGCTAATGGCGCGTCGGGCTCAAGCAGGTTTTTGAAGGGTAcc                                                         >  1:605497/1‑150 (MQ=255)
                                                                                             aCGTTGTGGCGGGCATACAAATAATGCATTTTCTTACACAACTTACGCCGGGGATAAACCGCGTAGCGCCCTTTCTTCT‑GTTTTTTGATTGAAGGGAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTTTAAGCGTAAcagc                                                      >  2:728294/1‑150 (MQ=255)
                                                                                                      cGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca                                             >  2:639069/1‑150 (MQ=255)
                                                                                                       gggCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaaca                                             >  2:728218/1‑149 (MQ=255)
                                                                                                             aCAAATCATGCATTTCCTGACACAACTGACGCAGGGAATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg                                     >  2:384425/1‑151 (MQ=255)
                                                                                                                      gCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTgg                             <  2:553449/150‑1 (MQ=255)
                                                                                                                      gCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGc                            >  1:708916/1‑151 (MQ=255)
                                                                                                                            ccTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAtc                      >  2:586470/1‑151 (MQ=255)
                                                                                                                             cTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAtct                     >  1:673654/1‑151 (MQ=255)
                                                                                                                                acacAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAtcttct                  <  1:10553/151‑1 (MQ=255)
                                                                                                                                 cacaACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTg                 <  2:525207/151‑1 (MQ=255)
                                                                                                                                        gACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGt          >  2:640777/1‑151 (MQ=255)
                                                                                                                                              gggTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTcagcag    <  2:299846/151‑1 (MQ=255)
                                                                                                                                                 tataACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGa  >  2:438653/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
TCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCT‑GTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGA  >  NC_000913/718545‑718839

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