Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,730	0	C	A	43.1%	95.4 / 24.4	52	N310K (AAC→AAA)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base C (6/23);  new base A (22/0);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 2.46e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGC  >  NC_000913/3585‑3878                                                                                                                                                   |                                                                                                                                                      ccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCGCG‑CGGACTGGGAAAATaaa                                                                                                                                                   >  1:426670/1‑150 (MQ=255)      ccGACCTTGTTCGCTCTGTGTGACAAGCCGGAGACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCG‑GGTAATGGGAAAATaaaaaaaa                                                                                                                                              >  1:83363/1‑145 (MQ=255)       cGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGGGCG‑CGGACTGGGAAAATaaaaaaaaa                                                                                                                                             >  1:1104462/1‑144 (MQ=255)       cGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAAct                                                                                                                                            <  1:576798/151‑1 (MQ=255)                tCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATc                                                                                                                                    <  2:1129351/150‑1 (MQ=255)                 cGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCt                                                                                                                                   <  1:47507/150‑1 (MQ=255)                      tgtgtgACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCACG‑AGGACTGGGAAAAAAAAAGAAAATCTACAAACTGaaa                                                                                                                               >  2:165709/1‑149 (MQ=255)                          tgACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑GGGACCGGAGAAATAAATGAAAATCTACAATCTGAAAAAtaa                                                                                                                          >  1:255858/1‑148 (MQ=255)                            aCAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACaa                                                                                                                       <  2:981440/151‑1 (MQ=255)                                   ggAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACC‑AGGACCGGGAAAATAAATGAAAATCTTCAATCTGaaaaaaaaaaa                                                                                                                       >  2:544694/1‑136 (MQ=255)                                          gCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCt                                                                                                        <  2:115230/151‑1 (MQ=255)                                          gCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGGACTGGGAAAATAAATGAAAATCTACAATCTGAAAGATCACAAACAGCAGGGCA‑GCt                                                                                                        >  2:1083663/1‑151 (MQ=255)                                          gCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGGACTGGGAAAAAAAAAAAAAATCTAAAAACTTAAAAAACACAAAAAGCAGGGCA‑g                                                                                                          >  1:587381/1‑149 (MQ=255)                                              aGCGCGTTGCCGACTTGTTGGGCAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCCCC‑AGGACTCGAGAAATAAATAAAAAACCACAAACTTTCAAAAAACAACACGCAGGTCA‑CCttt                                                                                                      >  1:132907/1‑149 (MQ=255)                                                     tGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCAAG‑AGGACTGG‑GGACTAAAAAAAAAACTACAAAATTAAAAATAACAACGAACAGGGCA‑GATTTTCGCaa                                                                                                >  1:916834/1‑147 (MQ=255)                                                      gCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTAATGGGAAAATAAATGAAAATCTTCAATCTTAAAGAACCCAAAGAGCAGGGCC‑GCCTTTTGCAAGcc                                                                                             >  2:573991/1‑150 (MQ=255)                                                       ccGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCCCG‑AGGACTTGGAAAAAAAAAAAAAATCTACAAACTGAAAAAACACAACCAGAAGGGCA‑GCTTTtcgccaacc                                                                                             >  1:263509/1‑144 (MQ=255)                                                         gACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCG‑AGGACTGGAAAAATAAATGAAACTCTACAAACTGAAAAAACACAACGAGAAGGGCA‑GCTTTTCGCAAGCCgca                                                                                          >  2:585407/1‑148 (MQ=255)                                                           cTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑CGGACTGGGAAAATAAAAGAAAAACTACAAACTGAAAAAACACAACAAGCAGGGCA‑GCTTTTCGCCAGCCCTAAcc                                                                                       >  2:254872/1‑151 (MQ=255)                                                             ggTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAAcc                                                                                       <  2:700420/149‑1 (MQ=255)                                                             ggTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCa                                                                                     <  2:179835/151‑1 (MQ=255)                                                              gTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCAGG‑AGTACTGGAAAACTAAATGAAGCCCTACAATCTGAAAGATCACAACTAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAg                                                                                    <  2:1109216/151‑1 (MQ=255)                                                                   gTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGtt                                                                               <  2:1037071/151‑1 (MQ=255)                                                                    tAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGGAATGGGAAAATAAATGAAAAACTACAAACTGAAAAATCACAACGAGCAGGGCA‑GCTTTTCGCAAACCGTAAACCCggggg                                                                                >  2:859617/1‑148 (MQ=255)                                                                        aaCTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGGAATGGGAAAATAAATGAAACTCTACAAACTGAAAGATCACAACGAGCAGGGCA‑GCTTTTCGCAAGCCGTAAACCCGGGGGTggg                                                                            >  2:468404/1‑149 (MQ=255)                                                                          cTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGGAATGGGAAAAAAAATGAAAATCTACAAACTGAAAAATCACAACGAGCAGGGCA‑GCTTTGCGCAAGCCGGAAcccccgggggggg                                                                            >  1:855710/1‑137 (MQ=255)                                                                           tACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCC‑AGGAATGGGAAAAAAAAAAAAAATATAAAAACTTAAAGATAACAAAAAACAGGGCA‑GATTTTCGCAAGCCCCAAACCCGGGGGTgggggaa                                                                        >  1:580875/1‑146 (MQ=255)                                                                            aCCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGTTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCaaa                                                                        <  1:163590/149‑1 (MQ=255)                                                                                     aTCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCACG‑AGGAATGGGAAAAAAAATAAAACTATACAAACTGAAAAAACAAAAACAACAGGGCAGGCTTTG‑GCAAGCCGGAACCCCGGGGGTGGGGGAAAAAAAggggg                                                              >  1:486691/1‑149 (MQ=255)                                                                                          gAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGGAATGGGAAAATAAATTAAAATCTACAAACTTAAAGAACACAAACAACAGGGCA‑GCTTTTTGCAAACCGTAACCCAGGGGGTGGGGAAAAAAAAGGGGGTGgttt                                                        >  1:724362/1‑151 (MQ=255)                                                                                               ttttGTTCATATTTGCCGGCTGGATACGGCGGGCGCGCGAAGTCTTGGAAAA‑TTAATGTAACTGCAAAAGGTGAACGATAAAAAACAGACGGTCA‑GCGCTTCGGCAGCCAGGACCCCGGGGGGGGGGAAAAATAAAGGGGCGttttttt                                                     >  1:639839/1‑148 (MQ=255)                                                                                               ttttGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAAATAAATTAAAATCTACAATCTTAAAGATCACAACGAGCCGGTCA‑GCTTTGCGCAAGCCGTAAACCAGGGGGTGGGGGAAAAACAGGGGGTGgtttttt                                                     >  2:814131/1‑148 (MQ=255)                                                                                                    ttCATATTTGCCGGCTGGATACGGCGGGGGCACG‑AGGAATTGGAAAAAAAAAAAAAATCTAAAAACTGAAAGAAAACAACGAGGAAGGAA‑GCTTTTCGCGAGCCGGAACCCAGGGGGGGGGGAAAAAAAAGGGGGGGTTTTTTTCGccc                                                >  1:282435/1‑147 (MQ=255)                                                                                                         atTTGCCGGCTGGATACGGCGGGCGCACG‑AGGACTGGGAAAATAAAAAAAAATCTACAAACTGAAAAATCACAAAGAGCAGGGCA‑GCTTTTCGCAAGCCGCAAACCACGGGGTGGGGGAAAAACAGGGGGTGtttttttccccc                                                >  2:142012/1‑138 (MQ=255)                                                                                                            tGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCg                                       <  1:902653/150‑1 (MQ=255)                                                                                                             gCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCgg                                      <  2:54442/150‑1 (MQ=255)                                                                                                                     gATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGcc                             <  2:946337/151‑1 (MQ=255)                                                                                                                      aTACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGcc                             >  1:192737/1‑150 (MQ=255)                                                                                                                          ggcgggcgCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCctgact                        <  1:213680/151‑1 (MQ=255)                                                                                                                            cgggcgCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCctgactga                      <  2:578337/151‑1 (MQ=255)                                                                                                                               gcgcACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAAt                   <  2:443545/151‑1 (MQ=255)                                                                                                                                 gcACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAAtt                  <  2:691784/150‑1 (MQ=255)                                                                                                                                   aCG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATtgat               >  2:498810/1‑151 (MQ=255)                                                                                                                                   aCG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATtgat               <  1:515003/151‑1 (MQ=255)                                                                                                                                       aGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATtgatg              >  1:446100/1‑149 (MQ=255)                                                                                                                                       aGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATgag            <  1:881341/151‑1 (MQ=255)                                                                                                                                           cTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGc        >  2:989344/1‑151 (MQ=255)                                                                                                                                             ggAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGc        <  2:364543/149‑1 (MQ=255)                                                                                                                                             ggAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGc        <  2:869595/149‑1 (MQ=255)                                                                                                                                                aaaCTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAg   <  1:823828/151‑1 (MQ=255)                                                                                                                                                aaaCTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAg   <  1:141991/151‑1 (MQ=255)                                                                                                                                                 aaCTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGc  >  2:46315/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACG‑AGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCA‑GCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGC  >  NC_000913/3585‑3878

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