Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,559,104	G→A	S141S (TCC→TCT)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,559,104	0	G	A	98.0%	156.8 / ‑6.8	49	S141S (TCC→TCT)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/23);  minor base T (1/0);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

ACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAA  >  NC_000913/4558865‑4559336                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         acaTCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATg                                                                                                                                                                                                                               >  1:217423/1‑251 (MQ=255)               cagcagGTCATAAATAGGCTGTAACGCCTTTTTACGGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACATCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCa                                                                                                                                                                                                                 <  2:265377/251‑1 (MQ=255)                agcgggTCATAAATAGGCTGAAACTCTTTTTTACGGTCGCCCAGGGGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTGTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAg                                                                                                                                                                                                                <  2:67358/247‑1 (MQ=255)                     gTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATcc                                                                                                                                                                                                           >  2:277957/1‑251 (MQ=255)                      tCATAACTAGGCTGTAACGCCTTTTTCCTGTCGACCATGTGGAACACGGTGACGCCAGCTTTACCCCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATcc                                                                                                                                                                                                           >  2:173538/1‑250 (MQ=255)                        aTAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAgg                                                                                                                                                                                                        >  1:311522/1‑251 (MQ=255)                           aaTAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGgcgc                                                                                                                                                                                                     >  1:295767/1‑251 (MQ=255)                                 cTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGc                                                                                                                                                                                               >  1:357836/1‑251 (MQ=255)                                   gTAACGCCTTTTTACTGTCGCCCATGTGGAACCCGGTGACGCAAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGc                                                                                                                                                                                               <  1:378558/249‑1 (MQ=255)                                      aCGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCACCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCcttct                                                                                                                                                                                          >  2:317509/1‑251 (MQ=255)                                        gCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCACTCACACGATCGATCATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACGTGATAAGCGCCGGTCAGCATCCAGGCGCTGCGGCcttcctc                                                                                                                                                                                        >  2:240913/1‑248 (MQ=255)                                        gCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCttcttc                                                                                                                                                                                        >  2:375279/1‑251 (MQ=255)                                            ttttACTGGCGCCCATTGGGAACCCGGTGACGCCATGTTTACCGCCGAGCAAACCGCCAGCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCAtt                                                                                                                                                                                     <  2:353631/250‑1 (MQ=255)                                             tttACTGTCGCCCATGTGGCACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCTGGCTTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCCTTTCACGCCAATCACCCGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATTGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGTCGCTGCTGCCttcttccttt                                                                                                                                                                                    >  2:290070/1‑246 (MQ=255)                                              ttACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGa                                                                                                                                                                                   <  1:300617/250‑1 (MQ=255)                                                 cTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAgcgc                                                                                                                                                                               >  2:289085/1‑251 (MQ=255)                                                      gCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACggg                                                                                                                                                                          >  1:213822/1‑251 (MQ=255)                                                         cATGTGGAACGCGGTGACGCCAGGTTTGCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGGGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCt                                                                                                                                                                       <  1:21310/251‑1 (MQ=255)                                                                aaCACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCg                                                                                                                                                                  >  2:350431/1‑249 (MQ=255)                                                                    gggTGACGCCGGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAgg                                                                                                                                                            <  2:269240/250‑1 (MQ=255)                                                                         aCGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTc                                                                                                                                                       >  2:173137/1‑251 (MQ=255)                                                                              aggtttcacgccgAGCACACCGACAACGCGGGTTGCTGCCGCCATATTGGCCAGGCGATAAACGTCCGGTGCGGCCGCACGGAGATCAGAGATCGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGCACCCGTAGTGGTGCGAGAAGGGACATGATAAGCGACGGTCAGCATCCAGGCGCTGAAGCCTTCGTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTTCgggg                                                                                                                                                  >  2:62467/9‑250 (MQ=255)                                                                               ggTTTCCCGCCGAGCAAACCGCAAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTg                                                                                                                                                 <  2:333493/251‑1 (MQ=255)                                                                                     ccgccgAGCAAACCGCCAACGCGGGATTCTGACGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCg                                                                                                                                              >  2:24785/1‑248 (MQ=255)                                                                                        ccgAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCgag                                                                                                                                            >  1:255826/1‑247 (MQ=255)                                                                                               aaCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAGTCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCGTCATTTAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCg                                                                                                                                 <  2:161586/251‑1 (MQ=255)                                                                                                   gCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGc                                                                                                                             >  2:268134/1‑251 (MQ=255)                                                                                                      aaCGCGGGATTCTGCAGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATACGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGCGATAGAGTCGGGGccc                                                                                                                           >  2:347633/1‑250 (MQ=255)                                                                                                         gcgGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCcagc                                                                                                                        >  1:250835/1‑250 (MQ=255)                                                                                                                    ccgccATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGTGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGCGCCCAGCAGACCAACCACt                                                                                                            >  2:261101/1‑251 (MQ=255)                                                                                                                        cATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGCTTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtgacg                                                                                                        >  2:303553/1‑251 (MQ=255)                                                                                                                          tatTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACtgacgtg                                                                                                      >  2:350016/1‑251 (MQ=255)                                                                                                                              ggCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAAACCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGc                                                                                                  >  2:307479/1‑251 (MQ=255)                                                                                                                                    gtgcgaaaCGTTCGGGGCGGCAGAACGGTGATCAGAGATGGCGCACTTCACGCCAATCACACGAGCGATAATCGCCACGGCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTGCATTGAGCGCACGGGTCTTGGCGAGCAGGGAGTCCGGGCGGGGTGCGATAGTGACGGTGCCCAGCAGACGATCCACtgaggcggcacccgcct                                                                                            <  2:246959/246‑15 (MQ=255)                                                                                                                                            cGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCAAATCACACGATAGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAgg                                                                                    >  1:214641/1‑251 (MQ=255)                                                                                                                                               ccGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGa                                                                                 <  1:277957/251‑1 (MQ=255)                                                                                                                                                        cAGAACGGTGATCAGAGATGGCGCAGTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTGAGCGcc                                                                        <  2:337152/251‑1 (MQ=255)                                                                                                                                                                    cAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGATCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCTACTTAGCGCCACTTCCGGCGTg                                                            <  2:123457/251‑1 (MQ=255)                                                                                                                                                                             cgcATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTg                                                   >  2:333949/1‑251 (MQ=255)                                                                                                                                                                                         cAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcg                                       <  1:362017/251‑1 (MQ=255)                                                                                                                                                                                          aaTCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgc                                      <  2:187911/251‑1 (MQ=255)                                                                                                                                                                                           aTCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACGTCCGGCGTGCGCGTCGTGGGACCTGCTTcgcc                                     <  2:143011/251‑1 (MQ=255)                                                                                                                                                                                                             gCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGGGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATg                   <  2:320533/251‑1 (MQ=255)                                                                                                                                                                                                               cACGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCGCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGa                 <  2:250835/251‑1 (MQ=255)                                                                                                                                                                                                                 cGTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACg               <  1:256297/251‑1 (MQ=255)                                                                                                                                                                                                                  gTCTTTTTCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACgt              <  1:63512/251‑1 (MQ=255)                                                                                                                                                                                                                       tttCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTg          <  1:160111/250‑1 (MQ=255)                                                                                                                                                                                                                       tttCCACGGAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTg          <  1:16447/250‑1 (MQ=255)                                                                                                                                                                                                                              ggAACCCGTAATGGTGCGAGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATaa  <  1:350016/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                         ACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAA  >  NC_000913/4558865‑4559336

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