Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,575	0	A	C	55.2%	67.5 / 93.4	58	V555G (GTG→GGG)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (24/2);  new base C (0/32);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 2.53e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-11

AAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGG  >  NC_000913/501428‑501720                                                                                                                                                     |                                                                                                                                                  aaGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCaccacc                                                                                                                                                 >  2:333860/1‑150 (MQ=255) aaGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCt                                                                                                                                                >  1:747652/1‑151 (MQ=255)        cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACccgc                                                                                                                                         >  2:547490/1‑151 (MQ=255)                    ggCCAGGTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCCCCGGTGGCGTTTc                                                                                                                             >  2:963890/1‑151 (MQ=255)                        aGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                                         >  2:351460/1‑151 (MQ=255)                           tAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGt                                                                                                                      >  2:1002744/1‑151 (MQ=255)                                  tccccgccgcATCTGGCTTTTGCCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATCCCTACATTTTTTTCATATTTTACACCCGGCAC‑CCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCa                                                                                                                  <  1:182538/144‑1 (MQ=37)                                     ccgccgcATTGGGTATTTGCCTCCATCTGCCTGATGCAGCGCGGCCGCATCTTATCATGCCACCATTTTTTTTAT‑TTTTACCTCCGGGAACCCACGTTTTCCCCCGTCCCCCCCACACCCGCCGGTGGCGTTTCCAGCAGTTcc                                                                                                                   <  2:520907/143‑1 (MQ=255)                                     aCGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                                            >  2:241986/1‑151 (MQ=255)                                     aCGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                                            >  1:811759/1‑151 (MQ=255)                                        ccgttgcTGGCATTTCTCTCTCGATACCGGTTGCGACGGTGGCGGGTCTTAACAAGCCTCCATATTTTACATATTTTTCATCCGGCAA‑CCTCCTTTTCCCCCGTCCCCCCCTCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                           <  1:1198479/143‑1 (MQ=255)                                        ccgcATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTCTTATCATGCCTACTTTTTTTTCATTTTTTACATCCGGCAA‑CCCCCTTTTCCCCCGTCACCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTAc                                                                                                         <  1:447647/151‑1 (MQ=38)                                           cATCTGTCTTTTGTCCCTCGTCGCCGGAGGCGCCGCCGGCGCGTTTTCTCGTCCCTAATTTTTTTTCATTTTTTTCATCCGCCAA‑CCACTGTTGCCCCTGCCGCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg                                                                                                            <  2:1201909/145‑1 (MQ=255)                                             tCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTTTTATCATCCCTACATTTTTTTCATTTTTTCCACCCGGCAA‑CCCCCTTTTACCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                                                                                                      <  2:1170954/149‑1 (MQ=255)                                              cTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTAACCCGTCACCACCTCACCCGCCGGCGGCGTTTCCAGCAGATCCAGCATGGTACGGGCg                                                                                                    >  2:792507/1‑150 (MQ=255)                                               tGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                                  >  1:362765/1‑151 (MQ=255)                                                ggCTTTTGGCTCTCTATCCTGGATGCGACGCTGGCCCGTTTTATCATGCCTACTTTTTTTCCATTTTTTACACCCGCCAC‑CCCCCGTTACCCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                                  <  2:79697/150‑1 (MQ=255)                                                ggCTTTTGGCTCTCGATGCCTGATGCGACGCGGGCGCTTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCAC‑CCACCTTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                                  <  1:845602/150‑1 (MQ=255)                                                ggCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                                                                                                 >  1:1224315/1‑151 (MQ=255)                                                    tttGTCTCCCGATGCCTTACGCGACGCGGGCGCGTCTTATCTGCCCTACATATTTTTAATATTTTACACCCGCCAA‑CCACCGTTTCCCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAttt                                                                                                <  1:1025726/148‑1 (MQ=255)                                                    tttGGCTCTCGTTGCCTGATGCGCCGCGGGCGCGTCTTATCATGCCTACATTTTTTTCATTTTTTCCACCCGGCAC‑CCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCa                                                                                              <  1:125289/150‑1 (MQ=255)                                                    tttGGCTCTCGATGCCTGATGCGCCGCGGGCGCTTCTTTTCATCCCTACATTTTTTTCATTTTTTACACCCGGCAC‑CCCCCTTTTCCCCCGCCCCCCCCGCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                               <  2:11358/149‑1 (MQ=255)                                                            tcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc                                                                                     >  2:963978/1‑151 (MQ=255)                                                                  ccTGATGCGACGCGGGCGCGTTTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCCCCTTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                                <  1:345796/150‑1 (MQ=255)                                                                    tGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                                             >  1:893082/1‑151 (MQ=255)                                                                        gCGACGCGGGGGGTTTTTATCATGCCTACATATTTTTCATATTTTACACCCGGCAC‑CCCCCGTTTCCCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg                                                                          <  2:793154/150‑1 (MQ=255)                                                                          gACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                                       >  1:1128184/1‑151 (MQ=255)                                                                                ggCGCTTTTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAA‑CCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGc                                                                    <  1:541485/148‑1 (MQ=255)                                                                                ggCGCGTTTTATCATCCCTACTTTTTTTTAATTTTTTACACCCGCCAC‑CCCCCTTTTACCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac                                                                  <  2:730483/150‑1 (MQ=255)                                                                                  cgctttttATCATGCCTACATATTTTTAATTTTTTCCATCCGCCAA‑CCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacca                                                                <  1:582016/144‑1 (MQ=255)                                                                                         tATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                                        >  1:743163/1‑151 (MQ=255)                                                                                         tATCATGCCTACATATTTTTCATATTTTACACCCGGCAC‑CCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc                                                           <  1:862260/148‑1 (MQ=255)                                                                                         tATCATGCCTACATATTTTTCATATTTTACACCCGGCAC‑CCCCCGTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTc                                                           <  2:468504/148‑1 (MQ=255)                                                                                          aTCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                                       >  1:39453/1‑151 (MQ=255)                                                                                           tCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                                       >  1:8469/1‑150 (MQ=255)                                                                                           tCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                                       >  1:534636/1‑150 (MQ=255)                                                                                            cATGCCTACATTTTTTTCATATTTTCCACCCGGCAC‑CCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                                        <  1:918004/148‑1 (MQ=255)                                                                                                ccTACATTTTTTTCATTTTTTACATCCGGCAC‑CCACCTTTTACCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata                                                    <  2:656535/148‑1 (MQ=255)                                                                                                   aCATTTTTTTCATATTTTACACCCGACAA‑CCCCCGTTTCCCCCTTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                 <  1:364875/148‑1 (MQ=255)                                                                                                    cATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                              >  2:582377/1‑150 (MQ=255)                                                                                                        tttttCATTTTTTACCTCGGGCAA‑CCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                                                <  1:911214/144‑1 (MQ=255)                                                                                                        tttttCATATTTTCCATCCGGCAA‑CCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                              <  1:675793/146‑1 (MQ=255)                                                                                                           ttctttttttccccccGGCAC‑CCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCAATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCAt                                       <  2:188934/136‑1 (MQ=255)                                                                                                              atatTTTACATCCGCCAC‑CCCCCTTTCCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt                                    <  1:1092294/150‑1 (MQ=255)                                                                                                               tatTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                                   >  1:537697/1‑150 (MQ=255)                                                                                                                 ttttACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                                 >  2:881282/1‑150 (MQ=255)                                                                                                                 tttaccacccGCCAC‑CCACCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                                   <  2:11179/140‑1 (MQ=255)                                                                                                                  tttACCCCCGGCAC‑CCCCCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                                <  2:288781/150‑1 (MQ=255)                                                                                                                    tACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                              >  2:174722/1‑150 (MQ=255)                                                                                                                    tACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                             >  1:730805/1‑151 (MQ=255)                                                                                                                    tACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                             >  2:278265/1‑151 (MQ=255)                                                                                                                         ccGGCAC‑CCCCCTTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc                        <  1:915026/151‑1 (MQ=255)                                                                                                                         ccGCCAC‑CCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg                         <  1:327015/150‑1 (MQ=255)                                                                                                                          cGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg                       >  1:539700/1‑151 (MQ=255)                                                                                                                              caccccccttttCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg                       <  2:741004/139‑1 (MQ=255)                                                                                                                               accccccttttCCCCCATCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc                    <  2:951021/142‑1 (MQ=255)                                                                                                                                  ccccGTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata                <  2:202919/149‑1 (MQ=255)                                                                                                                                           ccccGACCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGTGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCaa         <  2:76262/149‑1 (MQ=255)                                                                                                                                              cGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATcc    <  2:502055/151‑1 (MQ=255)                                                                                                                                                 caccacCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCgg  >  1:221516/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGG  >  NC_000913/501428‑501720

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