Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	15,436	0	A	T	44.7%	39.4 / 19.9	38	intergenic (+138/‑9)	dnaJ/insL1	chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
Reads supporting (aligned to +/- strand):  ref base A (4/17);  new base T (17/0);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 2.20e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCT  >  NC_000913/15294‑15536                                                                                                                                               |                                                                                                      gCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                                                               <  1:947674/151‑1 (MQ=255)    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                                                                                            <  1:500000/151‑1 (MQ=255)             aaaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTActct                                                                                  <  1:440070/151‑1 (MQ=255)              aaaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACtctc                                                                                 <  1:216064/151‑1 (MQ=255)                aGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCAATCCCATAAACGCTAAATTAAAGGGTGTGGGGTTATGCCTGATATTATTTTAAGTGCCGATG‑AATTACtctccc                                                                               >  2:41106/1‑149 (MQ=255)                     gtccGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATaa                                                                          <  2:672490/149‑1 (MQ=255)                         gTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGCCCACTTAAAGGGTGTGGGATTTCGCCCGGTATGATTTTAAGTGCCGGTG‑AATTAATCTCCCGaaaa                                                                          >  2:1041973/1‑147 (MQ=255)                              cAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTACTCCCCCAAACGCCTACTTTAAGGGTGGGGGATTTCGCCTGATATGATTTTAAGTGCCGGAG‑AAATAATCTCACGaaaaataggt                                                                    >  1:1047011/1‑142 (MQ=255)                               aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCa                                                                <  2:651967/151‑1 (MQ=255)                                  ccTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTTACACCCCCAAATGCTAAATTTAGGGGTGGGGGATATCGCCTGGTTTGTTTTTAATTGCCGATG‑AATTATTTTCACCATAACTTGTGAGCAAtt                                                             >  1:817952/1‑151 (MQ=255)                                  ccTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCATAAACGCTAACTTAAGGGGTGTGGGATTTAGCCTGATATGATTTTAAGTGCCGATG‑AATTACTTTCCCGATAACTGGGCAGCAAtt                                                             >  2:456299/1‑151 (MQ=255)                                    tGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAACGCCAACTTAAAGGGTGTGGGGTTACGCCTGATATTATTTTAAGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAAATCt                                                           >  2:766160/1‑151 (MQ=21)                                     gggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTg                                                          <  1:251979/151‑1 (MQ=255)                                      ggTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCCCCTAAACGCCAACTTTAGGGGTGTGGGTTTTTGCCTGATATGATTTTACGTGCCGGTG‑AAATACTCTCACGATAAATGGTGAGCAAATCt                                                           >  2:58758/1‑149 (MQ=255)                                       gTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                                        <  2:946638/151‑1 (MQ=255)                                         aaaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGAGAACACCCCCCTAAACGCCAAATTAAAGGGTGGGGGATTATGGCTGGTATGATTTTACGTGGCGGTGAAATTATTTTC‑CCATAAACTGTCAGCAAATCTTGccc                                                      >  1:421515/1‑151 (MQ=255)                                          aaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCa                                                     <  1:458862/151‑1 (MQ=255)                                           aaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCat                                                    <  2:1017719/151‑1 (MQ=255)                                             aTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGCTCACTTTAGGGGTGTGGTATTACGCCTGGTATGATTTTACGTGGCGATG‑AATTAATCTCCCGATAAATGGTCAGCAAATCTTGccc                                                      >  2:179610/1‑147 (MQ=255)                                                   tGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCATAAACGCTAACTTAAGGGGTGTGGGATTACGCCTGATATGATTTTAAGTGCCGATG‑AATTAATCTCACGATAAATGGTCAGCAATTCTTGCCCATATTgg                                               >  2:351909/1‑148 (MQ=255)                                                      gTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACACTCCCCCAAACGCCAACTTTAGGGGTGTGGGATTTAGCCTGATATGATTTTAACTGCCCATG‑AATTAATCTCCCGCTAAATGGGCAGGAATTCTTGCCCCTATTTGTaaagc                                         >  2:601377/1‑148 (MQ=255)                                                         gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCga                                      <  1:606155/151‑1 (MQ=255)                                                          aaGATATGCGAGCACCTGTAAAGTGGCGGGGGACACCCCCCTAAACGCCAACTTAAGGGGTGTGGGATTACGCCCGGTATGATTTTAAGTGGCGCTG‑AATTATTCTCCCCAGAAATGGTCAGCAAATTTTGCCCACAtttggaaaacc                                        >  1:388862/1‑139 (MQ=255)                                                          aaGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaa                                     <  1:764985/151‑1 (MQ=255)                                                           aGATATGCGAGCACCTGTAAAGTGGCGGGGATCACCCCCATAAACGCACACAAACAGGGTGTGGTGATACCCCTCATATTACTTTAGGTGCCGAAT‑AATTAACCCCCCTAGAAATGGTCAGCAATCTTGGccccaatttgggtgccc                                        >  2:693566/1‑133 (MQ=255)                                                             atatGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaagaa                                  <  1:94378/151‑1 (MQ=255)                                                              tatGCGAGCACCTGTAAAGTGGCGGGGGACACCCCCCTAAACGCTAACCTAAAGGGTGTGGGATTATGCCTGATATGATTTTATGTGCCGGTG‑AATTACTCTCACGATAAATGGTCAGCAATTCTCGGCCCTAtt                                                 >  1:282805/1‑135 (MQ=255)                                                                 gCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGCTAACTTAAAGGGTGTGGGATTAAGCCTGGTATGATTTTAAGGGCCGATG‑AAATAATCTCACGATAAATGGTCAGCAAATCTGGCCCCTATTTGTAAACCCGAAGAACAgg                              >  2:238438/1‑151 (MQ=255)                                                                    aGCACCTGTAAAGTGGCGGGGAACACTCCCATAAACGCTCAAATAAAGGGTGTGGGATTTCGCCCGGTATTATTTTACGTGCCGATG‑GAATACTCCCCCCCTAAACGGGCAGCAAATTTTGCCCCTATTTGTAAACCCgaagaaaagg                              >  1:549324/1‑144 (MQ=255)                                                                     gCACCTGTAAAGTGGCGGGGGTCACTCCCCTCAACCCCAACTTTAAGGGTGTGGGATTTCGCCTGATATGATTTTAAGTGGCGATG‑AAATATTCTCACGATAACTTGGTAGCCAATCTTGCCCCTAATTGTAAAGCCgaagaaatcg                              >  1:631851/1‑143 (MQ=255)                                                                           gTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGc                    <  1:999765/151‑1 (MQ=40)                                                                              aaGTGGCGGGGATCACTCCCCTAAACGCTAACTTTAGGGGTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAAATGGTCAGCAATTTTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGccc                  >  1:803827/1‑148 (MQ=14)                                                                              aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACg                 <  1:32221/151‑1 (MQ=38)                                                                                gTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa               <  2:593459/151‑1 (MQ=37)                                                                                   gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTaa              <  1:544523/149‑1 (MQ=34)                                                                                    cGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGc           >  2:246302/1‑151 (MQ=33)                                                                                           cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGct    >  1:252918/1‑151 (MQ=31)                                                                                             cTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGctct  >  1:255420/1‑151 (MQ=11)                                                                                                                                               |                                                                                                      GCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATG‑AATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCT  >  NC_000913/15294‑15536

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