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

GGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCA  >  NC_000913/501419‑501712                                                                                                                                                  |                                                                                                                                                     gggATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCg                                                                                                                                                 <  2:897473/151‑1 (MQ=255)    aTCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTca                                                                                                                                              <  1:28492/151‑1 (MQ=255)        ttAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcaccac                                                                                                                                          >  1:275609/1‑151 (MQ=255)           aGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTc                                                                                                                                       >  2:216846/1‑151 (MQ=255)            gTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCa                                                                                                                                      <  2:728815/151‑1 (MQ=255)                tcgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccg                                                                                                                                  >  1:584188/1‑151 (MQ=255)                  gtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgcc                                                                                                                                <  1:835826/151‑1 (MQ=255)                  gtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgcc                                                                                                                                >  1:649868/1‑151 (MQ=255)                  gtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccg                                                                                                                                  <  2:66358/149‑1 (MQ=255)                   tcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgcc                                                                                                                                >  2:19176/1‑150 (MQ=255)                     ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCTTTTGGCTCTCGATGCTTGTTGCGACGCGGGCGCTTCTTTTCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTCCCCCGCCCCCACCCCACccgccg                                                                                                                               <  2:137929/149‑1 (MQ=37)                        tCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGc                                                                                                                          >  2:449998/1‑151 (MQ=255)                            acgCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGAAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGttt                                                                                                                      >  2:371361/3‑151 (MQ=255)                             ggCCAGAAAGGGCGTCCACGCCGAATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTTTTATCATGCCTACATATTTTTCATATTTTCCACCCGGCACCCACCGTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCg                                                                                                                         <  2:673204/147‑1 (MQ=37)                                     caGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGtt                                                                                                             >  2:314112/2‑151 (MQ=255)                                     agggCTTTCCCGCCGCATCTGGCTTTTGCCTCTCGATGCCTGATGCGCCGCTGGCGCTTCTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCAACCCCCTTTTCCCCCGTCCCCCCTTCACCCGCCGGTGGCGTTTCCAGCAGtt                                                                                                             <  1:53424/149‑1 (MQ=255)                                                  cgcATCTGGTTTTTGGCTCTCGATGCTTGATGCGACGCGGGCGCGTTTTATCTTGCCTACTTTTTTTTCATTTTTTACACCCGGCAACCACCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACg                                                                                                <  2:673173/151‑1 (MQ=255)                                                      tCTGGCATTTGGCTCTCTATGCCTGATGCGCCGCGGGCGCTTTTTATCATCCCTCCATATTTTTCATATTTTACACCCGGCACCCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc                                                                                             <  2:63733/150‑1 (MQ=255)                                                       cTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCg                                                                                            >  2:920608/1‑150 (MQ=255)                                                         ggCTTTTGCCTCTCGTTGCCTGAGGCGCCGCTGGCGCGTCTTATCATGCCTACATTTTTTTAATTTTTTCCATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc                                                                                             <  2:21757/147‑1 (MQ=255)                                                                     tcGATGCCTTATGCGGCGCTGGCGCGTCTTCTCATGCCAACATATTTTTAATATTTTACCTCCGGCAACCACCGTTTCCCCCGTCCCCCCCCCAGCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCAc                                                                                     <  2:270574/143‑1 (MQ=255)                                                                     acGATGCTTGTTGCGACGCTGGCGCGCTTTATCATGCATACATATTTTCCATTTTTTCCACCGGACACCCGCCTTTTCCCCCGCCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc                                                                             <  1:385286/150‑1 (MQ=255)                                                                       gATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                                                                           >  2:484872/1‑151 (MQ=255)                                                                       gATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                                                                           >  2:159493/1‑151 (MQ=255)                                                                             tGATGCGACGCTGGCGCTTCTTATCATCCTTACATATTTTCCATTTTTTACACCCGGCAACCACCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                       <  1:43785/149‑1 (MQ=255)                                                                             tGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                                     >  1:67191/1‑151 (MQ=255)                                                                                tGCGACGCTGGCGCGGCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCt                                                                   >  2:294750/1‑150 (MQ=255)                                                                                 gCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg                                                                  >  1:924426/1‑150 (MQ=255)                                                                                  cGACGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGa                                                                 <  1:784753/150‑1 (MQ=255)                                                                                        tGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac                                                          >  2:928231/1‑151 (MQ=255)                                                                                        tGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcac                                                          >  2:17619/1‑151 (MQ=255)                                                                                         ggCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                              <  1:947682/146‑1 (MQ=255)                                                                                          gcgcGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc                                                         >  2:664737/1‑150 (MQ=255)                                                                                           cgcgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcaccac                                                       >  2:247912/1‑151 (MQ=255)                                                                                            gcgTTTTATCATGCCAACATATTTTTCATATTTTACACCCGGAAACCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcaccac                                                       <  1:458494/150‑1 (MQ=255)                                                                                               tCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCCACACCACCTCACCCGCCGGTGGCGTTTCCAGCAGCTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                    >  2:781900/1‑150 (MQ=255)                                                                                                cTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGCTTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                  >  2:616779/1‑151 (MQ=255)                                                                                                         ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                         >  1:272719/1‑151 (MQ=255)                                                                                                               tatTTTTCATATTTTACACCCGGCAACCCCCGTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata                                            <  1:8705/142‑1 (MQ=255)                                                                                                                atTTTTCATATTTTACATCCGGCAACCACCGTTTACCCGGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                                  >  1:209504/1‑151 (MQ=255)                                                                                                                  tttcaattttttCCACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTc                                 <  1:163396/142‑1 (MQ=255)                                                                                                                   tttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCAt                               >  2:152278/1‑151 (MQ=255)                                                                                                                     tCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGAAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCCGTGATATACGCCACTTCATCg                             >  1:937437/1‑151 (MQ=255)                                                                                                                        tatTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                           >  1:202120/1‑150 (MQ=255)                                                                                                                        tatTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                           >  2:637872/1‑150 (MQ=255)                                                                                                                           tttACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                        >  2:434696/1‑150 (MQ=255)                                                                                                                              acacCCGGCAACCACCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCATTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                      <  2:590289/149‑1 (MQ=255)                                                                                                                                  cgggCAACCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                    <  1:463519/145‑1 (MQ=255)                                                                                                                                  ccGGCACCCCCCTTTTCCCCCGCCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgc                  <  2:649491/149‑1 (MQ=255)                                                                                                                                  ccGCCACCCCCCGTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgc                  <  2:545917/149‑1 (MQ=255)                                                                                                                                   cGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCGTAATGGGcgcg               >  2:529870/1‑151 (MQ=255)                                                                                                                                       acccaccGTTTCCCCCGCCCCCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCa           <  1:827821/149‑1 (MQ=255)                                                                                                                                       aaCCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc            >  2:78607/1‑150 (MQ=255)                                                                                                                                         cccccGTTTCCCCCTCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCggg             <  2:277959/144‑1 (MQ=255)                                                                                                                                          ccccGTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata        <  1:244135/149‑1 (MQ=255)                                                                                                                                          ccccGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata        <  1:304330/149‑1 (MQ=255)                                                                                                                                            ccttttCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat      <  2:586884/148‑1 (MQ=255)                                                                                                                                            ccGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat      <  1:742198/151‑1 (MQ=255)                                                                                                                                             cGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc     >  1:694583/1‑151 (MQ=255)                                                                                                                                                 tACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCa  >  1:587981/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCA  >  NC_000913/501419‑501712

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