Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	15,174	G→A	R336K (AGG→AAG)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	15,174	0	G	A	92.3%	102.9 / ‑1.2	40	R336K (AGG→AAG)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (15/21);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTA  >  NC_000913/14936‑15405                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         gAGCGTGAAGGCAACAACCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAgct                                                                                                                                                                                                                              <  1:162003/251‑1 (MQ=255)       gatggCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGAGCTGGGGGGGGGAATCGAAGTACCGTCCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAAACCGGATAGCTATTCCGGATGGGCAGTAAAAGAGGCCAGTCTGGCCGCGGGGGGGCACCGGG‑GAGGTTGCTGGGCCGCGTTGTCGTCTAAACACCGGTAGGCCGGTACG‑AAAgggaggaaaaagtgcggc                                                                                                                                                                                                                         >  2:190544/4‑234 (MQ=255)             aacaacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTgctgcaagagct                                                                                                                                                                                                                  <  1:148881/251‑1 (MQ=255)             aacaacCTCTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAAGCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTgctgcaagagct                                                                                                                                                                                                                  <  1:95063/251‑1 (MQ=255)                              aGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTg                                                                                                                                                                                                 <  1:236094/251‑1 (MQ=255)                                                   ggcggcGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCacaaca                                                                                                                                                                            <  1:229422/251‑1 (MQ=255)                                                         gcTGGGTGGCGAAATCGAAGTACCGACCCATGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTACCGTATGCGCGGAAAAGGCGTCAAGTCGGTCCGCGGTGGAGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTAGAAACACCGGTGAGCCTGAACGTGCAGAC‑GCAACAGCGGCTGCAAATACTGCAAGAACGCTTCGGTGGCCCAACCgtagagcaaaaccgccagc                                                                                                                                                                      >  2:362496/1‑233 (MQ=255)                                                              gTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCa                                                                                                                                                                 <  1:182578/251‑1 (MQ=255)                                                              gTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCa                                                                                                                                                                 <  1:321470/251‑1 (MQ=255)                                                                     aaTCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGct                                                                                                                                                         <  1:315207/251‑1 (MQ=255)                                                                        cGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGGCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGcttc                                                                                                                                                       <  1:221192/250‑1 (MQ=255)                                                                        cGAAGTACCGAACCTTGATGGTCGCGGCAAACTGACAGTGCCTGGCGAAACCCAGACCGGTACGCTTTTCCGTATGCGCGGTACAGGCGTCAAGTCTGTCCGCGGTGCCCCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGCAGCGGAAACCGCTGCTGCAAGAGCCGCACGAAAGCTTCGGGGGCCCAACCGGAGAGCACCACACCCCGCGCTCC‑Aagacatact                                                                                                                                                      >  2:350922/1‑245 (MQ=255)                                                                             tACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTgctgca                                                                                                                                                                                                                        >  1:372293/1‑181 (MQ=255)                                                                              aCCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGAt                                                                                                                                                 <  1:239708/250‑1 (MQ=255)                                                                               ccggcgcTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTACGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGTCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGATGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCAAAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATgg                                                                                                                                               <  1:98661/245‑1 (MQ=255)                                                                                     cTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCAGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGgtgaga                                                                                                                                          >  1:239972/1‑247 (MQ=255)                                                                                             tCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑GGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGTTCA‑AAAAACTTATTTGATGGGGGgaagaacgatc                                                                                                                                 >  2:91296/1‑246 (MQ=255)                                                                                             tCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTgaagaa                                                                                                                                      <  1:91296/246‑1 (MQ=255)                                                                                              cgcttCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCCCAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAACCAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGGGGCCCAACCGGTGAGCACAAACGCCCGCGCTCA‑CAGAGCTTCTTTGATGGTGTGGATAAGttttt                                                                                                                                >  2:121477/5‑251 (MQ=255)                                                                                              cgcgTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGtttt                                                                                                                                 <  1:251252/250‑1 (MQ=255)                                                                                                       cTGAAAGTGCCTGGAGAAACCCAGACCGGTAAGCTTTTCCGGATGCGCGGCAAAGGCGCCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TGGCAGGGCCGCGTTGACGTCGGAACACCGGTCGACGTGAACG‑AGAACAAGAGAGAGCTGCTGCAAGAGATGCAAGAACGCTTCGGGGGGCCAACAGGAGAGCAAACCCGCACGCGC‑CACAAGACCATCTTTGATGGGGAGAAGACGATTTTTGACGACCt                                                                                                                       >  2:239963/1‑251 (MQ=255)                                                                                                                   ggCGAAACCCAGACCGGTAAGCTGTTCCGTATGCGCGGTAAAGGCGGCAAGTCTTACCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGGCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAAcc                                                                                                           <  1:259046/251‑1 (MQ=255)                                                                                                                               aCCGGTAAGCTTTTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGcc                                                                                               <  1:218698/251‑1 (MQ=255)                                                                                                                                       gcttttCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTgaa                                                                                                                                         >  1:110665/5‑201 (MQ=255)                                                                                                                                                 aTGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGcttg                                                                             >  1:11322/1‑248 (MQ=255)                                                                                                                                                   gcgcGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTggg                                                                           <  1:56903/251‑1 (MQ=255)                                                                                                                                                             ggcggcaggTCTGTCCGGGGTGGCGCTCAGGGTGAT‑TTGCTGTGCCGGGGTGTCGTCGAAACTCCGGTAGGCCCGAACG‑AAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGAAGACCTGACCCGCTAACCTCCCTAAAAGCCGGCCCGCGGGCAGGCCTGGGTAAAAATAgg                                                                 <  2:217816/243‑1 (MQ=255)                                                                                                                                                                   aaGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGCGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAAAAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGAAAAATTAGGGTgagt                                                           >  1:226369/1‑248 (MQ=255)                                                                                                                                                                       cTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGa                                                        >  1:314710/1‑250 (MQ=255)                                                                                                                                                                         atCCGCGGAGGCGCACAGGGAGAT‑TTGATGAGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑CGACGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCACCCGGTGAGCACAACAACCAGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAACAGATATTTGACGCCCAGACCCGCTAACCTACCCAAAAGCCGGCCCGTGGGAAGGCTTGGGTaaccttaaagcgcattgaggt                                                     >  2:263484/2‑232 (MQ=255)                                                                                                                                                                                 gggCGCACGGGGTGAT‑TTGGTAGGACGCGTTGTCGACAAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGTCCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTTCGTTGAAGATATGCGAg                                             <  2:166403/250‑1 (MQ=255)                                                                                                                                                                                    cgcACAGGGTGAT‑TTGCTGTTCCGCGTTGTCGTCGAAACACCGCTTGGCCCTAACG‑AGAAGCAGAAACAGCTGCTTCAAGCGCTGCAAGAAAGCTTCGGTGGCCCAATCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGACGGGGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGTGCAGGCTTGGGTAAAATTAGGGTGCATTGAAGATATGCGAGCAc                                          <  2:126426/251‑1 (MQ=255)                                                                                                                                                                                       acaGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTg                                       <  1:32788/251‑1 (MQ=255)                                                                                                                                                                                                         gcctcgTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACTCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCa                      <  1:141880/247‑1 (MQ=255)                                                                                                                                                                                                           cgagTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGCAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCTTGAGCAGGCTTGGGAAAAATTAGGGCACGTTGAAGATATTCGAGCACATGAAAAGGGGCggggtgacac                    >  2:21931/4‑245 (MQ=255)                                                                                                                                                                                                                   cgtcgAAACACCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCgggggagacacccctaa             >  1:382962/1‑237 (MQ=255)                                                                                                                                                                                                                   cgtcgAAACACCGGTAGGCCTGAACG‑AAAGGCAGAAACAGCCGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCTCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGTGATCACTCCCATAAg            <  1:190544/251‑1 (MQ=255)                                                                                                                                                                                                                            acCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCgggggacaatcccatcaacgcaaactt   >  1:351540/1‑228 (MQ=255)                                                                                                                                                                                                                            acCGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACtt   <  1:315028/251‑1 (MQ=255)                                                                                                                                                                                                                              cGGTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCAAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGGTGCGTTGAAGATATGCGAGCACCTGAAAAGTGGCGGGGaatacacccagacaaacgacatta  >  1:202854/1‑227 (MQ=255)                                                                                                                                                                                                                                gTAGGCCTGAACG‑AGAAGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGTGAGCACAACAGCACGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGACGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCTTGGGTAAAATTAGGATGCGTTGAAGATATGCGAGCAACTGAAAAGTGGCgggggatactactacaagagcaacaaga  >  1:335145/1‑224 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                         GAGCGTGAAGGCAACAACCTGTATTGCGAAGTCCCGATCAACTTCGCTATGGCGGCGCTGGGTGGCGAAATCGAAGTACCGACCCTTGATGGTCGCGTCAAACTGAAAGTGCCTGGCGAAACCCAGACCGGTAAGCTATTCCGTATGCGCGGTAAAGGCGTCAAGTCTGTCCGCGGTGGCGCACAGGGTGAT‑TTGCTGTGCCGCGTTGTCGTCGAAACACCGGTAGGCCTGAACG‑AAAGGCAGAAACAGCTGCTGCAAGAGCTGCAAGAAAGCTTCGGTGGCCCAACCGGCGAGCACAACAGCCCGCGCTCA‑AAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTA  >  NC_000913/14936‑15405

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