Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,556,917	A→G	E2G (GAA→GGA)	idlP →	iraD leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,556,917	0	A	G	100.0%	172.7 / NA	53	E2G (GAA→GGA)	idlP	iraD leader peptide
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/19);  total (34/19)

TAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGC  >  NC_000913/4556673‑4557163                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       taaataaTTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATg                                                                                                                                                                                                                                                   >  1:349027/1‑250 (MQ=255)  aaataaTTTCCCACCGGAAGTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAg                                                                                                                                                                                                                                                 <  1:25378/251‑1 (MQ=255)   aataatTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGc                                                                                                                                                                                                                                                >  1:510745/1‑251 (MQ=255)         ttCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAAt                                                                                                                                                                                                                                          >  1:338547/1‑251 (MQ=255)            ccACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATAc                                                                                                                                                                                                                                        >  1:518226/1‑250 (MQ=255)             cACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAAGGAGCATCAATACAg                                                                                                                                                                                                                                      <  1:499440/251‑1 (MQ=255)                      tAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTTTTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCGATTGGAAATGAGCATCAAGACAGAGGTGccc                                                                                                                                                                                                                              <  2:40546/250‑1 (MQ=255)                           agcagcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGtt                                                                                                                                                                                                                        <  2:351657/251‑1 (MQ=255)                            gcagcGGCGACAAGTATTTTCATCTACGTAACTATTCAGACTATTCAGAATATACTCGCGGTTTTTTTTTGAATTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAACAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGAATCAATACAGTGGTGCCCGGGGttt                                                                                                                                                                                                                       >  2:613779/1‑250 (MQ=255)                                   tGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCa                                                                                                                                                                                                                 >  1:353392/1‑250 (MQ=255)                                    gACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCa                                                                                                                                                                                                                 >  1:460840/1‑249 (MQ=255)                                    gACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTGTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTAAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAg                                                                                                                                                                                                               <  1:384454/251‑1 (MQ=255)                                          tATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACTGTGGTGCCCGGTGTTAAGGGCAAGCCGCat                                                                                                                                                                                                         >  2:168422/1‑251 (MQ=255)                                              ttCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGGGTTCAGGGCAAGCCGCATATGt                                                                                                                                                                                                     >  2:288490/1‑251 (MQ=255)                                                cATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAACTCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCCGGCCTGTATCAGACTGACATGGATGTAATGAAAAAAAGGGAATGGCTATGGGACATGAGCATCAATACAGTGGTGCCCGGGGTTAAGGGCAAGCCGAATATGTTg                                                                                                                                                                                                   >  1:83099/1‑251 (MQ=255)                                                                cGGAATATTCATATTATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATTAAATCATGACTTTTTTCCTGAAAAGTCAGGAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGAAATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGGTGCTAAACGTCAGGAGt                                                                                                                                                                                    <  2:424996/250‑1 (MQ=255)                                                                      atTCATAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCTCATATGTGGCTAAACGTCAGGAGTGCGCaa                                                                                                                                                                              <  1:509784/250‑1 (MQ=255)                                                                       ttCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCaaa                                                                                                                                                                             >  1:359147/1‑250 (MQ=255)                                                                         cAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGAATTTTTTCCTGAAAAGTCAGTAGATATTAATAGGCAAGTAAGGTTTTATACTTTGCTGGCAGGATTCAGGCCTGTCTCAGAATGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCAGCATATGTTGCTCAACGTCAGGAGTGCGCAAAatg                                                                                                                                                                          >  1:205322/1‑251 (MQ=255)                                                                                   tAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCAGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGGCAAt                                                                                                                                                                >  2:139975/1‑251 (MQ=255)                                                                                    aGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTATTTTACTGAAAAGTCAGTAGATAATAATCGCCAAGAAAGGATTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGAGTGACATGGATGAAATGAGCAAAAGGGAATGGCCATGGGTAATGCGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAGCGTCAGGAGTGCGCAAAATGATGCGACCATc                                                                                                                                                               >  1:51084/1‑251 (MQ=255)                                                                                        gTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTT‑CCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTGTATACTTTGCTGACAGGATTCAGGCCTGTCGCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCatat                                                                                                                                                                                                       <  1:368898/206‑1 (MQ=255)                                                                                              ttttGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCt                                                                                                                                                     >  1:302221/1‑251 (MQ=255)                                                                                                   aGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGtttt                                                                                                                                                >  2:111232/1‑251 (MQ=255)                                                                                                     tttGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGCGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTAc                                                                                                                                              >  2:156880/1‑251 (MQ=255)                                                                                                           gatTTTTATCAAATCATGACTTTTTTCCTGAAAAGGCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTGTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCAGCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGaa                                                                                                                                         <  1:127931/250‑1 (MQ=255)                                                                                                            atTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAAt                                                                                                                                       <  1:112354/251‑1 (MQ=255)                                                                                                               tttATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAAttt                                                                                                                                     <  1:288490/250‑1 (MQ=255)                                                                                                                      aaTCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAAt                                                                                                                             >  1:49897/1‑251 (MQ=255)                                                                                                                         cATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATaa                                                                                                                           >  1:451431/1‑250 (MQ=255)                                                                                                                           tGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTACTGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAACATGATGCGACAATCACTTCAGGCTGTTTTACCTGCAATTTCAGGCAATAAAAc                                                                                                                        >  1:427234/1‑251 (MQ=255)                                                                                                                                 ttttCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGACATTTCAGGCAATAAAACGTCATc                                                                                                                  >  1:22701/1‑251 (MQ=255)                                                                                                                                       tGAAAAGTCAGTAGATATTCATAGGCACGTAAGGATTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAACAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCGGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCg                                                                                                            >  2:587823/1‑251 (MQ=255)                                                                                                                                                             aGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGAtctc                                                                                      >  1:385027/1‑251 (MQ=255)                                                                                                                                                                              taCTTTGCTGCCAGAATTCAGGCCTGCATCAGACTGACATGGATGTAATGAACCAAAGGGACTGGATATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGATAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGCAATTTCAGGCAATAAAAAGTCATCACTGCGTAAATCGGTATGCAGCGATCTCATCACTCTTTTTAATtc                                                                     >  2:408483/1‑251 (MQ=255)                                                                                                                                                                                  ttGCTGACAGGATTCAGGCCTGTCTCAGACTGCCATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTctcct                                                                 >  1:212327/1‑251 (MQ=255)                                                                                                                                                                                    gCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCa                                                               <  2:158492/251‑1 (MQ=255)                                                                                                                                                                                        aCAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGATCAGGGCAAGCCGCATATGTTGCTAAACGACAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGACATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCg                                                           >  1:51287/1‑251 (MQ=255)                                                                                                                                                                                                cAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTAc                                                    <  2:376384/250‑1 (MQ=255)                                                                                                                                                                                                cAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACGCTTTTTAATTCTCCTCATTCGGCATTAc                                                    <  2:49897/250‑1 (MQ=255)                                                                                                                                                                                                         ctcAGACTGACATGGATGTAATGAACGAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATaaaa                                                                                                                         <  1:512494/172‑1 (MQ=255)                                                                                                                                                                                                         ctcAGACTGACATGGATGTAATGAACGAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATaaaa                                                                                                                         >  2:512494/1‑172 (MQ=255)                                                                                                                                                                                                           cAGACTGACATGGAGGTAATTAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGGCCGGTGTTCAGGGCAAGCCTCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATGTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTTCTCCCTCGTTTTAATCCTCCGCATTCGGCATTACCGTCATTaa                                           <  2:413914/248‑2 (MQ=255)                                                                                                                                                                                                                 gACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTc                                    >  2:151972/1‑249 (MQ=255)                                                                                                                                                                                                                      ggATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTGTACCTGAAATTTCAGGCATTAAATCGTCATCGCTGCGTAAATCGGTATGCAGCTATCTCCTCACTCTGTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCAt                              <  1:395185/250‑1 (MQ=255)                                                                                                                                                                                                                          gTAATGAACAAACGGGAATGGCTATGGGCAATGAGCATCAATACAGTGGTGCAAGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGACATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCg                         >  2:288575/1‑251 (MQ=255)                                                                                                                                                                                                                            aaTGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGATGCCCGGTGTTCAGGGCAAGCCGCATATGTTGATAAACGTCAGGAGTGCGCAAAATGATGCGCCAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAAAGTCATCGCTGCGTAAATCGGTCTGCAGAGATATCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTCCTTGTTTCAGGCATGCCCg                         >  1:416724/1‑249 (MQ=255)                                                                                                                                                                                                                             aTGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGt                      >  1:548796/1‑251 (MQ=255)                                                                                                                                                                                                                             aTGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGGTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGt                      <  2:212327/251‑1 (MQ=255)                                                                                                                                                                                                                                    aaaGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGAAATAAAACGTCATCGCTGCGTAAATCGGAATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGACTGCCCGAGTGGCAAgg               >  2:142527/1‑250 (MQ=255)                                                                                                                                                                                                                                     aaGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTg              >  1:616820/1‑251 (MQ=255)                                                                                                                                                                                                                                       gggAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCa            <  2:246537/251‑1 (MQ=255)                                                                                                                                                                                                                                                 aTGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTACATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGCGTGGCAAGTGCATAACCCGAGc  >  2:561017/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       TAAATAATTTCCCACCGGAATTAAGTAGCAGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGC  >  NC_000913/4556673‑4557163

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