Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,628,232	Δ1 bp	coding (918/1233 nt)	nadR →	trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,628,232	0	A	.	98.3%	272.3 / ‑1.3	60	coding (918/1233 nt)	nadR	trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (29/29);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGG  >  NC_000913/4628088‑4628374                                                                                                                                                 |                                                                                                                                               caGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCg                                                                                                                                         <  1:514992/151‑1 (MQ=255)  aGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGa                                                                                                                                        <  2:511581/151‑1 (MQ=255)       gTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCg                                                                                                                                   >  1:785581/1‑151 (MQ=255)         gggAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGa                                                                                                                                  <  2:88327/150‑1 (MQ=255)         gggAATATGGCCGCGATTATGTCTTTTCACACCTCCGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGa                                                                                                                                  <  2:470953/150‑1 (MQ=255)           gAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGAttt                                                                                                                               <  1:548760/151‑1 (MQ=255)            aaTATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGAtttt                                                                                                                              >  2:544418/1‑151 (MQ=255)                 ggCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTcac                                                                                                                         >  1:1311653/1‑151 (MQ=255)                    cgcgATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTcaccac                                                                                                                      <  2:1248850/151‑1 (MQ=255)                     gcgATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGGGGC‑TTTATCGATACCGATTTTGTcacc                                                                                                                        >  2:126662/1‑148 (MQ=255)                      cgATTATGTCTTTTCACACCTCGGCGGCGTTGAGATCGCATTGCAGTATTCTGTCTTAGCTAAAATCGTGCTGGGTCCCGCTCAATACATTTCTTTTGCGGTGAAAAATGCCAATTAAGGGGC‑TTTTTCCAGACCGATTTTGTca                                                                                                                          >  2:1009415/1‑145 (MQ=255)                       gATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCCATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTc                                                                                                                    >  2:1223649/1‑150 (MQ=255)                          tATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGc                                                                                                                >  1:419357/1‑151 (MQ=255)                           aTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGc                                                                                                                <  2:1037521/150‑1 (MQ=255)                            tGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGt                                                                                                              >  1:835394/1‑151 (MQ=255)                                 tttGACACCTCGCCGGGGATGAGATCGCATTTCAGTATTCCGGCTACGACAAAATCGCGCCGGGCCCCGCTCAAAACATTGATTTTGCGGCAAAATAAGCCAATAAAGTGGCCTTT‑TCGATACCGATTTTGTCACCAATCAGGCGTTCt                                                                                                           >  2:92337/1‑149 (MQ=255)                                  ttCACACCTCGGCGGTGATGGGATCGCATTGCCGGATTCTGACCACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATAAGCCAATAAAGTGGC‑TTTAACGATACCGATTTTGTCACCACTCAGGCGTTCTg                                                                                                          >  2:742193/1‑149 (MQ=255)                                    cacaCCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCaaa                                                                                                      <  2:129573/151‑1 (MQ=255)                                           cggcggTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTa                                                                                                 <  1:638695/149‑1 (MQ=255)                                            ggcggTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACg                                                                                               >  1:937555/1‑150 (MQ=255)                                              cggTGATGAGGTCGCATTGCAGTATGCTGACTACGGTAAAATCGCGCTGTGTGTCGCTCAATACATTGATTTTGCGGTGAAATATGCGAAGAAAGTTGC‑TTTATCGATACCGATTTTGTCACCACGCAGGTGTTCTGCaaaaactgctaa                                                                                             <  2:847668/150‑8 (MQ=255)                                               ggTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAgg                                                                                           <  2:1299618/151‑1 (MQ=255)                                                  gatgaGATCGCATTGCAGTATTCTGACTCCGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGACGGGCg                                                                                        <  2:1173636/151‑1 (MQ=255)                                                   atgaGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCgg                                                                                       >  1:1137027/1‑151 (MQ=255)                                                   atgaGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCgg                                                                                       >  1:542292/1‑151 (MQ=255)                                                            cATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATAGGCCAATAAAGTGGC‑TTTATAGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGAACGAAGAGCGGGAACATCCg                                                                              >  2:33660/1‑151 (MQ=255)                                                               tGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTc                                                                           >  1:252032/1‑151 (MQ=255)                                                               tGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTc                                                                           >  1:510713/1‑151 (MQ=255)                                                               tGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTc                                                                           >  2:912728/1‑151 (MQ=255)                                                               tGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGGCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGtt                                                                            <  2:1198120/150‑1 (MQ=255)                                                                gCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTc                                                                           <  1:1430322/150‑1 (MQ=255)                                                                gCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCg                                                                          >  1:988710/1‑151 (MQ=255)                                                                     aTTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAg                                                                     <  2:1231130/151‑1 (MQ=255)                                                                      ttCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAgg                                                                    >  2:394804/1‑151 (MQ=255)                                                                              aCGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtgatt                                                            <  1:1262823/151‑1 (MQ=255)                                                                              aCGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtgat                                                             >  1:1409630/1‑150 (MQ=255)                                                                               cGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtgattg                                                           >  1:308706/1‑151 (MQ=255)                                                                                  tAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATtgatg                                                        >  1:1072637/1‑151 (MQ=255)                                                                                        cgcgCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATAc                                                   <  1:986248/150‑1 (MQ=255)                                                                                         gcgcTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCg                                                 >  1:796892/1‑151 (MQ=255)                                                                                              gggCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCg                                            <  1:1093663/151‑1 (MQ=255)                                                                                                 ccACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATc                                         <  1:1025631/151‑1 (MQ=255)                                                                                                       tCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTg                                    >  1:1162218/1‑150 (MQ=255)                                                                                                        cAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGa                                   <  1:544595/150‑1 (MQ=255)                                                                                                                tGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCctgctg                           <  2:785355/150‑1 (MQ=255)                                                                                                                tGATTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCctgctg                           <  2:1159937/150‑1 (MQ=255)                                                                                                                  aTTTTGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTgg                          >  2:1251474/1‑149 (MQ=255)                                                                                                                    tttGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTGGTCCCCCCTCAGGCGTTCCGCAAAAAGTACGAAGGGCGGGAACCTCCGTCCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCctgctg                           <  1:419463/146‑1 (MQ=255)                                                                                                                      tGCGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGaac                     <  2:5942/150‑1 (MQ=255)                                                                                                                        cGGTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGaacaac                  >  2:196612/1‑151 (MQ=255)                                                                                                                          gTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAacac                <  2:510541/151‑1 (MQ=255)                                                                                                                          gTGAAATATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAacac                <  1:534196/151‑1 (MQ=255)                                                                                                                             aaaTATGCCAATAAAGAGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGGGATCCTGCTGGAGAACAACACGc              >  1:1089604/1‑150 (MQ=255)                                                                                                                             aaaTATGACAATAAAGTGG‑ATATATCGATACAGATTTTGTCACCACTCAGGCGATTTGCAAACAGTACGAAGGGCGGAAACATCCGATAGTGCAGCCGCTGATTGATGAATACCGTTTCGATCCGGTAATCCTGCTGGAGAACAACACGc              >  2:42977/1‑150 (MQ=255)                                                                                                                              aaTATGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGcc             >  2:244016/1‑150 (MQ=255)                                                                                                                               atatGCCAATAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTATGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCg            >  2:1143058/1‑150 (MQ=255)                                                                                                                                   gCCAATAAAGTTGC‑TTTACCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGTGCGGAAACATCCGTTCGTGCAGGCGCTGATTGATGAATGCCGTTTCGATCTGGTGATCCTGCTGGAGCACACCACGCCgtgagt       <  2:639374/151‑1 (MQ=255)                                                                                                                                      aaTAAAATGGC‑TTTATCGATACCGATTTTGTCACCACACAGGCGTTCTGCAAAGAGTACGAAGGGCTGGAGCAGCGGTTCGTGCAGGCTCTGATTGATGAATGCCGTTTCGAATTGGTGGTCCTGGTTGGGAACGACCCGACttgggtgg     <  2:210582/150‑1 (MQ=255)                                                                                                                                       aTAAAGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCg   >  2:313262/1‑151 (MQ=255)                                                                                                                                         aaaGTGGC‑TTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCgg  <  2:571197/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                               CAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGG  >  NC_000913/4628088‑4628374

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