Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,327,934	0	A	G	43.4%	63.6 / 23.6	53	D267G (GAT→GGT)	crfC	clamp‑binding sister replication fork colocalization protein, dynamin‑related
Reads supporting (aligned to +/- strand):  ref base A (4/26);  new base G (23/0);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 3.23e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATT  >  NC_000913/4327791‑4328082                                                                                                                                                |                                                                                                                                                         tGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TAc                                                                                                                                                   <  1:371278/150‑1 (MQ=255)   tCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGGGCTGGGTTT‑TACGc                                                                                                                                                 >  2:587143/1‑150 (MQ=255)    ccATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGATGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGtttt                                                                                                                                                      >  2:590311/1‑143 (MQ=255)      aTCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGtt                                                                                                                                                        >  2:379303/1‑141 (MQ=255)      aTCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCCGGGGTA‑TACGCAAc                                                                                                                                              >  2:913538/1‑150 (MQ=255)         tAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGGGCTGGGTta‑ttcacca                                                                                                                                               >  2:247156/1‑140 (MQ=255)          aGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑Gaaa                                                                                                                                        <  1:254328/151‑1 (MQ=255)          aGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑Ga                                                                                                                                          <  1:254898/149‑1 (MQ=255)           gCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAAt                                                                                                                                       <  2:791887/151‑1 (MQ=255)             ggggCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCg                                                                                                                                     <  1:244690/151‑1 (MQ=255)             ggggCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGCGGTGGAAta‑tattcacacgttaat                                                                                                                                       >  2:620473/1‑137 (MQ=255)                      gagTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGAt                                                                                                                            <  1:16524/151‑1 (MQ=255)                        gTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCg                                                                                                                              <  2:598842/147‑1 (MQ=255)                           aTCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAag                                                                                                                        <  2:704561/150‑1 (MQ=255)                               cGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTc                                                                                                                   <  1:199361/151‑1 (MQ=255)                                    aGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGTTT‑TAAGCAACT‑GAAAACGATCTCCCGTGAAGAGGGCC‑GTg                                                                                                              >  2:241465/1‑150 (MQ=255)                                     gTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGTTT‑TTAGCAACT‑GAAAACGAACTCCCCTGAAGAGGGcc                                                                                                                  >  2:154090/1‑146 (MQ=255)                                       tGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGTTA‑TACGCAAAT‑GAAAACGATCTCCGATGGAGAggggc                                                                                                                  >  2:666648/1‑141 (MQ=255)                                          cGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGGGCTGGGTTT‑TAAGCAAAC‑TAAAACGATCTCCGCTGGAAAGGTCC‑gggggg                                                                                                           >  1:864830/1‑142 (MQ=255)                                              aCTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGTTA‑TACGCAACT‑GAAAACGATCTCCCATGAAAAGGGCC‑GTGGGGGGGtttt                                                                                                    >  2:338482/1‑150 (MQ=255)                                              aCTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGGGCCGGGGTA‑TAAACCACT‑GAAAACCCACTCCCATGAAGAGGACC‑GCGGGGCGAtttt                                                                                                    >  1:949273/1‑150 (MQ=255)                                                  gATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGTTA‑TAAGCCACT‑GAAAACGATCTCCGCAGAAGAGGGCC‑GGGAGGGGGTTTTTGC‑gg                                                                                              >  2:885046/1‑151 (MQ=255)                                                   aTACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGt                                                                                             <  2:725961/151‑1 (MQ=255)                                                      cccccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGGGCTGGATTA‑TAAACAAATCGAAACCG‑TCTCCCCTGAAGAGGTTC‑GTGAGGGGGTTTTTGC‑gggggg                                                                                          >  1:661598/1‑151 (MQ=255)                                                            ggCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGGGCTGGGTTA‑TAACCAACT‑TAAAACGGTCTCCCCTGAAGAGGGCC‑GTGGGGGGGTTTTTGC‑Ggggggg                                                                                         >  1:325163/1‑146 (MQ=255)                                                                    gAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGGGGTGGGTTA‑TACGCAACT‑GAAAACGCACTCCGGTGAAGAGGGCC‑GGGGGGGGATTTTTGG‑GGGGGGGGAAACgg                                                                                  >  2:623496/1‑145 (MQ=255)                                                                       gCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGTTA‑TACGCAACT‑GAAAAAGATCTCCGCTGAAGAGGGCC‑GGGAGGGGGGTTTTGC‑GGGGGGGGaaaagg                                                                                  >  2:773109/1‑138 (MQ=255)                                                                       gCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTg                                                                          <  2:849771/150‑1 (MQ=255)                                                                        ccGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGTTT‑TAAGCCACT‑GAAAACGATCTCCCATGGAGAGGGGC‑GGGGGGGGATTTTGGG‑GGGGGGGGGATCGGGGCCGCCGtg                                                                        >  1:500228/1‑150 (MQ=255)                                                                                ccGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTTGGTTA‑TAAGCAACT‑TAAAACGATCTCCCGTGGAGAAGGTC‑GTGAGGCGGATTTTGC‑GGGGGGGGAATAGGGGGCGCttgaacgtctgg                                                                >  2:568074/1‑140 (MQ=255)                                                                                 cGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGt                                                               <  1:174568/151‑1 (MQ=255)                                                                                        gCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGTTA‑TACGCCACT‑GAAAACGATCTCCGATGAAGAGGGCC‑GTGGGGCGGTTTTTGG‑GGGGGGGGAATCGGGGCCGCTGTACGTGCTGGGCAATaa                                                         >  2:309080/1‑150 (MQ=255)                                                                                        gCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGTTA‑TAAGCCACT‑GAAATCGGACTCCGATGAAGAGGGCCGGTGGGGGGTTTTTGGG‑GGGGGGGC‑ATCGGGGGCGCTGTACGGGGTgggtaaaaa                                                         >  2:914868/1‑143 (MQ=255)                                                                                        gCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAg                                                        <  2:359784/151‑1 (MQ=255)                                                                                        gCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAg                                                        <  1:462619/151‑1 (MQ=255)                                                                                            aaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCg                                                    <  1:416010/151‑1 (MQ=255)                                                                                            aaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCg                                                    <  1:587046/151‑1 (MQ=255)                                                                                              aaTGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTTGGTTA‑TAACCCAAT‑GAAAAAGATCTTCCGTGGAGAGGGCC‑GGGGGGGGGTTTTTGCGGGGGGGGCAACCGGGGCGGC‑GTCCGGGCTGGTccaaaaaatt                                                     >  1:857537/1‑139 (MQ=255)                                                                                                       ccAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGAt                                         <  1:929474/151‑1 (MQ=255)                                                                                                           cagcTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGTTA‑TACGCCAAT‑GAAAACGATCTCCGATGAAGAGGGCC‑GTGGGGGGAGTTTTGC‑GGGGGGGGAAACGGGGGCGCTGTACGGGCCGGGCCATAAGtttggttaaaa                                             >  2:881907/1‑134 (MQ=255)                                                                                                           cagcTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTa                                     <  2:9494/151‑1 (MQ=255)                                                                                                             gcTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTaa                                    <  2:361395/150‑1 (MQ=255)                                                                                                               tGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGGTTA‑TACGCAAAT‑GAAAAAGAACTCCCATGGAAAGGGCC‑GGGGGGGGATTTTTGG‑GGGGGGGGAAACGGGGGCCCCGTTAGTGCTGGGCAAAAAATTTGGTCaaaaaggacggaa                                    >  2:475710/1‑137 (MQ=255)                                                                                                               tGGCACGCGCCTCGGCGGTACTGGCGGTGCCGGGTTATTAGCCAATT‑GAAA‑CGATCTCCCGTGGAGAGGGGC‑GGGGGGGGGTTTTTGG‑GGGGGGGGGAACGGTGGCGCCGTGCGGGGTTGGCAAAAAATTTGTTCAACAgggacggaaaa                                  >  1:798849/1‑141 (MQ=255)                                                                                                                ggCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACa                                  <  1:570049/149‑1 (MQ=255)                                                                                                                  cACGCGCTTCGGCGGTACTGGCGATGCTGGATTA‑TACGCAGCT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTga                              <  1:118565/151‑1 (MQ=255)                                                                                                                              cggTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCa                  <  2:760478/151‑1 (MQ=255)                                                                                                                               ggTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAg                 <  1:289731/151‑1 (MQ=255)                                                                                                                                   cTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGc             >  1:910089/1‑151 (MQ=255)                                                                                                                                     ggCGGTGCTGGATTA‑TACGCAACT‑GAAAGCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCg            <  2:485158/150‑1 (MQ=255)                                                                                                                                       cGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGc         >  1:154228/1‑151 (MQ=255)                                                                                                                                        ggTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCa        <  2:8925/151‑1 (MQ=255)                                                                                                                                              ggATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGAGCAATCGGTGCCGCTGTACGTGCTGGTCCATAAGTTCGATCAACATGATCGTAACAGTGACGACGCCGACCAGGTGGGGGCCTTGAtt  <  1:518205/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                         TGTCCATCTAGCGGGGCTGGAGAGTTATCCCGGTCAGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTA‑TACGCAACT‑GAAATCGATCTCCGATGAAGAGGTCC‑GTGAGGCGATTTTGGC‑GGTGGGGCAATCGGTGCCGCTGTACGTGCTGGTCAATAAGTTCGATCAACAGGATCGTAACAGTGACGACGCCGACCAGGTGCGGGCACTGATT  >  NC_000913/4327791‑4328082

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