Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,622,250	0	A	G	41.9%	54.1 / 12.5	45	E161G (GAG→GGG)	yjjJ	putative protein kinase
Reads supporting (aligned to +/- strand):  ref base A (6/19);  new base G (18/0);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 2.28e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACC  >  NC_000913/4622102‑4622398                                                                                                                                                      |                                                                                                                                                         ccGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCggggagg                                                                                                                                                        <  2:389520/151‑1 (MQ=255)      gggCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGaa                                                                                                                                                    <  2:354054/150‑1 (MQ=255)          tttttAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTTGTTGGTCggggggggggaat                                                                                                                                                  >  2:568172/1‑144 (MQ=255)             ttAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGGTGGGCGGGGGGGGGAATTAt                                                                                                                                                >  2:497056/1‑147 (MQ=255)               aGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCg                                                                                                                                          <  2:27881/150‑1 (MQ=255)                         tGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTTGTTGGTCGGGGGGGGGGATTATCA‑GCGATGG‑Gtttc                                                                                                                                >  1:50276/1‑147 (MQ=255)                         tGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGGTGGGCGGGGGGGGGGATTATCA‑GCGATGG‑GTTAc                                                                                                                                >  1:18271/1‑149 (MQ=255)                               aGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGGATTATCA‑GCGATGG‑GTTACTGCacaa                                                                                                                         >  1:682804/1‑150 (MQ=255)                               aGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCacaaa                                                                                                                        <  1:138993/151‑2 (MQ=255)                                     ttAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGGCGGGGGGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAc‑cc                                                                                                                   >  1:678635/1‑149 (MQ=255)                                       aGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGc                                                                                                                <  1:347084/150‑1 (MQ=255)                                       aGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TAAGTAGGTGTGGGGAGGGAAATAGCA‑‑CTATGGCATTACTGAacaaatttc                                                                                                                    >  2:428892/1‑142 (MQ=255)                                              cAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTc                                                                                                        <  1:568447/151‑1 (MQ=255)                                                aCTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCct                                                                                                      <  2:684689/151‑1 (MQ=255)                                                       ctgactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTTGGTTGTGGGGGGGGGGGATTATCA‑GGGATGG‑AATAATTCAACAAAC‑CCTGCGGGAATTTcttttggt                                                                                                >  2:337198/1‑144 (MQ=255)                                                         gactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCa                                                                                             <  1:30757/150‑1 (MQ=255)                                                          actgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCaaa                                                                                           <  2:440888/151‑1 (MQ=255)                                                               tgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑CTGGTTGGGCGGGGGGGGGAATTATCA‑GCGATGG‑GGTACTGCACAAAAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTc                                                                                      >  1:376499/1‑151 (MQ=255)                                                                gatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTc                                                                                      <  2:676640/150‑1 (MQ=255)                                                                    taCGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGGCGGGGGGGGGGATTATCA‑GCGATGG‑GTTTCTTCACAAAAC‑CCTGCGGGAATTCCTCT‑GGATCAAAAACTCAcccc                                                                                 >  2:297430/1‑150 (MQ=255)                                                                       gTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTa                                                                              <  1:145947/151‑1 (MQ=255)                                                                          tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGGTGGGCGGGGGGGGGGATTATCA‑GCGATGG‑GTTACTGCACAACAC‑CCTGCGGGAATTCCTCT‑GGATCAAAAACTCACCCAtt                                                                               >  1:274650/1‑147 (MQ=255)                                                                          tctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑GTGGGTGGGCGGGGGGGGGAATTATCA‑GCGATGG‑GTTACTGCACAACAC‑CCTGCGGAAATTTCTCT‑GGATCAAAAACTCACCCATTACg                                                                            >  1:218819/1‑150 (MQ=255)                                                                             ggCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGa                                                                           <  2:210932/148‑1 (MQ=255)                                                                                  gaagaaGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGGGGGGAATTATTA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGGAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTgg                                                                   >  2:550152/1‑151 (MQ=255)                                                                                      aagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTTGTCCGGGGGGGGGATTATTA‑TCGATGG‑AATACTGTACAACAC‑ACTGCCGGAATTCCTGC‑GGATCAAAAAATCACCCATTACAAACAGCTGGCa                                                                 >  2:194311/1‑149 (MQ=255)                                                                                         aTGTGCTCTACGCCCTGACCGTATTTAACGTTGAAGATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTg                                                             <  2:46120/150‑1 (MQ=255)                                                                                             gCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGGTGGTCGGGGGGGGGAAATATTA‑GCGGTGG‑GTTACTGCACAAAAC‑CCTGCGGGAATTCCTCT‑GGATCAAAAAATCAACCATTACGAACAGCTGGCAAGTGa                                                            >  1:242983/1‑147 (MQ=255)                                                                                              ctctACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTTGGTGGTCGGGGCGGGGAATTATCA‑AAGGTGG‑ATTACTGCCCAACAA‑CCCGCGGGAATTCCTCT‑GGATCAAAAAACCACCCATTACGAACAGCCgg                                                                   >  2:671241/1‑139 (MQ=255)                                                                                              ctctACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGGATTATTA‑GCGATGG‑ATTAATTCACAACAA‑CCTGCGGAAAATCCTCT‑TGATAAAAAAATCAACCCTTAAGAACAACTGGGaa                                                                >  2:339860/1‑142 (MQ=255)                                                                                                   cGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGGGGGGAAATATCA‑GGGCTGG‑ATTACTGCACAAAAA‑CCCGCGGAAATTCCCCT‑GGATCAAAAAATCACCCATTACGAAAAGCTGGGAAGGGGTGCAc                                                       >  2:2841/1‑146 (MQ=255)                                                                                                   cGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑GTGGGTGGTCGGGGGGGGGTATTATCAGGCGA‑GG‑ATTTTTGCGCAAAAC‑CCTTCCGGAAATCCTCT‑GGATCAAAAACCCAACCACTACCAACAACTCGCAAGTGAAGCCCTgg                                                    >  1:376957/1‑149 (MQ=255)                                                                                                      ccTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGGTGGTCGGGGGGGGGGATTATCA‑GGGATGG‑ATTACTGCACAACAA‑CCTGCGGAAAATTCTCT‑GGGTCCAAAAATTACCCATTACGAACAACCGGCAAATGATGCACTGGCCgg                                                >  2:370646/1‑150 (MQ=255)                                                                                                            cGTATTTAACGGTGAATATACTGGCGGCTTGGTGGGGGGGGATGGGGAATT‑TCA‑ATGATGG‑ATTATGGTGAAATACACACGCGG‑AATTTCTCTAGCATCAAAAA‑TGACCCATTACGAACAAATGGCAAGCGCTGAACTGGGAGGAGaa                                            >  2:73495/1‑148 (MQ=255)                                                                                                             gTATTTAACGGTGAATATACTGGCGG‑GTTGGTGGTCGGGGAGGGGAAATATCA‑GCGATAG‑ATTACTGCACAACAC‑CCTGCGGGAAATCCTCT‑GGAACAAAAAATCACCCATTACGCAAAACTGGCAAGTGATGCACTGGCCGGAGgaa                                           >  2:648619/1‑148 (MQ=255)                                                                                                               aTTTAACGGTGAATATACTGGCGG‑TTGGGTGGGTGGGGGGGGGGATTATTA‑GCGGTGG‑ATTTCTGCACAACAA‑CCCGCGGGAATTTCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGGGATGCACTGGCAGGAGAAAAtgt                                       >  2:658855/1‑150 (MQ=255)                                                                                                                 ttAACGGTGAATATACTGGCGG‑TTGGGTGGGCGGGGAGGGGGATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCCGCGGGAATTTCTCT‑GGGTCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACCGGCAGGGGAAAATGTggg                                    >  2:319188/1‑151 (MQ=255)                                                                                                                   aaCGGCGAATATACTGGCGG‑CTTGTTGGTTGGGGAGGGGGATTATTA‑GTGATGG‑ATTTATGGACAACAC‑CCCGTGGGAATTCCTCT‑GGGTCCAAAACTCACCCATTACGAACAGCCGGCAAATGATGCACTGGCAGGAGAAATTGAggg                                    >  2:376795/1‑149 (MQ=255)                                                                                                                     cGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGGGGGGAATTATTA‑GCGATGG‑ATTTCTGCACAACAC‑CCCGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGGGATGCACTGGCAGGAGAAATTGTGGGgtct                                >  1:100087/1‑151 (MQ=255)                                                                                                                         gAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTg                            <  1:335357/151‑1 (MQ=255)                                                                                                                          aaTATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTtcttct                             >  1:209696/1‑149 (MQ=255)                                                                                                                          aaTATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgc                           <  2:109702/151‑1 (MQ=255)                                                                                                                               aCTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgcgggc                       <  2:397629/150‑1 (MQ=255)                                                                                                                                  ggCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCGCAGCGC‑TCTGCGGAAATTCCTCT‑GGATCAAAAGCTGACCGATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgcgggcgg                     <  2:162938/149‑1 (MQ=255)                                                                                                                                      g‑ttggttgGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGTCGagca               <  2:558894/151‑1 (MQ=255)                                                                                                                                              ggtCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCaaaa        <  1:550423/151‑1 (MQ=255)                                                                                                                                                 cGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAAttt     <  1:352784/151‑1 (MQ=255)                                                                                                                                                    ggaggggaATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTAcc  <  1:667937/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                         CCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGG‑TTGGTTGGTCGGGGAGGGGAATTATCA‑GCGATGG‑ATTACTGCACAACAC‑CCTGCGGAAATTCCTCT‑GGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACC  >  NC_000913/4622102‑4622398

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