Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,895,122	0	A	T	22.9%	124.3 / 36.0	61	A106A (GCA→GCT)	pabB	aminodeoxychorismate synthase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base T (11/3);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 3.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATC  >  NC_000913/1895027‑1895218                                                                                                |                                                                                                 aCCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTgg                                                                                              >  2:293629/1‑100 (MQ=255) aCCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTgg                                                                                              >  1:22925/1‑100 (MQ=255)  ccACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTggg                                                                                             >  2:82821/1‑100 (MQ=255)   cACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTgggg                                                                                            <  2:1036647/100‑1 (MQ=255)   cACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTgggg                                                                                            <  1:789341/100‑1 (MQ=255)    aCTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGgt                                                                                           >  1:552653/1‑100 (MQ=255)     cTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGgtt                                                                                          >  1:177304/1‑100 (MQ=255)      tgatgaCCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCTCTGGGgttg                                                                                         <  1:707207/100‑1 (MQ=255)         tgaCCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGttt                                                                                      <  1:771436/100‑1 (MQ=255)                     ggTGCTCCAGCAGGTGCAGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTg                                                                          >  2:195928/1‑100 (MQ=255)                       tGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTggg                                                                        <  1:1279768/100‑1 (MQ=255)                       tGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTggg                                                                        >  2:398850/1‑100 (MQ=255)                       tGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTggg                                                                        <  2:241767/100‑1 (MQ=255)                        gCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGgc                                                                       >  2:526093/1‑100 (MQ=255)                         cTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGgcc                                                                      <  1:843840/100‑1 (MQ=255)                                 ggTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGtttt                                                              <  1:125350/100‑1 (MQ=255)                                  gTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCTCTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTg                                                             >  1:974825/1‑100 (MQ=255)                                  gTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCTCTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTg                                                             >  2:641678/1‑100 (MQ=255)                                  gTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTg                                                             <  1:193700/100‑1 (MQ=255)                                    gCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCTCTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                           <  2:58110/100‑1 (MQ=255)                                    gCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAg                                                           <  2:627328/100‑1 (MQ=255)                                     cTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGt                                                          >  1:83046/1‑100 (MQ=255)                                      tGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTc                                                         >  1:1148285/1‑100 (MQ=255)                                         aTCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCTCTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACt                                                      >  1:12456/1‑100 (MQ=255)                                         aTCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCTCTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACt                                                      >  1:12588/1‑100 (MQ=255)                                         aTCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACt                                                      >  2:150091/1‑100 (MQ=255)                                         aTCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACt                                                      <  1:950323/100‑1 (MQ=255)                                          tCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTg                                                     <  1:237170/100‑1 (MQ=255)                                           cgcgCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGc                                                    <  1:990164/100‑1 (MQ=255)                                           cgcgCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGc                                                    <  2:3024/100‑1 (MQ=255)                                            gcgcAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGcc                                                   >  1:534043/1‑100 (MQ=255)                                             cgcAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCTCTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCa                                                  >  2:1033679/1‑100 (MQ=255)                                             cgcAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCa                                                  >  2:28885/1‑100 (MQ=255)                                               cAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCTCTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGa                                                >  1:305383/1‑100 (MQ=255)                                               cAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGTCGTTTTGAGTCACTGCCAGa                                                <  2:799467/100‑1 (MQ=255)                                               cAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGa                                                >  1:1224878/1‑100 (MQ=255)                                                   cATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAAtt                                            >  2:912611/1‑100 (MQ=255)                                                   cATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAAtt                                            <  2:840372/100‑1 (MQ=255)                                                       cGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCATTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCgg                                        <  2:913505/100‑1 (MQ=255)                                                         cacAACGCATAACTAAGATTTGCCATGTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCgg                                        <  1:566068/96‑1 (MQ=255)                                                            aaCGCATAACGAAGATTTGCCATTTCAGGGCGGCGCTCTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACaa                                   <  1:999602/100‑1 (MQ=255)                                                             aCGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAg                                  >  1:394258/1‑100 (MQ=255)                                                                cATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACACGat                                >  2:524968/1‑99 (MQ=255)                                                                cATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGata                               >  2:599855/1‑100 (MQ=255)                                                                     cGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGtt                          <  1:644231/100‑1 (MQ=255)                                                                     cGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGtt                          >  1:640375/1‑100 (MQ=255)                                                                        aGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTg                       <  2:331036/100‑1 (MQ=255)                                                                         gATTTGCCATTTCAGGGCGGCGCTCTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGc                      >  1:764597/1‑100 (MQ=255)                                                                          aTTTGCCATTTCAGGGCGGCGCTCTGGGGTTGTTTGGCTACGATCTGGGCAGCCGATTTGAGTCACTGCCAGAGACTGAGGAACAAGATATAGTTATGcc                     >  1:849163/1‑100 (MQ=255)                                                                           tttGCCATTTCAGGGCGGCGCTCTGGGGTTGTTTGGCTACGATCTGGGCCGGCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCg                    >  1:32710/1‑100 (MQ=255)                                                                            ttGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCgg                   >  1:353602/1‑100 (MQ=255)                                                                              gCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGat                 <  1:622890/100‑1 (MQ=255)                                                                                  tttCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATgg             <  2:183958/100‑1 (MQ=255)                                                                                   ttCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGc            <  2:22925/100‑1 (MQ=255)                                                                                     cAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAg          <  1:1129935/100‑1 (MQ=255)                                                                                       gggCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCCGAAATTGCGGCACAAGATATCGTTCTGCCGGATATGGCAgag        <  1:878860/100‑3 (MQ=255)                                                                                          cggcgCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGt     >  2:5056/1‑100 (MQ=255)                                                                                           ggcgCTCTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCTGATATTGCGGAACAAGCTCTCGTTCTCCCGGATATGGCCGTTGGt     >  2:835338/1‑99 (MQ=255)                                                                                           ggcgCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTa    <  1:192562/100‑1 (MQ=255)                                                                                            gcgcTCTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTAt   >  1:822964/1‑100 (MQ=255)                                                                                             cgcACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATc  <  2:83044/100‑1 (MQ=255)                                                                                                |                                                                                                 ACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCCCAACGCATAACGAAGATTTGCCATTTCAGGGCGGCGCACTGGGGTTGTTTGGCTACGATCTGGGCCGCCGTTTTGAGTCACTGCCAGAAATTGCGGAACAAGATATCGTTCTGCCGGATATGGCAGTGGGTATC  >  NC_000913/1895027‑1895218

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑