Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,820,036	C→T	D379N (GAT→AAT)	ligB ←	DNA ligase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,820,036	0	C	T	94.6%	111.8 / ‑3.4	37	D379N (GAT→AAT)	ligB	DNA ligase B
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (22/13);  total (23/14)
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

TAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGC  >  NC_000913/3819794‑3820282                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       tAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGCTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAGCCGGCCTCACCAATGCCATCCAGCCCAAGACCCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGc                                                                                                                                                                                                                                                >  1:248079/1‑251 (MQ=255)  aaTTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGCTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAGCCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCg                                                                                                                                                                                                                                               <  2:373594/251‑1 (MQ=255)           ggCGTTAATAAAAGCCAGGAAAAGATATGCTCAAAGCGATGAGTATGATGCAGCGCGCGCCAGCCGGCATCACCAATGCCATCCAGACCAAGAAACTGTTTTGCCCCCAGCCCGACTAAGCGGGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGAGACGGTTTTGTACGTTCTGCACCGCGCCACACCAAATTATCAATGCGAGGAATAcc                                                                                                                                                                                                                                      >  1:419968/1‑251 (MQ=255)           ggCGTTAATAAAAGCCAGGAAAAGATATGCACAAAGCGATGAGTCTGATGCAGCGCGCGCCAGCCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATa                                                                                                                                                                                                                                        <  2:139950/249‑1 (MQ=255)                      aaGCCAGGAAAAGATATGCTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAGCCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCa                                                                                                                                                                                                                            >  1:310807/1‑250 (MQ=255)                          cAGGAAAAGATATGCTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAGCCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCt                                                                                                                                                                                                                       <  1:26090/251‑1 (MQ=255)                             gAAAAGATATGCTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAGCCGGCATCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCCCACCACATTATCAATGCGAGGAAAACCCTGACCGGCTAGGCTggc                                                                                                                                                                                                                    >  2:128912/1‑249 (MQ=255)                                         tCAAAGCGATGAGTCTGATGCAGCGCGCGCCAGCCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGACCTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGCGAATCTGATc                                                                                                                                                                                                        >  1:432836/1‑251 (MQ=255)                                           aaaGCGATGAGTCTGATGCAGCGCGCGCCAGCCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCAc                                                                                                                                                                                                      >  2:399682/1‑251 (MQ=255)                                             aGCGATGAGTCTGATGCAGCGCGCGCCAGCCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCAcc                                                                                                                                                                                                     >  1:209181/1‑250 (MQ=255)                                                           tGCAGCGCGCGCCAGCCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGGTTTGTACGTTCTGCACCGCGCCACACCACCTTATCAATGCGAGGAATACCCTGCCCGGCAAGGCTGACGAGAATCTGATCACCCGGCGCAATATcccc                                                                                                                                                                                      >  2:123927/1‑250 (MQ=255)                                                             cAGCGCGCGCCAGCCGGCCGCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACt                                                                                                                                                                                    <  1:49494/251‑1 (MQ=255)                                                                    cgcCAGCCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGcc                                                                                                                                                                              >  2:348286/1‑250 (MQ=255)                                                                         gccGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGcc                                                                                                                                                                        <  2:310807/250‑1 (MQ=255)                                                                                    cAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTTGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCa                                                                                                                                                             <  1:31602/251‑1 (MQ=255)                                                                                        gCCAACCAGCCCAAGAACCTGTCTTGCCCACAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCCGCGCCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAAGATCCCACTCCTGCCAGCGCCTGACGGAACCAatat                                                                                                                                                         <  2:98106/251‑1 (MQ=255)                                                                                              cAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCAccc                                                                                                                                                   >  1:327961/1‑251 (MQ=255)                                                                                                        aCCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGCGTTCAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACAACCCATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATAACCCAATCCTGCCAGCGACTGACGGAACCAATATTCACCCGCTGGCCttt                                                                                                                                         >  1:428372/1‑251 (MQ=255)                                                                                                                 gCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCCAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTatcatc                                                                                                                                >  2:73586/1‑251 (MQ=255)                                                                                                                      cagccagACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAAGCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCa                                                                                                                           <  1:430080/251‑1 (MQ=255)                                                                                                                                      gAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGa                                                                                                           <  1:271328/251‑1 (MQ=255)                                                                                                                                          tGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTATGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCCATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGCTCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCCCCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGc                                                                                                        >  2:425041/1‑250 (MQ=255)                                                                                                                                            aaCTGTTCCTGACAAAAATCAGAAGCAAAGTAGCCGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATCTTCACCCGCTGGCCTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAac                                                                                                     >  2:83113/1‑251 (MQ=255)                                                                                                                                             aCTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAAGGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCCCCCGCTGGACTTTTTTATCATCCAGCAGGACAGGTGCGAGTGACGCAacc                                                                                                    <  2:294002/251‑1 (MQ=255)                                                                                                                                                   ccTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGat                                                                                              >  1:120145/1‑251 (MQ=255)                                                                                                                                                      gACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCGGCGCCt                                                                                                                                                                       >  1:106937/1‑175 (MQ=255)                                                                                                                                                      gACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCt                                                                                                                                                                       <  2:106937/175‑1 (MQ=255)                                                                                                                                                               cAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGct                                                                                   <  1:348286/250‑1 (MQ=255)                                                                                                                                                               cAGAAGCAAAGGAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGctc                                                                                  >  1:114908/1‑251 (MQ=255)                                                                                                                                                                        aGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGCCTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCAcc                                                                         >  2:228991/1‑251 (MQ=255)                                                                                                                                                                                      gAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCCAGGCTGACGCGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGCACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAAtcgc                                                           >  1:428522/1‑248 (MQ=255)                                                                                                                                                                                       aGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAAAGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCAAGCGCCTGACGAAACCAATATTCACCCGCTGGCCTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCACCCACCGCTATTTTACAGCTCTTACCCACCTCAAACTGAATCGcc                                                          >  1:284662/1‑251 (MQ=255)                                                                                                                                                                                        gTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATAACCTGACCGGCAAGGCTGACGAGAATCTGATCCACAGGCGCAATATCCCACTCCTGCCAGCGCCTGAAGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAAGTGCGAGTGCCGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATCGCCt                                                         >  1:270313/1‑251 (MQ=255)                                                                                                                                                                                              gCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATCGCCTTAACtt                                                   >  1:40964/1‑251 (MQ=255)                                                                                                                                                                                                                       gTACGTTCTGCACCGCGCCACACCACATTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATCGCCTTCACTTCGGCAACCTGAGCTACAGGTTGat                           <  2:364422/250‑1 (MQ=255)                                                                                                                                                                                                                                gCACCGCGCCACACCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGAAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAgg                  <  2:120145/250‑1 (MQ=255)                                                                                                                                                                                                                                               acaTTATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATCGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGc  >  1:403223/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       TAATTGCTCTGGCGTTAATAAAAGCCAGGAAAAGATATGTTCAAAGCGATGAGTCTGATGCAGCGCGCGCCAACCGGCCTCACCAATGCCATCCAGCCCAAGAACCTGTTTTGCCCCCAGCCAGACTAAGCGTGAAATGAACTGTTCCTGACAAACATCAGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGC  >  NC_000913/3819794‑3820282

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