Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,123,823	2 bp→AG	coding (231‑232/648 nt)	grxB ←	reduced glutaredoxin 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,123,823	0	G	A	98.1%	188.8 / NA	52	L78L (CTG→TTG)	grxB	reduced glutaredoxin 2
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (23/28);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	1,123,824	0	T	G	98.1%	194.9 / ‑3.1	52	coding (231/648 nt)	grxB	reduced glutaredoxin 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (23/28);  minor base . (1/0);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAGTAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACGT  >  NC_000913/1123580‑1124056                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                         cGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGt                                                                                                                                                                                                                                    >  1:580132/1‑251 (MQ=255)   aTCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGttt                                                                                                                                                                                                                                  <  2:413900/251‑1 (MQ=255)      aGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGcc                                                                                                                                                                                                                               <  1:196373/251‑1 (MQ=255)            cccaggagtcGTAAATCAGCGCTGAAATTCTTAATCAGACCGTCAGGGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCTTCTTTCTTGGCGCCGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGGCCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAg                                                                                                                                                                                                                         <  2:307412/242‑1 (MQ=255)             ccAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGCGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAACCTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAAAGATTGCCGG‑TCAAGAACGGTTTGCCggcgcgt                                                                                                                                                                                                                        >  1:137905/1‑245 (MQ=255)                  gcACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTCCCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATc                                                                                                                                                                                                                   <  2:478083/251‑1 (MQ=255)                         aaTCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGGAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCACCAGCAGTTTGTTGGCGAAGCCCTTGACCTTGCGCCGCCACTCTTCACTGGCAGGGGAACGTTGGCCGG‑TCAAGGACGGGTTGCCGTCGAGGTTATCGACATa                                                                                                                                                                                                             >  2:505773/1‑250 (MQ=255)                          atcatgGCTGATATTCTTAATCAGACCGTCAGAGTGGGCAAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGAAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTGGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGt                                                                                                                                                                                                           <  2:339049/251‑1 (MQ=255)                            catcGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTg                                                                                                                                                                                                          <  1:563015/250‑1 (MQ=255)                                           tAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACCCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATg                                                                                                                                                                                             <  2:442489/248‑1 (MQ=255)                                              tCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGGGTAGCCATTGACCTTGCGCAGCCACTCTTCCATGGCAGGGGAGAGTTTGCCGGGCCAG‑AACGGTTTGCCGGCGAGTTTATTGACCTAGTGAACGATGttgatgccttct                                                                                                                                                                                       >  2:315718/1‑240 (MQ=255)                                                        agagTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAAAAGCAGTTTGTTGGCGTAGCCCTTGAACTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGCACGATGTCCATGCTTTCTGGCATCTAGc                                                                                                                                                                             >  1:234600/1‑251 (MQ=255)                                                                  agcagGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATc                                                                                                                                                                   >  1:396975/1‑251 (MQ=255)                                                                        tCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCtttttt                                                                                                                                                             >  2:172372/1‑250 (MQ=255)                                                                            cAAAATTCCCCGGGCTGGCCTCTTTCTTTGCGACGAAATATTTGCGCGCGGCGGGATTAGAAATCTCATCAAAGGCCGATTTGGCAAACCGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGa                                                                                                                                                         <  2:299740/251‑1 (MQ=255)                                                                                    cccGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAAAGATGTACATGCTTTCTGGCATATAGCGGCTGTCCTCTTTTTGCAGAATGGGAAc                                                                                                                                                 >  2:413209/1‑251 (MQ=255)                                                                                       gcgcTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTg                                                                                                                                              <  1:34980/251‑1 (MQ=255)                                                                                               ctctTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCCTATAGCGGCTGTCATCTTTTTGCAGAATGGGGACCTGTTTTTgcc                                                                                                                                      >  2:176106/1‑249 (MQ=255)                                                                                                   ttcttGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTgacc                                                                                                                                     <  2:91029/248‑1 (MQ=255)                                                                                                          cgacgaAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGTAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACTATGTCCATGCTTTCTGGCATATAGCGTCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCg                                                                                                                            <  1:363468/250‑1 (MQ=255)                                                                                                          cgacgaAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCAATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCgt                                                                                                                           <  2:304629/251‑2 (MQ=255)                                                                                                                        cgcgcgGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGc                                                                                                             <  1:151229/251‑1 (MQ=255)                                                                                                                         gcgcgGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGcg                                                                                                            <  1:441329/251‑1 (MQ=255)                                                                                                                              gcggGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGTAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGcgtcgtc                                                                                                       >  1:14169/1‑251 (MQ=255)                                                                                                                                 ggAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGtt                                                                                                    >  1:195879/1‑251 (MQ=255)                                                                                                                                    gTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTg                                                                                                   >  1:134120/1‑249 (MQ=255)                                                                                                                                         aaaCTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGa                                                                                             >  1:226252/1‑250 (MQ=255)                                                                                                                                         aaaCTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACAGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGa                                                                                             >  1:127010/1‑250 (MQ=255)                                                                                                                                              catcaAATGCCGATTGGGCAAAACGCGGCAACAGCAGTTGGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGt                                                                                       <  2:432387/251‑1 (MQ=255)                                                                                                                                              catcaAATGCAGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGt                                                                                       >  1:144044/1‑251 (MQ=255)                                                                                                                                               atcaAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTa                                                                                      <  1:550368/251‑1 (MQ=255)                                                                                                                                                 caAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCCTCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAAt                                                                                    >  2:306376/1‑251 (MQ=255)                                                                                                                                                 caAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTaa                                                                                     <  1:496876/250‑1 (MQ=255)                                                                                                                                                    aTGCCGATTTGGCAAAACGCGGCATCAGCAGTTTTTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGT‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCg                                                                                 <  1:367367/251‑1 (MQ=255)                                                                                                                                                             tGGCAAAACGCGGCAACGGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGat                                                                         <  2:433214/250‑1 (MQ=255)                                                                                                                                                                aaaaaCGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATAttt                                                                     <  2:344552/250‑1 (MQ=255)                                                                                                                                                                        ggCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGcc                                                             >  2:411429/1‑251 (MQ=255)                                                                                                                                                                        ggCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGcc                                                             <  1:210903/251‑1 (MQ=255)                                                                                                                                                                              agcagTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCCTCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGGTTCTGAGTCGTCGTTGAGCAGAACATGTAAATCGACGGGGATATTTTTCAGGCCGAAAAt                                                       >  2:322639/1‑251 (MQ=255)                                                                                                                                                                                 agTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCAt                                                    >  2:471649/1‑251 (MQ=255)                                                                                                                                                                                          gCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCttt                                           <  1:251272/251‑1 (MQ=255)                                                                                                                                                                                               gCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTGCTGCGTCGTCGTTGAGCAGAACATGTAAGTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGc                                      <  1:213600/251‑1 (MQ=255)                                                                                                                                                                                                     gACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAg                                <  1:84740/251‑1 (MQ=255)                                                                                                                                                                                                             gcAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTg                         >  1:94759/1‑250 (MQ=255)                                                                                                                                                                                                                    ctctTCAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTaa                  >  1:461858/1‑250 (MQ=255)                                                                                                                                                                                                                         cAATGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGta             <  1:172372/250‑1 (MQ=255)                                                                                                                                                                                                                          aaTGGCCGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGTGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGGGATCTTAAATGtata           <  2:533121/251‑1 (MQ=255)                                                                                                                                                                                                                          aaTGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGtata           <  2:388292/251‑1 (MQ=255)                                                                                                                                                                                                                            tGGCAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGTCCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGc         >  2:400096/1‑251 (MQ=255)                                                                                                                                                                                                                               cAGGGGAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTc      >  2:125621/1‑251 (MQ=255)                                                                                                                                                                                                                                  gggAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGCCGGGGATATTTTTCAGGCCGAACATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACg   >  2:35169/1‑251 (MQ=255)                                                                                                                                                                                                                                   ggAACGTTTGCCGG‑TCAAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACGt  <  1:214798/251‑1 (MQ=255)                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                         CGATCAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGG‑TCAGTAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGGGTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACGT  >  NC_000913/1123580‑1124056

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