Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,564	0	A	C	44.4%	122.0 / 22.5	72	*559E (TAA→GAA)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (35/5);  new base C (0/32);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.55e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCT  >  NC_000913/501418‑501710                                                                                                                                                   |                                                                                                                                                    tGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑cccc                                                                                                                                                >  1:716013/1‑151 (MQ=255)    gATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTc                                                                                                                                             <  1:735193/151‑1 (MQ=255)       cGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTcacc                                                                                                                                          >  1:329091/1‑151 (MQ=255)        gTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTcacc                                                                                                                                          >  2:310475/1‑150 (MQ=255)           aaGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCt                                                                                                                                      <  1:450591/151‑1 (MQ=255)                    tcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACccgcc                                                                                                                              <  1:1117573/150‑1 (MQ=255)                      ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCgg                                                                                                                            <  2:1237923/150‑1 (MQ=255)                              ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                                   >  1:421622/1‑151 (MQ=255)                                ccAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTcc                                                                                                                  >  1:884091/1‑150 (MQ=255)                                 cAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTcc                                                                                                                  >  2:90871/1‑149 (MQ=255)                                  aGATAGGGCGTTCCCGCCGCATTTGGTTTTTGGCTCTCGATCCCTGATGCGCCGCTGGCGCGTCTTATCATCCCTACATTTTTTTCATATTTTCCACCCGGAACCCCCCTTTTC‑CCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTcc                                                                                                                  <  1:1119601/148‑1 (MQ=255)                                  aGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                               >  1:1094251/1‑151 (MQ=255)                                   gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                               >  2:1105609/1‑150 (MQ=255)                                     tAAGGCTTCCCCCCCTCACCTGTCATTTGGCTCTCGATGCCTGATGCCCCGCGGGCGCGTCTTATCATGCCTACATTTTTTTCTTATTTTCCACCCGGCACCCCCCTTTTC‑CCCCGCCCCCCCCACCCCCGCCGGTGGCGTTTCca                                                                                                                 <  1:683391/146‑1 (MQ=255)                                     tAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGt                                                                                                            >  2:603859/1‑151 (MQ=255)                                     tAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGt                                                                                                            >  1:70047/1‑151 (MQ=255)                                      aaGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGtt                                                                                                           >  1:700243/1‑151 (MQ=255)                                       aGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTc                                                                                                          >  2:587561/1‑151 (MQ=255)                                                cgccgcTTCGGGCTTTTGCCTCTCGATGCCTGATGCGACGCTGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACCCCCGCCACCCCCCGTTTC‑CCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                                                                                    <  1:296140/148‑1 (MQ=255)                                                cgccccATCTGGATTTTGGCTCTCGATCCCTGATGCGCCGCGGGCGCTTTTTTTCATCCCTACATATTTTTCATTTTTTACACCCGCCACCCCCCTTTTC‑CCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                 <  1:178236/151‑1 (MQ=255)                                                  ccgcATCTGGCTTTTGGCTCTCGATGCCTGAGGCGACGCTGGCGCTTTTTATCATGCCTAATTTTTTTTCATATTTTACACCCGCCACCCCCCTTTTC‑CCCCTCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                                  <  2:69116/148‑1 (MQ=255)                                                  ccgcATCTGGCTTTTGCCTCTCGATGCCTGTGGCAACGCGGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACCTCCGGCACCCCCTGTTTC‑CCCCTTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                 <  1:382616/149‑1 (MQ=255)                                                  ccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTAc                                                                                               >  2:464999/1‑151 (MQ=255)                                                   cgcATCTGGTTTTTGGCTCTCGATCCCTGATGCGCCGCGGGCGCTTTTTATCATCCCTAATTTTTTTTCATTTTTTCCACCGGGCACCCCCCGTTTC‑CCCCTCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTa                                                                                                <  1:782838/149‑1 (MQ=255)                                                   cccttctgctttttGGCTCTCGATGCCTGATGCGCCGCGGGCGCGTTTTTTCATCCCTACATATTTTTCATATTTTACCCCCGGCACCCCCCTTTTC‑CCCCTCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGc                                                                                                      <  1:477340/132‑1 (MQ=255)                                                     cATCTGGCTTTTGGCTCTCGATGCCTGAGGCGACGCTGGCGCTTCTTATCATCCCTACATATTTTTCATATTTTACATCCGGCAACCCCCGTTTC‑CCCCGTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACgg                                                                                             <  2:924184/150‑1 (MQ=255)                                                        cTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCg                                                                                          >  1:668936/1‑150 (MQ=255)                                                              tttGGCTCTCGATGCCTGATGCGACGCGGGGGCTTCTTTTCATGCTTACTTTTTTTTCATTTTTTACACCCGGCACCCACCTTTTC‑CCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAttt                                                                                      <  1:981847/148‑1 (MQ=255)                                                                tGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTTTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAACCACCGTTTC‑CCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACg                                                                                  <  1:424366/150‑1 (MQ=255)                                                                   ctcccgagcccTGTTGCGCCGCGGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGCCACCCCCCTTTTC‑CCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                                <  1:520281/140‑1 (MQ=255)                                                                       cGATCCCTGAGGCGACGCGGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                             <  1:1050816/148‑1 (MQ=255)                                                                             cTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                     >  2:770464/1‑150 (MQ=255)                                                                               gATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                                                                  >  2:69399/1‑151 (MQ=255)                                                                                aTGCGCCGCTGGCGCTTCTTATCATCCCTACATATTTTCCATTTTTTCCACCCGGCACCCCCCTTTTC‑CCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                     <  2:1042315/147‑1 (MQ=255)                                                                                    gccGCTGGCGCTTTTTATCTTCCCTACATTTTTTTCTTATTTTCCATCCGGCAACCCCCGTTTC‑CCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                              <  2:1142080/148‑1 (MQ=255)                                                                                     aCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                            >  1:661883/1‑151 (MQ=255)                                                                                     aCGCTGGCGCGTCTTATCATCCCTACATATTTTTCATTTTTTACACCCGGCAACCCCCGTTTC‑CCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                             <  1:699131/150‑1 (MQ=255)                                                                                         tGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac                                                        >  1:1026880/1‑151 (MQ=255)                                                                                          ggCGCGTTTTATCATGCCTACATTTTTTTCATTTTTTCCCCCCGGCACCCCCCTTTTC‑CCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGca                                                         <  2:228275/149‑1 (MQ=255)                                                                                            cgcgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCTCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcaccac                                                     >  2:833707/1‑151 (MQ=255)                                                                                             gcgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg                                                    >  1:871390/1‑151 (MQ=255)                                                                                              cgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt                                                   >  1:425747/1‑151 (MQ=255)                                                                                               gTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                  >  1:1249296/1‑151 (MQ=255)                                                                                                 cTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                >  2:37797/1‑151 (MQ=255)                                                                                                       aTGCCTACATTTTTTTCATATTTTACACCCGGCAACCCCCTTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                           <  2:120253/150‑1 (MQ=255)                                                                                                       aTGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                           >  2:573737/1‑150 (MQ=255)                                                                                                         gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                         >  2:561460/1‑150 (MQ=255)                                                                                                          ccTACATTTTTTCCATATTTTACACCCGCAAACCCCCTTTTC‑CCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGa                                            <  1:1170347/146‑1 (MQ=255)                                                                                                             acttttttttCATATTTTCCACCCGCCAACCCCCTTTTC‑CCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                               <  1:1025645/135‑1 (MQ=255)                                                                                                                tttttttCATTTTTTACACCCGGCACCCACCTTTTC‑CCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCAc                                  <  2:255181/148‑1 (MQ=255)                                                                                                                tatTTTTCATATTTTACCCCCGCCACCCCCCGTTTC‑CCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCAc                                  <  1:768635/150‑1 (MQ=255)                                                                                                                 atTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                                >  1:721295/1‑151 (MQ=255)                                                                                                                  tttttCATATTTTACACCCGGCAACCACCGTTTC‑CCCCCCAACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACt                                 <  2:87692/149‑1 (MQ=255)                                                                                                                  tttttCATATTTTACACCCGGCAACCACCGTTTA‑CCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACt                                 <  2:23877/149‑1 (MQ=255)                                                                                                                    ttttATTTTTTCCATCCGCCACCCCCCGTTCC‑CCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                              <  1:750223/150‑1 (MQ=255)                                                                                                                        attttttACACCCGGCACCCCCCGTTTC‑CCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATc                            <  1:748104/145‑1 (MQ=255)                                                                                                                          atTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                       >  2:1091537/1‑151 (MQ=255)                                                                                                                           ttttACACCCGCCACCCCCCGTTTC‑CCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                       <  1:1177645/150‑1 (MQ=255)                                                                                                                            tttACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                      >  2:689990/1‑150 (MQ=255)                                                                                                                            tttACACCCGCCACCCCCCTTTTC‑CCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                      <  1:630519/150‑1 (MQ=255)                                                                                                                                cATCCGGCACCCCCCTTTTC‑CCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                    <  1:409689/148‑1 (MQ=255)                                                                                                                                cATCCGGCACCCCCCGTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCACACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                  <  2:884104/150‑1 (MQ=255)                                                                                                                                cATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATggg                 >  1:908784/1‑151 (MQ=255)                                                                                                                                 cacccgCAACCCCCCTTTCCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                    <  2:224459/143‑1 (MQ=255)                                                                                                                                 aTCCGGCACCCCCCGTTTC‑CCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                    <  2:37781/147‑1 (MQ=255)                                                                                                                                    cGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg             >  2:605062/1‑151 (MQ=255)                                                                                                                                     ggCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg            >  2:1246497/1‑151 (MQ=255)                                                                                                                                          cccccGTTTC‑CCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGAGCaa        <  1:11844/147‑1 (MQ=255)                                                                                                                                           caccGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaat       >  2:41539/1‑150 (MQ=255)                                                                                                                                           caccGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaat       >  1:1128005/1‑150 (MQ=255)                                                                                                                                             ccGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat    >  2:1091566/1‑151 (MQ=255)                                                                                                                                               gTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct  >  2:503229/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTA‑CCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCT  >  NC_000913/501418‑501710

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