Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,995,457	0	T	G	41.8%	112.4 / 15.7	68	G232G (GGT→GGG)	dapF	diaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base T (7/32);  new base G (28/0);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.15e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAA  >  NC_000913/3995307‑3995603                                                                                                                                                       |                                                                                                                                                    cGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGt                                                                                                                                                     <  2:1077069/151‑1 (MQ=255)     cGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGc                                                                                                                                                 <  1:608165/151‑1 (MQ=255)        tttCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGGTTGCt                                                                                                                                                >  2:633264/1‑149 (MQ=255)              gAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTTAGCAAGGGTTTCTGGcc                                                                                                                                            >  1:541796/1‑147 (MQ=255)                gcgcgcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCgaagaa                                                                                                                                      <  1:1066434/151‑1 (MQ=255)                gcgcgcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCgaagaa                                                                                                                                      <  2:564763/151‑1 (MQ=255)                gcgcgcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTTCTGGCCgaaga                                                                                                                                       >  2:1073964/1‑150 (MQ=255)                  gcgcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGTTGGGATTCAGCAAAGGTTGGTGGCCggagaa                                                                                                                                      >  1:105869/1‑149 (MQ=255)                  gcgcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGt                                                                                                                                    <  1:250159/151‑1 (MQ=255)                    gcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTAc                                                                                                                                  <  2:413624/151‑1 (MQ=255)                     ccAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGCGGAAGGGTTTAGCAAGGTTTGCTGGCCggagaaag                                                                                                                                    >  1:1215934/1‑146 (MQ=255)                      cAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCTGTTGGGGCTCAACGAAGGTTGCCGGGCCAAAAAGTat                                                                                                                                  >  2:148521/1‑148 (MQ=255)                      cAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGGTTGCTGGCCGaagaagga                                                                                                                                   >  1:892081/1‑146 (MQ=255)                       aaTATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTTAGCAAGGGTTGCTGGCCCGAGAAGGAcgcg                                                                                                                               >  1:1263103/1‑151 (MQ=255)                        gtatCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGGGTGCGGCGGTTTCGGGATGGGGTTAGGAAGGATTTCCGGGCGAAGGAGTAccggg                                                                                                                              >  1:274535/2‑147 (MQ=255)                          atCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTg                                                                                                                             <  2:151598/150‑1 (MQ=255)                            cGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGaa                                                                                                                          <  1:1271808/151‑1 (MQ=255)                             ggTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGTAGGGGTTCAGCCAGGGTTTCTGGCCGGAGAAGGAAGCGTgg                                                                                                                            >  1:1194659/1‑148 (MQ=255)                             ggTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGGAGGGGTTCAGCAAGGGTTGGTGGCCGGAGAAGGACGCgggg                                                                                                                            >  1:158315/1‑148 (MQ=255)                             ggTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAAc                                                                                                                         <  2:571446/151‑1 (MQ=255)                                tttATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTc                                                                                                                       <  2:1085215/150‑1 (MQ=255)                                    tGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGTTCAGCAAGGGTTGCTGGCCGAAGAAGTACGCGGGGGACTcccccg                                                                                                                  >  2:380103/1‑149 (MQ=255)                                     gCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCcggc                                                                                                                 <  2:1279678/151‑1 (MQ=255)                                         gTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCCAGGGTTGCTGGCCGAAGAAGTACGCGGGGGACTCCCcggcgg                                                                                                               >  2:403121/1‑149 (MQ=255)                                          tGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAATGGGGATTTCGCCAGGGTTGATGGCCGGAGAAGAACGCGGGGGACTcccccgggg                                                                                                               >  1:1296094/1‑143 (MQ=255)                                           ggTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGgtcgt                                                                                                           <  2:1252095/151‑1 (MQ=255)                                            gTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGgtcgtc                                                                                                          <  2:278275/151‑1 (MQ=255)                                              tAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGAGGAGGGGTTCAGCAAGGGTTTCTTGCCGAAGAAGGAAGCGTGGGAATCCCCCGCGGGCGTCtt                                                                                                        >  1:59981/1‑151 (MQ=255)                                                  cgcgAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGata                                                                                                    <  2:1108011/151‑1 (MQ=255)                                                   gcgAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGatat                                                                                                   <  2:990260/151‑1 (MQ=255)                                                        cATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGAATCAGCGAGGGTTGCCGGCCGAAGAAGGACGCGGGGGAGTCCCCCGGGGTGGTCTTTCTATccc                                                                                                >  1:383201/1‑147 (MQ=255)                                                         atatTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTGCTGGGCGAAGAAGGACGCGGGGAAATCCCCGGCGGGCGGCTTGGTATTG‑cc                                                                                              >  2:880348/1‑149 (MQ=255)                                                           atTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGTTTAGCAAGGTTTGCTGGCCGGAGAAGGACGCGGGGGACTCCCCGGCGGGCGTCTTTGTTTCG‑cc                                                                                              >  1:825466/1‑147 (MQ=255)                                                              cGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGaaa                                                                                        <  1:767189/150‑1 (MQ=255)                                                               gTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTTCAGTTGGGATTCAGCAAGGGTTGCTGGCCGAAGAAAGAAGCGGGGGACTCCCCGGCGGGCGCCTTGTTATCG‑CCCGGgaa                                                                                        >  2:538143/1‑149 (MQ=255)                                                               gTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTAAGCGTGGAACTCCCCGGCGGGCGGCTTGATATTG‑CCCGGAAAgg                                                                                      >  2:211098/1‑151 (MQ=255)                                                               gTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGAATCAGCAAGGGTTGCTGGCCGAAGAAGGACGCGTGGAACTCCCCCGCGGGCGGCTTGATATCG‑CCTGGaaa                                                                                        >  1:568231/1‑149 (MQ=255)                                                                  tACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGAGGAGGGATTCAGCAAGGGTTGCTGGCCCGAGAAAGACGCGGGGGACTCCCCCGCGGGCCTCTTTATATCG‑CCCGGGAAGGGcc                                                                                   >  1:427686/1‑151 (MQ=255)                                                                  tACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGCTGCAGTAGGGAGTTAGCAAGGGTTGCTGGCCGCAAAAGAACGCGTGGAACTCCCCCGCGGGCGTCTTTTTATTG‑CCTGGgaacg                                                                                      >  1:823377/1‑146 (MQ=255)                                                                      cgTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGt                                                                               <  1:967603/151‑1 (MQ=255)                                                                         ttATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGAGGAGGGGGTCAGCAAGGGTTTGTGGCCGAAGAAGGACGCGGGGAACTCCCCCGGGGGCGGCTTTGTATTG‑CCCGGGAAAGGCCgggggac                                                                            >  1:745162/1‑147 (MQ=255)                                                                         ttATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCa                                                                             <  1:67058/150‑1 (MQ=255)                                                                          tATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTGCTGGCCGAAGAAAGAAGCGTGGGACTCCCCCGCGGGCGGCTTTTTATCG‑CCTGGGAAAGGCCggggcccc                                                                           >  1:1153126/1‑146 (MQ=255)                                                                           aTGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTTCTGGCCGAAGAAGTACGCGTGGGACTCCCCGGCGGGCGGCTTGATATCG‑CCTGGGAAGGGCCGGGTCAccc                                                                          >  2:314492/1‑151 (MQ=255)                                                                               gCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTta                                                                      <  1:764438/151‑1 (MQ=255)                                                                                   ggggCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGAGTAGGGATTTAGCAAGGGTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCCGCGGGCGTCTTGGTATTG‑CCTGGAAAGGTCCGGGGCACCCGtt                                                                       >  1:515157/1‑146 (MQ=255)                                                                                       catgagaAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTTCTGGCCGGAGAAATACGCGGGGGAATCCCCCGCGGGCGTCTTTGTATCG‑CCCGGGAAGGGCCCGGGCACCCGTTatatat                                                                 >  2:238638/4‑148 (MQ=255)                                                                                       cAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGAGGAGGGATTCAGCAAGGGTTTCTGGCCCGAGAAGGACGCGGGGGACTCCCCCGCGGGCGGCTTGGTTTCGCCCTGGAAAGG‑CCGGGGCAccccgt                                                                       >  1:721820/1‑139 (MQ=255)                                                                                        aGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGTTAAACAAGGGTTGCTGGCCCAAGAAGGACGCGTGGGAACCCCCGGCGGGCGTCTTTATATTG‑CCAGGGAAAGTCCCGGTTGCCCGTtttattt                                                                 >  2:565219/1‑145 (MQ=255)                                                                                         ggAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTg                                                            <  2:943055/151‑1 (MQ=255)                                                                                         ggAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGGTTGCTGGCCGAAGAAGGAAGCGGGGGACTCCCCGGCGGTCGTCTTTATATTG‑CCTGGAAAGGGCCGGGGCCCCCCTtatat                                                                   >  2:1167960/1‑144 (MQ=255)                                                                                               aCCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCcgg                                                      <  1:252157/151‑1 (MQ=255)                                                                                                cccAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCcggc                                                     <  2:826459/151‑1 (MQ=255)                                                                                                      ccTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACa                                               <  2:24577/151‑1 (MQ=255)                                                                                                         gCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACAt                                              <  1:519301/149‑1 (MQ=255)                                                                                                         gCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATg                                             <  1:245897/150‑1 (MQ=255)                                                                                                             cAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTac                                        <  2:1119991/151‑1 (MQ=255)                                                                                                                 ggcgCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAAGGTTGCTGGCCGAAGAAGGACGCGTGGAACTCCCCGGCGGGCGTCTTTATATCG‑CCCGGAAAAGGCCGGGGCACCCGTTATATATTACTGGCCCGGCGGGACAAGTTTacgac                                     >  2:557087/1‑150 (MQ=255)                                                                                                                  gcgcGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACgg                                   <  1:1225762/151‑1 (MQ=255)                                                                                                                    gcgTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGAt                                 <  2:2357/151‑1 (MQ=255)                                                                                                                      gtgtGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGGCACCCGTTATATATGACTGGCCCGGCGGGACATGTCTACGACGGAttt                               >  2:1059834/1‑151 (MQ=255)                                                                                                                          gcggcgGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTAtt                            <  1:547952/150‑1 (MQ=255)                                                                                                                              cggTTGCAGTAGGGATTCGGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCAt                         >  2:689077/1‑149 (MQ=255)                                                                                                                              cggTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCt                       <  2:576225/151‑1 (MQ=255)                                                                                                                                  tGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATg                    >  2:151053/1‑150 (MQ=255)                                                                                                                                    cAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCTGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGa                   <  1:613667/149‑1 (MQ=255)                                                                                                                                    cAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGa                   >  1:323021/1‑149 (MQ=255)                                                                                                                                              tCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAggg        <  1:674059/150‑1 (MQ=255)                                                                                                                                                     ggTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGgaagaa  <  2:833427/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                    CGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCG‑CCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGAAGAA  >  NC_000913/3995307‑3995603

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: ‑