Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,435,124	C→T	W481* (TGG→TAG)	cpdB ←	2',3'‑cyclic‑nucleotide 2'‑phosphodiesterase/3'‑nucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,435,124	0	C	T	97.4%	106.0 / NA	38	W481* (TGG→TAG)	cpdB	2',3'‑cyclic‑nucleotide 2'‑phosphodiesterase/3'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/18);  minor base G (0/1);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCCAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCACCGGC  >  NC_000913/4434879‑4435366                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     aaGCAAAAGCGATAGGGCTGGCGCCCGTACCGGCAAATTTGCCGCCGTTAGCGCTATAGTTATTGGTGGCAACCAGGAACATGGCGTGCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTtg                                                                                                                                                                                                                                                 <  1:317199/251‑1 (MQ=255)        agagATATGGCTGTCGCCCGTACCGTCAAATTCGCCTCCGTAAGCGCTATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGAGCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAAGCATACTAGTTGATGAGTg                                                                                                                                                                                                                                          <  1:164395/248‑1 (MQ=255)                                  aaaaTTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCa                                                                                                                                                                                                                <  1:558707/250‑1 (MQ=255)                                                    cgcgATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTGCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑ccc                                                                                                                                                                                             <  1:506258/251‑1 (MQ=255)                                                     gcgATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCcg                                                                                                                                                                                            >  1:521715/1‑251 (MQ=255)                                                                          cAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCtt                                                                                                                                                                      >  1:455851/1‑251 (MQ=255)                                                                           tggAACAAGGTGTTCGGATCACTAGGCTTGACATTAAAGGTCAGGCTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACCTCAATCTGATAATTCACACCATCAATATAATCCAAGTTATAAGTGCGGAAACCATACTAGTTGATGAGTGTCTGTGGTTTAGTGCGGTTAGGATCCATCTGGTTAAACTGT‑ACCGGCGAG‑ACGTCCAGCCACTCttt                                                                                                                                                                     >  2:79829/2‑251 (MQ=255)                                                                              aaCATGACGTTCGGGTCAATCGGCTTGCCATTAAAGGTCAGGATCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATCACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACACTCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGATTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCAAGCCAATCTTTCAc                                                                                                                                                                  >  2:268040/1‑251 (MQ=255)                                                                                       ttgggATCAATCGCCTTGCCATTAAAGGTCAGGTTTTTACTCCTTCCCTCATTGGCTTTAATCATCTGGCACTCGCCGTTATAACGGGCGGTTTGGGTACCATCAAGCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCTTGCTGTTGGGATCAATCTGGTTAAACTTC‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTAc                                                                                                                                                          <  2:449049/247‑1 (MQ=255)                                                                                                          cATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACAATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATc                                                                                                                                      >  1:533279/1‑251 (MQ=255)                                                                                                               aaGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGt                                                                                                                                 >  1:189242/1‑251 (MQ=255)                                                                                                                            tAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATag                                                                                                                       <  1:495740/248‑1 (MQ=255)                                                                                                                             attCCTTCCCGCATGGGCATTAATCATCGGGCACTCGCCGTCATACCGGTCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAATTTATAAGTGCGGAAACTAGTCTAGTGGATGATTGCCGTTGGTTTCGTGCGGTTGGTATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCGTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAgaga                                                                                                                    <  2:310855/248‑1 (MQ=255)                                                                                                                              aTCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGAt                                                                                                                   >  1:193933/1‑250 (MQ=255)                                                                                                                                ccTTACCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGACTGGGTAACATCAATCTTATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGATTGACTGTGTTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCTC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATaaa                                                                                                                <  1:7257/251‑1 (MQ=255)                                                                                                                                      cGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATaa                                                                                                                 >  1:350356/1‑244 (MQ=255)                                                                                                                                                          ggCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCa                                                                                        <  1:76926/249‑1 (MQ=255)                                                                                                                                                                 gCCGTCAAAACGGGTGGGCTGGGTAACATCAACCTGATAATTCACCCCATCAATAGCATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGATTGACGGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTTGCCtt                                                                               <  2:567013/251‑1 (MQ=255)                                                                                                                                                                         aaCGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCAcc                                                                       <  1:383159/251‑1 (MQ=255)                                                                                                                                                                           cgggcgggCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGACGGACAACAGGCCTTTTTCCACCTc                                                                     >  2:140015/1‑251 (MQ=255)                                                                                                                                                                              gcgggcTGGGTAACATCAATCTGATAATTCACACCATCAAAAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCcac                                                                  <  1:99671/251‑1 (MQ=255)                                                                                                                                                                                      ggTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTg                                                          >  1:176888/1‑251 (MQ=255)                                                                                                                                                                                       gTAAAATCAATTTGATAAGTCACACCATCAATAAAATCAAAGTTATAAGTGCGGAAACCATCCTATTTGATGAGTGACTGTGGTTTCGTGCTGTTGGTATTAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGc                                                         <  2:350356/251‑1 (MQ=255)                                                                                                                                                                                              cAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGTATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGcc                                                                                                                                                                            <  2:261198/129‑1 (MQ=255)                                                                                                                                                                                              cAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGcc                                                                                                                                                                            >  1:261198/1‑129 (MQ=255)                                                                                                                                                                                                         ttCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCa                                                                                                                                                                               >  1:336098/1‑115 (MQ=255)                                                                                                                                                                                                         ttCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACGGTGGTTTCTTGCTGTTGGGATCAATCTGGTTAAACGGT‑CCCGCGGAGCAC‑TCCa                                                                                                                                                                               <  2:336098/115‑1 (MQ=255)                                                                                                                                                                                                           cacaCCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCCCCTCCACCTAGATTGCCGGGTCATTTTTGCcaccac                                     >  2:211184/1‑251 (MQ=255)                                                                                                                                                                                                           cacaCCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGaccac                                     >  1:381239/1‑251 (MQ=255)                                                                                                                                                                                                             acccATCAAAAACTTAAAAGTTATAATGGCGGAAACCATCCGAGTTGATGGTTGACTGTGGTTTTGTGCTGTTGGGATCATTCTGTTTAAACTGT‑CCCGCGGACCAC‑TCCAGCCTCTCTTTCACCTCTTTACCGCTGGCTTTCACCGCAATCAGCGTATTGGGATAGATATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCCCATAGCTTGCCGGGTCATTTTTGCGacca                                      <  2:479293/246‑1 (MQ=255)                                                                                                                                                                                                              acCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCgc                                  >  1:298436/1‑250 (MQ=255)                                                                                                                                                                                                                 aTCAATAACATCAAAGTTATAAGAGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCCCCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACtt                               >  1:375065/1‑251 (MQ=255)                                                                                                                                                                                                                 aTCAATAACATCAAAGATATAAGTGCGGGAACCATCCTAGTTGATGAGAGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACATCCACATCGCTTGACGGGTCATTTTTGCGACCACCGTCtt                               >  2:426994/1‑251 (MQ=255)                                                                                                                                                                                                                       aaCATCAAAGTTATAAGTGCGGAAACCATCCTAGTTGATGGGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTTCGACCACCGACTTTAAACg                         >  1:202380/1‑251 (MQ=255)                                                                                                                                                                                                                             aaaGTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGcc                                                                                 <  1:546805/189‑1 (MQ=255)                                                                                                                                                                                                                                gTTATAAGTGCGGAAACCATCCTAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCCCCGACTTTACACGGTGCGGCac                >  2:242440/1‑250 (MQ=255)                                                                                                                                                                                                                                           ggAAACCATCCTGGTTGAGGAGTGCTTGTGGTTTCGTGTGGTTTGGATCAATCTGGTTAAAC‑GTGCCCGAGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCACCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCAcc     <  2:521715/251‑1 (MQ=255)                                                                                                                                                                                                                                              aaCCATCCTAGTTGATGAGTGACTGTGGTTTCGGGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCACCGGc  <  1:417330/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     AAGCAAAAGCGATATGGCTGTCGCCCGTACCGGCAAATTTGCCGCCGTAAGCGCGATAGTTATTGGTGGCAACCAGGAACATGGCGTTCGGATCAATCGGCTTGCCATTAAAGGTCAGGTTCTTAATCCTTTCCGCATTGGCATTAATCATCTGGCACTCGCCGTCATAACGGGCGGGCTGGGTAACATCAATCTGATAATTCACACCATCAATAACATCAAAGTTATAAGTGCGGAAACCATCCCAGTTGATGAGTGACTGTGGTTTCGTGCTGTTGGGATCAATCTGGTTAAACTGT‑CCCGCGGAGCAC‑TCCAGCCACTCTTTCACCTCTTTACCGCTGGCTTTCACCACAATCAGCGTATTGGGATAGAGATAAAGATCGGCGGCATTACGGAAGGTCAACTGGCCTTTTTCCACCTCCACATAGCTTGCCGGGTCATTTTTGCGACCACCGACTTTAAACGGTGCGGCAGCTGAAAGCACCGGC  >  NC_000913/4434879‑4435366

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