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

TCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGA  >  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)
                             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)
                                                                                                                             cTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAtct                     >  1:673654/1‑151 (MQ=255)
                                                                                                                                acacAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATAtcttct                  <  1:10553/151‑1 (MQ=255)
                                                                                                                                 cacaACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTg                 <  2:525207/151‑1 (MQ=255)
                                                                                                                                        gACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGt          >  2:640777/1‑151 (MQ=255)
                                                                                                                                              gggTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTcagcag    <  2:299846/151‑1 (MQ=255)
                                                                                                                                                 tataACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGa  >  2:438653/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
TCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATACAAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGAAGTTTGGCCATATCTTCTGCCGGAGTCAGCAGGA  >  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: