Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,367,038	G→A	I430I (ATC→ATT)	aspA ←	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,367,038	0	G	A	98.6%	239.1 / ‑3.2	71	coding (1290/1437 nt)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (38/32);  minor base . (0/1);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTA  >  NC_000913/4366792‑4367278                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 ttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCAGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAAttt                                                                                                                                                                                                                                              >  1:338902/1‑251 (MQ=255)         ttGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACg                                                                                                                                                                                                                                      >  2:519385/1‑251 (MQ=255)         ttGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACg                                                                                                                                                                                                                                      >  1:226010/1‑251 (MQ=255)             tatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGt                                                                                                                                                                                                                                  >  2:94685/1‑251 (MQ=255)                 aGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCAcc                                                                                                                                                                                                                              <  2:350416/251‑1 (MQ=255)                  gATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGACCGACTTC‑ACGTA‑CACTCTTACCGTTTT‑CGGCACAAATTTTACCCACGATGTCACCg                                                                                                                                                                                                                             >  2:610147/1‑251 (MQ=255)                  gATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCg                                                                                                                                                                                                                             <  1:373753/251‑1 (MQ=255)                   aTTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGt                                                                                                                                                                                                                            <  2:412888/251‑1 (MQ=255)                                    aCTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTCCCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAAc                                                                                                                                                                                                           <  1:489911/251‑1 (MQ=255)                                     cTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACg                                                                                                                                                                                                          >  2:135314/1‑251 (MQ=255)                                          cacaAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGtt                                                                                                                                                                                                      <  1:571518/250‑1 (MQ=255)                                          cacaAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGGTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGtt                                                                                                                                                                                                      <  2:435651/250‑1 (MQ=255)                                            caAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGTGTTc                                                                                                                                                                                                     >  1:645882/1‑249 (MQ=255)                                                 aaaaaaGCACGTAATCTGACGTGCCTTTTGTATTTGTACTACCCTGTAAGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCACCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGATCCAGAACGACTTC‑ACGAA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCAACGTTGTGGTGACGTATGAACGGGTGCAGgta                                                                                                                                                                                                >  2:620824/1‑249 (MQ=255)                                                  aaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGgtaagt                                                                                                                                                                                             <  1:527943/251‑1 (MQ=255)                                                    aaGGCACGTCATATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGgtaagta                                                                                                                                                                                            >  2:155354/1‑250 (MQ=255)                                                     aGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAAc                                                                                                                                                                                          >  1:536162/1‑251 (MQ=255)                                                     aGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAAc                                                                                                                                                                                          >  2:131689/1‑251 (MQ=255)                                                         aCGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGAt                                                                                                                                                                                       >  1:543692/1‑250 (MQ=255)                                                             cATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATAAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGGCAACAGACCGCGTTCCAGAACTACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTTCCCACGAGGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATCCCGa                                                                                                                                                                                  <  1:442285/251‑1 (MQ=255)                                                               tCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGAt                                                                                                                                                                                 >  1:368522/1‑250 (MQ=255)                                                                  gACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTGCCGGTTT‑CGGCACAAATTTTACCCACGTTGGCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATagag                                                                                                                                                                             <  1:373751/251‑1 (MQ=255)                                                                       gCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTa                                                                                                                                                                        <  1:238168/251‑1 (MQ=255)                                                                        ccTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTaa                                                                                                                                                                       <  1:284989/251‑1 (MQ=255)                                                                           tttttATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAAc                                                                                                                                                                     <  2:421800/250‑1 (MQ=255)                                                                             tttATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGt                                                                                                                                                                   >  2:166008/1‑250 (MQ=255)                                                                              ttATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTa                                                                                                                                                                  >  1:512761/1‑250 (MQ=255)                                                                                aTTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAAcc                                                                                                                                                               >  1:385885/1‑251 (MQ=255)                                                                                              gtacaattccttttggttttCATCAGTGTGCCGTTTTACTTTTTAAGCCGGGTCCATGAGATCCTGTACG‑AGAAAAATATCGTCA‑GTTCCCCTTTAAGTCAACAGACCGCGTTCCAGA‑CGGCTTCCAC‑TATCCCCCTT‑GCGGTTTACGCCACA‑ATTTTTCCCACGATGTCACCGTTGTGGTGACCGATGAACGGTTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCttt                                                                                                                                                <  1:321573/234‑1 (MQ=255)                                                                                              gTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACActtc                                                                                                                                                   <  1:595784/249‑1 (MQ=255)                                                                                                aCGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTg                                                                                                                                               <  1:310943/251‑1 (MQ=255)                                                                                                aCGAGTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCttt                                                                                                                                                <  1:136937/250‑1 (MQ=255)                                                                                                         gTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAg                                                                                                                                        <  2:100445/249‑1 (MQ=255)                                                                                                          ttCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTg                                                                                                                                      >  1:128966/1‑250 (MQ=255)                                                                                                              cTTTCATCAGTATAGAGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGcc                                                                                                                                 >  1:509716/1‑251 (MQ=255)                                                                                                                 tcatccGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGtt                                                                                                                              >  2:629229/1‑251 (MQ=255)                                                                                                                 tcatcaGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCGT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATGTTACCCACGATGTCACCGTTGTGGTGACCGATGAACTGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCTCACACTTCTTTGTTAGCAGTGATGCCGt                                                                                                                               <  2:322281/250‑1 (MQ=255)                                                                                                                   atcaGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTa                                                                                                                             >  1:359339/1‑250 (MQ=255)                                                                                                                     caGTATAGCGGTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATg                                                                                                                          <  1:371517/251‑1 (MQ=255)                                                                                                                         ataGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCAtt                                                                                                                      <  1:139494/251‑1 (MQ=255)                                                                                                                                  ttcTTTGTAAGCCGTGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCcagc                                                                                                              <  1:384964/248‑1 (MQ=255)                                                                                                                                              cGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGc                                                                                                                                                          >  1:164326/1‑194 (MQ=255)                                                                                                                                              cGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGc                                                                                                                                                          <  2:164326/194‑1 (MQ=255)                                                                                                                                                    cATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTgg                                                                                           >  1:260773/1‑251 (MQ=255)                                                                                                                                                       cAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTc                                                                                         >  1:434835/1‑250 (MQ=255)                                                                                                                                                        aGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAg                                                                                       <  2:226010/251‑1 (MQ=255)                                                                                                                                                              tGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAAtgtg                                                                                 >  1:378566/1‑251 (MQ=255)                                                                                                                                                                 aCGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAAc                                                                              <  2:645882/251‑1 (MQ=255)                                                                                                                                                                 aCGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAAc                                                                              >  2:311991/1‑251 (MQ=255)                                                                                                                                                                 aCGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAAc                                                                              >  2:345124/1‑251 (MQ=255)                                                                                                                                                                  cGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACg                                                                             >  2:509748/1‑251 (MQ=255)                                                                                                                                                                    gagaAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACgg                                                                            >  1:446879/1‑250 (MQ=255)                                                                                                                                                                       aaaa‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTc                                                                        <  1:347801/251‑1 (MQ=255)                                                                                                                                                                               cGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACAtg                                                                 <  1:659838/251‑1 (MQ=255)                                                                                                                                                                                gTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACAtgg                                                                >  1:397951/1‑251 (MQ=255)                                                                                                                                                                                                 gTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGTCCAATGACCgg                                                >  1:552283/1‑251 (MQ=255)                                                                                                                                                                                                    aaCAGACCGCGTTCCAGAACGACTTC‑ACGCA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTc                                             >  1:551660/1‑251 (MQ=255)                                                                                                                                                                                                     aCAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑AGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTcc                                            >  1:259422/1‑251 (MQ=255)                                                                                                                                                                                                             cgTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGAc                                                   >  2:317325/1‑236 (MQ=255)                                                                                                                                                                                                             cgTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGAc                                                   <  1:317325/236‑1 (MQ=255)                                                                                                                                                                                                                       cGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGc                          <  2:207026/251‑1 (MQ=255)                                                                                                                                                                                                                        gACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCt                         <  2:359339/251‑1 (MQ=255)                                                                                                                                                                                                                         aCTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTg                        >  2:541131/1‑251 (MQ=255)                                                                                                                                                                                                                                 gTA‑CACTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCt                 >  1:600455/1‑251 (MQ=255)                                                                                                                                                                                                                                     cacTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTct              <  2:543692/251‑1 (MQ=255)                                                                                                                                                                                                                                     cacTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTct              <  1:155354/251‑1 (MQ=255)                                                                                                                                                                                                                                      acTCTTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTctg             >  2:450389/1‑251 (MQ=255)                                                                                                                                                                                                                                       ctctTACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTctgc            >  2:262555/1‑251 (MQ=255)                                                                                                                                                                                                                                           tACCGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTctgct           <  1:166008/248‑1 (MQ=255)                                                                                                                                                                                                                                              cGGTTT‑CGGCACAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGccctg     >  2:192520/1‑248 (MQ=255)                                                                                                                                                                                                                                                  tt‑CGGCCCAAATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTa  <  2:552197/250‑1 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 TTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAA‑TATCGTCAAGTTCCGCTT‑CAGTCAACAGACCGCGTTCCAGAACGACTTC‑ACGTA‑CACTCTTACCGGTTT‑CGGCACAGATTTTACCCACGATGTCACCGTTGTGGTGACCGATGAACGGGTTCAGGTAAGTAACGATACCGATAGAGTTGTAAACGTAACCTTCGCACACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTA  >  NC_000913/4366792‑4367278

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