Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,504,799	0	A	G	38.5%	52.3 / 12.1	39	G361G (GGA→GGG)	insQ	IS609 transposase B
Reads supporting (aligned to +/- strand):  ref base A (3/21);  new base G (15/0);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 3.25e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACA  >  NC_000913/1504662‑1504949                                                                                                                                             |                                                                                                                                                       tACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGGTCTT‑TCCTGTGGGGAGATGGTg                                                                                                                                               >  2:166519/1‑147 (MQ=255) tACGGTGTTCAAACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCa                                                                                                                                             <  2:365332/149‑1 (MQ=255)       ggTAATACGGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAgg                                                                                                                                       <  2:660870/149‑1 (MQ=255)            tcccggaaaaGGGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTc                                                                                                                                          <  2:109479/134‑1 (MQ=255)            tacgccGAAAAAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGACCGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGgccg                                                                                                                                    <  1:299095/142‑1 (MQ=255)              cAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑TCCTGTGGGGAGATGGTGGAGGCAGGGCGCCCgct                                                                                                                              >  2:674361/1‑149 (MQ=255)                  gAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGCATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGGGC‑ACCGTTCTT‑CGCCGTGGGGAGAAGGGGCGGCTAGGGCGccccgtgt                                                                                                                            >  1:362919/1‑144 (MQ=255)                    aaGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGGGAGAAGGTGCAGGCAGGCCGCCCGGTgga                                                                                                                           >  2:74183/1‑146 (MQ=255)                    aaGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCa                                                                                                                        <  2:26428/151‑1 (MQ=255)                      gagaACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑AGGAA‑CCCGGTCTT‑GCCCGGGGGGGGGACGTGCAGTCAGGccccccccttaagaagg                                                                                                                      >  1:156737/1‑140 (MQ=255)                       agaACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGGAC‑GCCGGTCTT‑GCCTGTGGGGAGATGGGGGAGGCAGGCCGCCCGTTGAAGCAgg                                                                                                                      >  2:329536/1‑150 (MQ=255)                       agaACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGGCAGGCCGCCCGTTGAAGCAGGa                                                                                                                     <  1:515829/151‑1 (MQ=255)                                    acaAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGGACGGCCGTCCTT‑‑CCTGGGGGGGGATGGGGGAGGCCGGGCGCCCGCGGAAGCCGGAAACCACCGaaaa                                                                                                         >  1:94194/1‑149 (MQ=255)                                    acaAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGGAC‑GCCGCTCTT‑TCCCTTGGGGGGATGGGGGAGGGAGGGCCCCCCCTGAAGCAAGAAccccccccaa                                                                                                          >  2:284881/1‑145 (MQ=255)                                     caAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGGAC‑GCCCCTCCT‑TCCCGTGGGGGGATGGGGCAATCACGGCGCCCGTTGGAGCAAGGACCCACCCaaa                                                                                                          >  2:260798/1‑148 (MQ=255)                                         gTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGGAC‑GCCGGTCTT‑GCCTGGGGGGGGATGGGGCAGGCAGGCCGCCCGTTGTAACAGGAACCCACCGAAATAAtt                                                                                                     >  2:58543/1‑149 (MQ=255)                                                cAGATGCCAGGTATGTGGATATCCAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCTTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACa                                                                                            <  1:537570/151‑1 (MQ=255)                                                     gCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGcat                                                                                        <  1:264837/150‑3 (MQ=255)                                                       cAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCCTTCTT‑TCCTGTGGAGGGATGGTGCAGGCAGGCCGCCCCGTGAAGCAGGAACCCACCGAAAAGATTTAGGCGACAGcc                                                                                         >  2:456180/1‑147 (MQ=255)                                                               tgGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCCTTTTT‑GTCTGTGGGGAGATGGGGGAGTCAGGGCCCCCGTTGAATCAGGGAACCACCGAAAAAATTCAGGGGACAGGCTGAACGTAg                                                                                >  2:484651/1‑148 (MQ=255)                                                                      aCAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCCTTCTT‑TCCTGTGGAGAGATGGGGCAGTCAGGCCGCCCGTTTAAGCAGGAACCCACCGAAAAGATTTAGGCGGCAGGCTGAACGGAGCAGGGAGcc                                                                       >  2:330158/1‑150 (MQ=255)                                                                        aGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACg                                                                    <  2:678201/151‑1 (MQ=255)                                                                           gAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGAGACCC‑GCC‑TTCTTCGCTTTAGGGGGGGTGG‑GCAGTCAACCCCCCCCCTGAAGCACACGACCACCCAAACGATTAAGGGGGCAGCCTGGACCTAGAAggggatcaccaac                                                                 >  1:660561/1‑141 (MQ=255)                                                                           gAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTcc                                                                 <  2:108603/151‑1 (MQ=255)                                                                            aaCGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCt                                                                <  2:25075/151‑1 (MQ=255)                                                                              cGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTc                                                              <  2:591353/151‑1 (MQ=255)                                                                               gCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCa                                                             <  1:546108/151‑1 (MQ=255)                                                                                  gATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑CCCGGTCTT‑GCCTGTGGGGAGATGGGGCAGGCAGGCCGCCCGGTGGAAGAGGAACCCACCCAAAAGATTTAAGCGGCAGCCCGAAAGGAGCAGGGGACCACGTCCCTCCggg                                                          >  2:669912/1‑151 (MQ=255)                                                                                    tGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGATAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCATGGc                                                         <  1:74128/150‑1 (MQ=255)                                                                                      tAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCCTTTTT‑TCCTGTGGGGAGATGGTGGAGTCAGGGCGCCCGCTGAAGCAGGAACCCCCCGAAAAAATTTAAGCGACAGCCTGAACGTAGCAGGGGTCCCCGTCCCTCAGAGCggg                                                      >  1:491781/1‑149 (MQ=255)                                                                                       aaaTGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTgg                                                     <  2:94229/151‑1 (MQ=255)                                                                                            gcgcTCGCAACATTTTAGCGGCGG‑GGGAC‑GCCGTTCTT‑TCCTGTGGAGAGATGGTGCAGGCAGGGCGCCCGTTGAAGCAGGAACCCACCCAAATGATTTAGGCGGCAGCCCGAACGGAGCAGGGGTCCACCTCCTTCAGGGCGgggggg                                                  >  2:485062/1‑149 (MQ=255)                                                                                              gcTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCCTTTTT‑TCCCGTGGGGAGATGGTGGAGGCAGGCCGCCCGGTGGAGGAGGAAACCACCGAAAAGATTTAAGGGGCAGCCCGGACGGAGCAGGGGGCCACCTCCTTCCggggggggggg                                                  >  2:301003/1‑139 (MQ=255)                                                                                                 cGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGACCCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCaa                                           <  1:120336/151‑1 (MQ=255)                                                                                                    aaCATTGCAGCGGCGG‑GGCAC‑GCCGTTGTT‑GCCTGTGGAGAGATGGTGCAGTCATGCCTCCCGTTGACGCAGGATCCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTTGAGGATGGCAata                                         <  1:637507/150‑1 (MQ=255)                                                                                                           tAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGtt                                   <  1:606945/149‑1 (MQ=255)                                                                                                                   g‑ggCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCg                         <  1:619034/151‑1 (MQ=255)                                                                                                                            ccGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAAt                  <  2:179689/151‑1 (MQ=255)                                                                                                                                         gtgGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAAcaccac        <  2:259101/149‑1 (MQ=255)                                                                                                                                             agagaTGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAAc   >  2:389668/1‑150 (MQ=255)                                                                                                                                             agagaTGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACa  <  1:342076/151‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                       TACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGG‑GGCAC‑GCCGTTCTT‑GCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACA  >  NC_000913/1504662‑1504949

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