Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,215,670	G→T	D65Y (GAT→TAT)	aceB →	malate synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,215,670	0	G	T	94.9%	170.2 / 1.1	60	D65Y (GAT→TAT)	aceB	malate synthase A
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (31/25);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

GTTCACAGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTGATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAGTTACA  >  NC_000913/4215433‑4215916                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         gTTCACAGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑Cgg                                                                                                                                                                                                                                             >  1:648621/1‑250 (MQ=255)           gggAAGTTTTCGGATCCATGACGAGAAGCTGCACGCTGACTGAACAGGCAACAACACCCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCGGCGCGCATTCAGACGCAGCAAGATCTTGATCACGGAACGTTGCCTTATTTTATTT‑CTGAAACAGCTTcc                                                                                                                                                                                                                                  >  2:327163/1‑251 (MQ=255)                    tCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGAt                                                                                                                                                                                                                         >  1:120370/1‑251 (MQ=255)                        aTCCATGACGCGGAGCTGCACGATGACTGAACAGGCAACAACAACCGCTGAACTGGCTTTCACAAGGCCGTATAGCGAGCAGAAGAAGCAAATTCTTACTGCCGAAGCGGAAGAATTACTGACTTAGCTGGTGCCGCATTTTACGCCACAAAGCAATAAACTTCTGTCAGAGAGCGTTAAGCAGCAGCAGGATATTGTTAACGGAACGTGGCCT‑ATTTTATTTACCAAAACACCTTCCTTTCGCGATgttg                                                                                                                                                                                                                     >  2:370685/1‑248 (MQ=255)                         cccATGACAAGGCGCTTCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCCCAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCTGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACTCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATAGTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGa                                                                                                                                                                                                                    <  2:31164/250‑1 (MQ=255)                                        ttcACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTAGGGCGAGAAGGAGAAGCAAATTCTTACTGCCGAAACGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTATGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCATCCATTAGCGAGGCTGTTTGGAAAATTCGCgg                                                                                                                                                                                                     >  2:355633/3‑251 (MQ=255)                                             gATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTc                                                                                                                                                                                               >  1:338151/1‑251 (MQ=255)                                                gACTGAACAGGCAACAACAACCGATGAACTGGCTGTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTg                                                                                                                                                                                            <  2:331618/251‑1 (MQ=255)                                                 aCTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTAGCGG‑GATTCCTGc                                                                                                                                                                                           >  2:487619/1‑251 (MQ=255)                                                       cAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACtt                                                                                                                                                                                     <  1:549623/251‑1 (MQ=255)                                                        atgCAACAACAAACGATGAACTGGCTTTCACAAGGCCGTAGGGAGAGCAGGAGAAACAAATTCTTACTGCCGAAGAGGTAGAATGTCTGACTGAGCGGGTGACGCATTTTACGCCACAACGGAATAATATTCTAGCAGCGCGCATTCCGCAGCAGCAAGATATTGATAACGGAACGTTGACTTATTTTATTT‑CGGAAACAGCTTACCTTCGCGATGATTATAGGAAAATTAGCGGTGAT‑CCTGCGGACttt                                                                                                                                                                                    >  2:188659/3‑250 (MQ=255)                                                               aacaacCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAAAAGCAACTTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGATGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCCGCGCGCATTCAGCAGCAGCAAGATATTGATAAAGGAACGTTGCCTTATTTTATTT‑CGGACACCGCTTCCATTCGCGATGCTGCTTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAACAcc                                                                                                                                                                             >  2:48167/1‑251 (MQ=255)                                                                  tacCGATGAACTGGCTTCCACAAGGCCGTATGGCGAGCAGGAGACGCAAATTCTTACTGCCGAAGAGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGAGATGCTGATTGGAAAATTCGCGG‑GATTACTGCTGCCTTAGAAGAcagcc                                                                                                                                                                          >  2:214853/2‑251 (MQ=255)                                                                  aacCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGCGAAGCAAATTCTTACTGACGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAACCAGCTTCCATTCGCGATGATGATTGGAAATTTCGCGG‑GCTTCCTGCGGACTTAGAAGAccgcc                                                                                                                                                                          >  2:258024/1‑251 (MQ=255)                                                                         gAACTGGCTTTCACAAAGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGCTAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGt                                                                                                                                                                      >  2:72624/1‑248 (MQ=255)                                                                           aCTGGCTTTCACAAGGCAGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGAAGAATTTCTGACTGAGCTGGAGACGCATTTTACGACACAACGCAATAAACTTCTGGCAGCGCGCATTCAACAGCAGCAAGGTATTGATAACGGAACGTAGCCTTATTTTATTT‑CGGAAACAGCTTCCGTTCGCGATGCTGTTTAGAAAATTAGCGG‑GATTCCGCCGGACTTAGAAGACAGCCGGGTaggg                                                                                                                                                                  >  2:468964/1‑248 (MQ=255)                                                                              ggCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACGGCCGAAGCGGTAGAATTTCTGGCTGAGCTGGTGACGCATTTTACGCCACAACGCAATAATCTTCTGGTAGCGCGCATTCAGCAGCAGCAAGAGATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGTTTCCAGTCGCGATGTTGATTGTAAAATTCGCGG‑GATTCCTGCGGACTTagaagaacggggcttagaga                                                                                                                                                                 <  2:651776/248‑15 (MQ=255)                                                                                    cacaATGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGcc                                                                                                                                                        <  1:302026/251‑1 (MQ=255)                                                                                          gCCGTATGGCAAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGACTTTCTGACTGAGCTGGTGACGCCTTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCACGATATTGCTAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCGGCGGACTTAGACGACCAACGCGTAGAGATAACTGGcacgata                                                                                                                                                   >  2:510687/1‑244 (MQ=255)                                                                                             gTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTGATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAgc                                                                                                                                               <  1:91318/251‑1 (MQ=255)                                                                                                                       ttACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGCTTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGACGACCGCCGCGTAGAGATAACTGGCCCg                                                                                                                                                      >  2:188105/1‑218 (MQ=255)                                                                                                                       ttACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCg                                                                                                                                                      <  1:188105/218‑1 (MQ=255)                                                                                                                          cTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATACCTGGCCCGGTAGAGCGCAAGCTGGTGATCAACGCGCTCAACGcc                                                                                                                  >  2:323478/1‑251 (MQ=255)                                                                                                                                     gTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTCTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGACAATGTGAAAGt                                                                                                       >  1:617907/1‑251 (MQ=255)                                                                                                                                       aGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCt                                                                                                     <  1:292398/251‑1 (MQ=255)                                                                                                                                       aGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGAAAGATATTGATAACGGAACGTTGCCTGATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCATGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCAGGTAGAGCGCAAGATGGTGATCAACGAGCTCAACGCTAATGTGAAAgact                                                                                                     >  2:304704/1‑248 (MQ=255)                                                                                                                                       aGAATTTCTGACTGAGCAGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGCCTTAGAAGACCGCCGAGTAGAGATAACTGGCCCGGGAGAGCGCAAGTTGGTGATCAACGCGCTCAACGACAAGGTGAAAGTCt                                                                                                     >  2:292398/1‑251 (MQ=255)                                                                                                                                          attgctgactgaGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTG‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTa                                                                                                  <  1:48167/247‑1 (MQ=255)                                                                                                                                                 actgaGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGa                                                                                           <  1:487619/251‑1 (MQ=255)                                                                                                                                                    gaGCTGGTGACGCATTTTACGCCACAACGAAATAAACTTCTGGCAGCGCGCATTCAGAAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCGGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTAGGGCCGAt                                                                                          <  2:288081/249‑1 (MQ=255)                                                                                                                                                     aGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTc                                                                                       <  1:327163/251‑1 (MQ=255)                                                                                                                                                         ggTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGAc                                                                                                                                                                              >  2:23027/1‑160 (MQ=255)                                                                                                                                                         ggTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGAc                                                                                                                                                                              <  1:23027/160‑1 (MQ=255)                                                                                                                                                           ggccgaaTTTTACGCACCAACGCAATAACCTTCTGGTAGCGCGCATTCAGCAGCCGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGTTTCCATTCGCGATGTTGATTGGAAAATTCGCGG‑GATTCCTGCGGAATTAGAAGACCGCCGCGTAGAGATTACTGGCCCGGTAGATCGCAAGAGGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCGTTATGTCCGATTTCGAAGAt                                                                                 <  2:414782/245‑1 (MQ=255)                                                                                                                                                            gACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAAcgc                                                                                                                            <  2:331119/207‑1 (MQ=255)                                                                                                                                                            gACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAAcgc                                                                                                                            >  1:331119/1‑207 (MQ=255)                                                                                                                                                             aCGCATTTTACGAACCAAAGCAATAACATTCTTGAAGCGCGCATTAAGCAGCAGCAAGATATTGATCACGGAACGTTGCCTTATTTTATTT‑CGGACACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCAGCGAAGAGATAACTGGCACGGTAGAGCGAAAGATGGTGATCAACGCGCTCAACGCCACTGTGAAAGTCTTTATGGCCGATTACGAAAATTc                                                                               >  2:437526/1‑251 (MQ=255)                                                                                                                                                                cATTTTACGCCGCAACGCAATAAACTTCTGGCAGCGGGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTGATTTTATTG‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCTCTAAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACt                                                                            <  2:272216/251‑1 (MQ=255)                                                                                                                                                                 aTTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGCTATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGAc                                                                                                                                                                              >  2:414258/1‑152 (MQ=255)                                                                                                                                                                 aTTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGAc                                                                                                                                                                              <  1:414258/152‑1 (MQ=255)                                                                                                                                                                          cacaACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCATCAAGATAGTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCTCTGGCACCAGAc                                                                  <  2:316366/251‑1 (MQ=255)                                                                                                                                                                          cacaACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATGCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGAc                                                                  <  1:583483/251‑1 (MQ=255)                                                                                                                                                                          cacaACGCAATAAACTTCTGGCAGCGCGAATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGAc                                                                  >  1:23869/1‑251 (MQ=255)                                                                                                                                                                             aaCGCAATAAACTTCTGGAAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGAACGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGCAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTgg                                                               >  1:328115/1‑251 (MQ=255)                                                                                                                                                                                     aaaCTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CTGAAACAGCTTCCAt                                                                                                                                                                                                                                <  2:111228/83‑1 (MQ=255)                                                                                                                                                                                     aaaCTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCAt                                                                                                                                                                                                                                >  1:111228/1‑83 (MQ=255)                                                                                                                                                                                      aaCTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTg                                                      >  1:394729/1‑251 (MQ=255)                                                                                                                                                                                               cAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAGCAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACggg                                             >  1:327717/1‑251 (MQ=255)                                                                                                                                                                                                 gcgcgcATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGc                                            >  1:467936/1‑250 (MQ=255)                                                                                                                                                                                                     gcattatgcagcagcaAGATATTTATAACGGAACTTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTTAAAGTCTTTATTGCCGATTTCGAAGATTCACTGGCCCCAGACTGGAACAAAGTGATCGACGGGCAAAtt                                       <  2:328115/244‑1 (MQ=255)                                                                                                                                                                                                      cATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATGCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTa                                      <  2:338151/251‑1 (MQ=255)                                                                                                                                                                                                            gcagcagcaAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGCCCGCCGCGTAGCGATAACTGGCCCGGTAGAGCGAAAGATGGTGCTCAACGCGCTCAACGCCAATGCGCAAGTCTTTATGGCAGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAATGATCGCCGAGAAACTTAACCTGc                                >  2:350042/1‑251 (MQ=255)                                                                                                                                                                                                                 agcaAGATATTGATAACGGAACGTTGCATTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGAt                           <  2:467936/251‑1 (MQ=255)                                                                                                                                                                                                                  gcaAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTGTGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTgctg                          >  1:57108/1‑248 (MQ=255)                                                                                                                                                                                                                  gcaAGATATTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTGTGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATg                          <  2:57108/251‑1 (MQ=255)                                                                                                                                                                                                                       atatTGATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGtt                     <  2:120135/251‑1 (MQ=255)                                                                                                                                                                                                                            gATAACGGAACGTTGCCTTATTTTATTT‑CGGACACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGCCTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGGTTTAGAAGTTTCAATGGCACCAGTCTGGAACAAGGTGATAGACGGGAAAATTAACCTGCTTGATGCGGTTAACgg                >  2:406450/1‑251 (MQ=255)                                                                                                                                                                                                                            gATAACGGAACGTTGCCTTATTTTATTT‑CGGAAACAGCTTCCATGCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACg                 <  2:477237/250‑1 (MQ=255)                                                                                                                                                                                                                                          gCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAGTTaca  <  1:178169/251‑1 (MQ=255)                                                                                                                                                                                                                                          gCCTTATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAGTTac   >  2:178169/1‑250 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         GTTCACAGTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTGATTTTATTT‑CGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGG‑GATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAGTTACA  >  NC_000913/4215433‑4215916

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