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	95.6%	128.5 / ‑3.6	45	E2G (GAA→GGA)	idlP	iraD leader peptide
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (32/11);  total (32/13)
Fisher's exact test for biased strand distribution p-value = 7.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

AGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGC  >  NC_000913/4556702‑4557163                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        agcGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTc                                                                                                                                                                                                                        >  1:527929/1‑249 (MQ=255)    ggTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATCTACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATAATGACTTTTTTCCTGAAAAGTCAATAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGAATGACATGGCTGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCAACAATAAAGTGGTGCCCGGTGTTCAggt                                                                                                                                                                                                                    >  2:598770/1‑249 (MQ=255)     gTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGc                                                                                                                                                                                                                   >  1:67253/1‑250 (MQ=255)          aaGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTGTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCg                                                                                                                                                                                                             <  2:588851/251‑1 (MQ=255)                     tCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTAATGATTTTTATCAAATCCTGACTTTTTTCCTTAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGCACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCCGGGCCAGCCGCATATGATGCt                                                                                                                                                                                                  >  1:632854/1‑251 (MQ=255)                         aggaaCTATTCGGAATAATCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATGTTTATCAAATCATGATTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCGGTATCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGGGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAAc                                                                                                                                                                                              <  2:311472/248‑1 (MQ=255)                                   cGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCTGGAGTg                                                                                                                                                                                    >  2:638716/1‑251 (MQ=255)                                     gAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTgcg                                                                                                                                                                                  <  1:255697/251‑1 (MQ=255)                                                        gCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATc                                                                                                                                                                >  1:116008/1‑250 (MQ=255)                                                                        tttGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAACTGAGCATCAATACGGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGGGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTAc                                                                                                                                               >  2:511448/1‑251 (MQ=255)                                                                               atTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAAt                                                                                                                                        >  1:265424/1‑251 (MQ=255)                                                                                     aTCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAg                                                                                                                                   <  1:220223/250‑1 (MQ=255)                                                                                      tCAAATCATGACTTTTTTCCTGAAAAGTAAGTAGATATTCATAGGAAATTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAgg                                                                                                                                  <  2:116008/250‑1 (MQ=255)                                                                                                    ttttCCTGAAAAGTCAGTAGATATTCATAGGCCAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCCATACAGTGGTGCCCGGTGTTCAGGGc                                                                                                                                                                                                                   >  2:469864/1‑155 (MQ=255)                                                                                                    ttttCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGc                                                                                                                                                                                                                   <  1:469864/155‑1 (MQ=255)                                                                                                      ttCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGCCAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGc                                                                                                                 >  2:336487/1‑251 (MQ=255)                                                                                                       tcctgccaagtccttctctatTCATAGCCAAATAAGGTTTTATACTTTGCTCACAGGATTCAGGCCTGTCTCAGACTGACATGTATGTAATGAACAAAAGGGAGTGGCTATAAGAAATGAGCATCAATACAGTGGGGCCCGGTGTTCAGAACAAGCCGCATATGTTGCTAGACGTCGGGATTGCGCGAAATAATGCGACAATCCCTTCAGGCTGTTTTACCTGAACTTTCAGGCAATAAACCGTCATCGCt                                                                                                                >  1:428662/19‑251 (MQ=255)                                                                                                                      aGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGCCAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGCGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGtctg                                                                                                 >  2:399144/1‑248 (MQ=255)                                                                                                                        atatTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCata                                                                                                                                                                                                         <  1:644951/145‑1 (MQ=255)                                                                                                                        atatTCATAGGCAAATAAGGTTTTATACTTTGCTGACAGGATTCAGGCATGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCata                                                                                                                                                                                                         >  2:644951/1‑145 (MQ=255)                                                                                                                             cATATGCAAGTAAGTTTTTATACTTTGCTGACATGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTCCCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGAt                                                                                         <  2:519201/251‑1 (MQ=255)                                                                                                                             cATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGAt                                                                                         <  1:150571/251‑1 (MQ=255)                                                                                                                                      gTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAgg                                                                                                                                  >  1:86744/1‑202 (MQ=255)                                                                                                                                       tgaagTTTTATACTTTGCTCACAGGATTCAGGCCTGTCACAAACAGACATGGATGTAATGTACAAAAGAGAATGGATATGGGAAATGAGCATCAATACATTGGAGCCCGAAGTGCAGGACAAGAAGAATATATTGCTAAACGACAGGAGTGCGCAAAAGCATGCGACAATCACTTCAGGCTGTTTTACCTAAAATTTCAGGCAATCAAACATCATCTATGCGAAAATCGGTTT‑GCATCGAtaacatcact                                                                                >  2:434823/5‑241 (MQ=255)                                                                                                                                       tAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGAGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCActc                                                                               >  1:263547/1‑251 (MQ=255)                                                                                                                                            ttttATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCttttt                                                                          >  1:655498/1‑251 (MQ=255)                                                                                                                                                        cTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAGGCGGCATCGGTTGTTAACGGTCAGGAGTGCGCAAAATGTTGCGCCAGTCAGTTCGGGCGGTTTTCCTTGAATTTTCGGGCATTAAAACTTCATCCCTCCTTAATTCG‑TATTGCAGGGATCCCTTCACTTTTTTTAATtctcttctt                                                              >  2:268067/1‑246 (MQ=255)                                                                                                                                                        cTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCAt                                                              >  2:232715/1‑251 (MQ=255)                                                                                                                                                           aCAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAACCGTCAGGAGTGCGCAAACTGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATATCCTCACTCTTTTTAATTCTCCTCATTCg                                                           >  1:328355/1‑251 (MQ=255)                                                                                                                                                           aCAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCg                                                           >  1:30251/1‑251 (MQ=255)                                                                                                                                                                    aGGCCTGTCTCAGACTGACATGGATGTAATTAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGATCAGGGCAAGCCGCATATGTTGCTAAACGTCCGGATGGCGCAAAATGATGCGACAATCACTTCAGGCAGTTTTACCTGAAATTTCAGGCAATAAAACATAATCGCTGCGTAAATAGGTAT‑GCAGAGATCTCCTCACTCTTTTTATTTCTCTTCATTCGGAATTaacg                                                  >  2:362839/1‑248 (MQ=255)                                                                                                                                                                        cTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCAt                                              >  2:333780/1‑251 (MQ=255)                                                                                                                                                                            ctcAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGATCAGGGCAAGCAGCATATGTTGCTAAAAGTCAGGAGTGCGCAAAATGATGCGACAAACACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGAAAATCGGAAT‑GAAGCGATCTCCACACTCTTTTTAATTATCCTCATTCGGCATTACCGTCAttcct                                          >  2:357179/1‑248 (MQ=255)                                                                                                                                                                               aGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGt                                       >  1:153387/1‑251 (MQ=255)                                                                                                                                                                                          gATGTAAGGAACAAAAGGGAATGGCTTGGGAAAATGTGCAGAAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGGTGCTAAACGTCAGGAGTGAGCAAAATGATGCGACAATCACTTCAGGCTGGTTGACCTGAAATTTCAGGCAATAAAACGGCATCGCTGCGTAAATCGGTAT‑GCAGTGATCTCCTCACTCTTTTTAAGTCTCCTCATTCGGCAGTACCGTCATTACTTGTTTCAGGCaggc                            <  2:246623/251‑4 (MQ=255)                                                                                                                                                                                           aTGTAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGcc                           >  1:318027/1‑251 (MQ=255)                                                                                                                                                                                              tAATGAACAAAAGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCg                         >  1:345384/1‑250 (MQ=255)                                                                                                                                                                                                      aaaaGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAg                <  1:516367/251‑1 (MQ=255)                                                                                                                                                                                                      aaaaGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCACGAGTGGCa                  >  1:478790/1‑249 (MQ=255)                                                                                                                                                                                                       aaaGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGt               >  1:444354/1‑251 (MQ=255)                                                                                                                                                                                                        aaGGGAATGGCTATGGGAAATGAGCATCAATTCAGTGGTGCCCGGTTTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTg              <  2:265424/251‑1 (MQ=255)                                                                                                                                                                                                         aGGGAATGGCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGc             >  1:315621/1‑251 (MQ=255)                                                                                                                                                                                                                ggCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAACATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGAAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGtt                                      >  2:248511/1‑219 (MQ=255)                                                                                                                                                                                                                ggCTATGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGtt                                      <  1:248511/219‑1 (MQ=255)                                                                                                                                                                                                                    aTGGGAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGc  >  2:640949/1‑251 (MQ=255)                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        AGCGGTGACAAGTATTTTCATCTACGTAACTATTCGGAATATTCAGAATATACTAGCGGTTTTTTTTTGAGTTTGATGATTTTTATCAAATCATGACTTTTTTCCTGAAAAGTCAGTAGATATTCATAGGCAAGTAAGGTTTTATACTTTGCTGACAGGATTCAGGCCTGTCTCAGACTGACATGGATGTAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAT‑GCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAGC  >  NC_000913/4556702‑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: ‑