Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,157,116	Δ1 bp	coding (267/1374 nt)	argH →	argininosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,157,116	0	A	.	81.8%	73.5 / 22.6	34	coding (267/1374 nt)	argH	argininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base . (7/20);  total (11/23)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

TACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGAAGGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAACCGCACAAAACAATCAGGACGCGGTAATGCCAGGTTACACTCACCTGCAACGC‑GCCCAGCCGGTGACGTTCGCGCACTGGTGCCTGGCCTATG  >  NC_000913/4156873‑4157366                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            ttcccAGGCAGCACATCAACGGTTCAAAAAATGCAACGACTCAATGCGATTTGATTCCCGTGTGGCGGAGTAGGATATTGTTGGTGCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGTGCTGAGCGTGTTGCTGGAAGATGTTCGCGCCAGGCCACCACAAATCCTTGAAATCGATGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAAc                                                                                                                                                                                                                                                      <  2:149090/248‑1 (MQ=255)         cagcagATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGAc                                                                                                                                                                                                                                              >  1:159936/1‑250 (MQ=255)                             aaaTTCAACGACTCCATGCGCTTTGATGTCCGTCTGGCTGAGCAGGATATTGTTGGCTCTGTGTCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAAGCCTTGAAAGCGACGTCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCaaa                                                                                                                                                                                                                         <  2:237174/250‑1 (MQ=255)                               aTTCAACGTCTCACTGCTCTTTGATTACCTTAGGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAATAGCAGGCGCAACTTGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAATCGACGCCGAAGATATCCATAGCTTGTTGGA‑GGTAAACAGATCGACAAAGTGCGCCAGTTAGGCaaaaa                                                                                                                                                                                                                       <  2:34900/251‑1 (MQ=255)                                        cTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCATGGTAAAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGCCGCCGAAGATATCCATAGCTGGGTGGC‑GGCAAACTGATAGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATAc                                                                                                                                                                                                              >  1:215773/1‑251 (MQ=255)                                              gcgcTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGc                                                                                                                                                                                                         <  1:152841/250‑1 (MQ=255)                                                             tCTGGCGGAGCAGGATTTTGGTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAAAAGAGCAGGCGCACCTGGAAGAGACGCTGAACGTGTTGCTGGCAGATGTTCGCGCCAGGCCACAACACATCCTTGAAAGCGACGCCGGACAGATCCCAAGCTGGGTGGA‑GGACATCTGTTCGTCCAAGTGGGCCAGTTAGGCCAAAAGCTGCACACAGAGCGTCGCCGGAATGagtc                                                                                                                                                                                         >  2:76696/1‑248 (MQ=255)                                                                ggctgAGCAGGATATTGTTGGCTCTGTGGCCTGGTCAAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACGGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGCAGCCGTAATGATCAg                                                                                                                                                                                        <  1:244415/245‑1 (MQ=255)                                                                       aaGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCTCGGTAGGCGTGTTAACCGAAGAAGAGCAGGCGCAACTGGAAGAGGAGCTGACCGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCTATGAAAGCGACGCCGAAGATATCCATAGCTGGTTGGA‑GGCAAACTGATCGACAAAGTGGTCCAGTTAGGCAAAAAGCTGCATACCGGGAGTAGCCGTAATGATCAGGTAGCGACt                                                                                                                                                                               <  2:72301/250‑1 (MQ=255)                                                                                ggtgGCTCTGGGGGCTGGGCCAAAGGCCTGTGCACGGGAGGCGTGGTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCGGGAAGATTGTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGTGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGaa                                                                                                                                                                       <  2:38558/248‑1 (MQ=255)                                                                                       ctGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACtgtg                                                                                                                                                                 <  1:147356/249‑1 (MQ=255)                                                                                         gtgGCCTTGTACCAAGCCCTGGGCACGGTCGGCGTGTTAACCGCTGTAGAGCAGGCGCAACTGGAAGAGGCGATGAACGGGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCAGTGTAAGCGATGTCGAAGATATCCATAGCTGGG‑GGGAGGCAAACTGGCCTACAAAGTGGGACAGTTAGGCGAAAAGCTGCATACCGGGCGTATCCGTAAGGATCAGGTAGCGACTGACCTGAAACTGTGGTGc                                                                                                                                                             <  2:128004/251‑1 (MQ=255)                                                                                         gtgGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTgtgg                                                                                                                                                                <  1:81508/248‑1 (MQ=255)                                                                                                   ccAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATAcc                                                                                                                                                    <  1:235671/250‑1 (MQ=255)                                                                                                                         cGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTa                                                                                                                              >  2:48462/1‑250 (MQ=255)                                                                                                                                aCCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGAAGGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCa                                                                                                                       <  1:155373/251‑1 (MQ=255)                                                                                                                                        agagCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATc                                                                                                              <  1:307959/251‑1 (MQ=255)                                                                                                                                         gagCATGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGACGATGTTCGCGTCAGGCCAAAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGG‑GGGAGGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCg                                                                                                             <  2:74560/251‑1 (MQ=255)                                                                                                                                          agCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGAAGGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCg                                                                                                             >  2:31211/1‑251 (MQ=255)                                                                                                                                                       ggaagcgacgcTGGACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTAGAACGCGACGCCGAAGCTATCCCTAGCTGGGTG‑TCGGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAACAGCTGCATACCGGGAGAAGCCGTAATGATCAGGTAGCGACTGACATGAAAATGTGGTGCAAAGATACCGATAGCGAGTTACTGGCGGCTAACCGGCAACAGACATCGGCGCTGGTGGAAAc                                                                                               >  2:166434/9‑251 (MQ=255)                                                                                                                                                        gAAGAGGGTCTGAACGTTTTGCGGGAAGCTGTTCGCGCAAGGCCACACCAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGACCGACAAAGTTGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGCCCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAAcc                                                                                              <  2:277770/251‑1 (MQ=255)                                                                                                                                                                      cGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAACCGCACAAAACAATCa                                                                                <  1:48462/251‑1 (MQ=255)                                                                                                                                                                                    aTGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAACCGCACAAAACAATCAGGACGCGGTAATGc                                                                  >  1:3844/1‑251 (MQ=255)                                                                                                                                                                                      gtgtgtgacgggCCACAACAAATCCTTGAACGCGCCGTCGTAGATATCCATAGCTGGGTGGA‑GGCAAAGGGATAGACAAAGTGGGTCAGTTAGGCAAAAAGCTGCATACCGTGTGTAGCCGTAATGATCAGGTAGCGATTGACCTGAAACTGTGGTGCAAAGATGCCGTTTGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGATGGGTTAAACCGCACAAAATAATCAGGACGCTGTAATGCCa                                                                <  2:180881/241‑1 (MQ=255)                                                                                                                                                                                           cgcCAGGCCTCAACAAATACTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGACAGTTAGGCAAAAAGCTGCATAACGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAACCGCACAAAACAATCAGGACGCGGTAATGCCAGGTTa                                                           <  1:89533/251‑1 (MQ=255)                                                                                                                                                                                           cgcCAGGCCACAACAAATCCTTGAAAGCGACGCCGACGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAACCGCACAAAACAATCAGGACGCGGTAATGCCAGGtt                                                            >  1:282670/1‑250 (MQ=255)                                                                                                                                                                                           cgcCAGGACACAAAAAAGCCTTGAAAGCGACGACGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGGAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAACCGCACAAAACAATCAGGTCGCGGTAATGCCAGGtt                                                            <  2:339474/250‑1 (MQ=255)                                                                                                                                                                                              cAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGAAGGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGAATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGCCCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCTGCAGCTGCAATCGGCGCTGGTGGAAACCGCACAAAACAATCAGGACGCGGTAATGCCAGGTTac                                                          >  1:50895/1‑250 (MQ=255)                                                                                                                                                                                              cAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGA‑GTCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAGCCGCACAAAACAAGCAGGACGCGGTAATGCCAGGTTacac                                                        <  1:92309/251‑1 (MQ=255)                                                                                                                                                                                                        caaATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGCGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGCGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAACGGACCAAAACAATCAGGACGAGGGAATGCCGTGTTACACTCACCTGCac                                              >  2:173678/1‑250 (MQ=255)                                                                                                                                                                                                                       gACGCCGAAGATATCCATAGCTGGGTGGA‑GGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATTCCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAGCCGGCAGCTGCAATCGGCGCTGGTGGAAACCGCACAAAACAATCAGGTCGCGGTAATGCCAGGTTACACTCACCTGCAACGC‑GCCCAGCCGGTg                              <  2:161679/251‑1 (MQ=255)                                                                                                                                                                                                                        aagcagaagattaccATAGCTGGGTGGG‑GGCAAACTGATCGACAAAGTGGGCCAGTGAGGCAAAAAGCTGCATCCCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAACCGCACAAAACAATCAGGACGCGGTAATGCCAGGTTACACTCACCTGCAACGC‑GCCCAGCCGGTGa                             <  2:83308/238‑1 (MQ=255)                                                                                                                                                                                                                                       aTAGCTGGTTGGA‑GGCAAACTGATCGACTAATTGGGCCAGTTAGGCAAAAAGCTGCATGCCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGTCGCTGGTGTAAACCGCACAAAACAATCAGGACGCGGTAATGCCAGGTTACACTCACCTGCAACGC‑TCCCAGCCGGTGACGTTCGAGCACTGGt              <  1:259431/251‑1 (MQ=255)                                                                                                                                                                                                                                                 ggAAGGCAAACTGATCGACCCAGTGGGCCAGTTAGGAAAAAAGCTGCATACCCGCCGTAGCAGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATCCCGTTAGCGCGTTACTGACGGCTAACCGACAGATGCAATCGGCGCAGGTGGAAACCGCACAAAACACCCCGGCCGCGGTAATGCCAAGTTACACTCACCTGCAACGC‑GCCCAGCCGGTGACGTTCCAGCACTGGAGACTGGCCt     >  2:166290/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gaGGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAACCGCACAAAACAATCAGGACGCGGTAATGCCAGGTTACACTCACCTGCAACGC‑GCCCAGCCGGTGACGTTCGCGCACTGGTGCCTGGCCTAt   <  1:30825/250‑1 (MQ=255)                                                                                                                                                                                                                                                    aGGCCAACTGATCGACAACGTGGGCCAGTTCGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGCGCAGGTAGCGACTGAACTGAAACTGTGGCGCAAAAATACCGTTAACGAGTTACTAACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAACCGCACAAAATAATCAAGACGCGGTAATGCCCGGTTCCACCCACCTGCAACGCCACCCA‑CCGGCGACGGTCACGCACTGGTGCCTGGCCTATg  >  2:259808/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            TACCCAGGCAGCAGATCAACGGTTCAAACAATTCAACGACTCACTGCGCTTTGATTACCGTCTGGCGGAGCAGGATATTGTTGGCTCTGTGGCCTGGTCCAAAGCCCTGGTCACGGTAGGCGTGTTAACCGCAGAAGAGCAGGCGCAACTGGAAGAGGCGCTGAACGTGTTGCTGGAAGATGTTCGCGCCAGGCCACAACAAATCCTTGAAAGCGACGCCGAAGATATCCATAGCTGGGTGGAAGGCAAACTGATCGACAAAGTGGGCCAGTTAGGCAAAAAGCTGCATACCGGGCGTAGCCGTAATGATCAGGTAGCGACTGACCTGAAACTGTGGTGCAAAGATACCGTTAGCGAGTTACTGACGGCTAACCGGCAGCTGCAATCGGCGCTGGTGGAAACCGCACAAAACAATCAGGACGCGGTAATGCCAGGTTACACTCACCTGCAACGC‑GCCCAGCCGGTGACGTTCGCGCACTGGTGCCTGGCCTATG  >  NC_000913/4156873‑4157366

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