Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,790,821	A→T	I175I (ATT→ATA)	tdh ←	threonine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,790,821	0	A	T	100.0%	150.2 / NA	44	I175I (ATT→ATA)	tdh	threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (26/18);  total (26/18)

GCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAACTCCGCCATCACGTCATTGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACG  >  NC_000913/3790578‑3791048                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    gccgTGGTTCATGGTGTCGAGAATGGTACGAAACGCTGGCGCCGCACCGGACATTTCCAGACCGACATCAAAGCCTTCGGTCATGCCCAGTTCTGCCATCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAAAGTGTTTCGCCACCGCCGCAGCCATCATACCTATCggac                                                                                                                                                                                                                              >  2:227453/1‑249 (MQ=255)        ttCATGGTGTCGAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAGCCTTCGGTCATGCCCAGTTCTGCCATCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCAcc                                                                                                                                                                                                                       >  2:272148/1‑251 (MQ=255)                gtcgaGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAGCCTTCGGTCATGCCCAATTCTGCCATCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGAcc                                                                                                                                                                                                                >  1:294377/5‑250 (MQ=255)                              gAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAGCCTTCGGTCATGCCCAGTTCTGCCATCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGcc                                                                                                                                                                                                  >  1:409294/1‑250 (MQ=255)                                        cggcgCCCCGGACATTTCCAGACCGACATCAAAGCCTTCGGTCATGCCCAGTTCTGCCATCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATcga                                                                                                                                                                                       <  1:323839/251‑3 (MQ=255)                                        cggcgCACCGGACATTTCCAGACCGACATCAAAGCCTTCGGTCATGCCCAGTTCTGCCATCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATcg                                                                                                                                                                                        >  2:74472/1‑249 (MQ=255)                                                         ccAGACCGACATCAAAGCCTTCGGTCATGCCCTGTTCTGCCATCACATCAGTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCATTGCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGc                                                                                                                                                                      <  1:308876/251‑1 (MQ=255)                                                                      aaagccTTCGGTCATGCCCAGTTCTGCCATCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCg                                                                                                                                                          >  1:342062/5‑250 (MQ=255)                                                                         gccTTCGGTCATGCCCAGTTCTGCCATCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAgg                                                                                                                                                      >  1:221625/2‑251 (MQ=255)                                                                         cccTTCGGTCTTGCCCATTTCTGCCATCACATCATTAAGGTTTTCTTGAGCGACGTGGACCGCCCGGGTGATACCCATTTGGCGCGCCAGTTCGAGACGATATTCATTAACACCAGTGATCACCACATTGCGTGCCCCACCGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAgg                                                                                                                                                      <  1:240831/250‑1 (MQ=255)                                                                             tCGGTCATGCCCAGTTCTGCCATCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGCTCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCCTAATACCTATCGGACCTGCCCCAGAGCCCCGCACATCCTCGCCAACCAGATCGACCGGCAGCAAGGTATGCACGGCGTTACCGAACGGATc                                                                                                                                                  >  1:177353/1‑251 (MQ=255)                                                                                    tgcccagttctgCCATCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAAtg                                                                                                                                           >  2:72950/12‑250 (MQ=255)                                                                                     gcccagttctgCCATCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCTTTTTGCGCGCCATTTCGAGACGATATTCATTAACCTCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAAtgg                                                                                                                                          <  1:37800/241‑3 (MQ=255)                                                                                       ccagttctgCCATCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAAtggaa                                                                                                                                        <  1:23713/243‑5 (MQ=255)                                                                                                 catcacatcaTTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAAtggaagccaggtcgt                                                                                                                              >  2:396238/1‑237 (MQ=255)                                                                                                  aTCACATCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGAGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTACCATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCAGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGACCGACAGCGCGGTATGCACGGCGTTACCGAAGGGGTCAAAAAtggaagccaggtcgtc                                                                                                                             >  1:438507/1‑236 (MQ=255)                                                                                                       atCATTAAGGTTTTCTTTAGCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAAtggaagccaggtcgtcgg                                                                                                                           >  1:353733/2‑231 (MQ=255)                                                                                                                        tagCGACGTTGACCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGAtttt                                                                                                        <  1:332170/248‑1 (MQ=255)                                                                                                                                  gaCCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGa                                                                                             <  1:286917/250‑1 (MQ=255)                                                                                                                                   aCCGCACGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCCCCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGaa                                                                                            <  1:72950/251‑1 (MQ=255)                                                                                                                                        aCGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCg                                                                                       >  1:316260/1‑251 (MQ=255)                                                                                                                                        aCGGGTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCg                                                                                       >  1:237303/1‑251 (MQ=255)                                                                                                                                          gggTGATACCCATTTTGCGCGCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCggg                                                                                     >  1:329362/1‑251 (MQ=255)                                                                                                                                            gggacaCCCATTTTGCGCGCCAGTTCGATACGATATTCATTAACATCAGGGACCACCACATTCCGTGCCCCAACGGGTTCCGCCTCCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGAt                                                                                   <  2:309999/246‑1 (MQ=255)                                                                                                                                                         tgcgcgCCAGTTCGAGACGATATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCg                                                                      >  1:290646/3‑251 (MQ=255)                                                                                                                                                         tgcccgccagttcgagacgttATTCATTAACATCAGTGATCACCCCATTGCGTGCCCCAACGTGTTTCGCCCCCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTc                                                                       <  2:147274/230‑1 (MQ=255)                                                                                                                                                             cgccagttcgagacgatATTCATTAACATCCGTGCTCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGTATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGcg                                                                   >  1:141177/17‑249 (MQ=255)                                                                                                                                                                     cgagacgatATTCATTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGccc                                                          >  1:96581/9‑251 (MQ=255)                                                                                                                                                                     cgagacgatATTCATTAACATCAGGGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGcc                                                           <  1:272148/242‑1 (MQ=255)                                                                                                                                                                                 cattAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAacg                                              >  2:318795/3‑250 (MQ=255)                                                                                                                                                                                     accaacAGTGATCACCACATTGCGTCCGTCATCGTGTTTCGCCTCCGCCGCAGCCATAATACCTATCGGATCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAacgccga                                          <  1:151002/246‑6 (MQ=255)                                                                                                                                                                                           aGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCACGGACGGTTAACGCCGACGccg                                     >  1:332659/1‑249 (MQ=255)                                                                                                                                                                                            gTGATCACCACAGTGCGTGCACCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGccgat                                   <  1:396238/251‑4 (MQ=255)                                                                                                                                                                                               aTCACCACATTGCGTGCCCCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCTGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCgt                                <  1:74472/251‑1 (MQ=255)                                                                                                                                                                                               aTCACCACATTGCGTGCACCCAAGTGATTCGCAACCGCCGCAGCAATCATACCTACCGGACCTGCACCAGAGCCCAGCACATCCTCGCCACCCAGCTCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGCTCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCgt                                >  2:383074/1‑251 (MQ=255)                                                                                                                                                                                                                ccaaaCGGGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCGTGTTGCGGCACAGATggg               <  2:187316/248‑1 (MQ=255)                                                                                                                                                                                                                aCCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATc                                  <  1:45921/232‑1 (MQ=255)                                                                                                                                                                                                                aCCAACGTGTTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATc                                  >  2:45921/1‑232 (MQ=255)                                                                                                                                                                                                                    aCGTGTTCCCCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCGTGTTGCGGCACAGATGGGTACg           <  2:96581/251‑1 (MQ=255)                                                                                                                                                                                                                      gtgtTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAagca                                                              <  1:79029/198‑1 (MQ=255)                                                                                                                                                                                                                      gtgtTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAagca                                                              >  2:79029/1‑198 (MQ=255)                                                                                                                                                                                                                      gtgtTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCGTGTTGCGGCACAGATGGGTACGac         >  2:121527/1‑251 (MQ=255)                                                                                                                                                                                                                      gtgtTTCGCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCGTGTTGCGGCACAGATGGGTACGac         >  1:226537/1‑251 (MQ=255)                                                                                                                                                                                                                             gCCACCGCCGCAGCCATAATACCTATCGGACCTGCACCAGAGACCAGCACATCCTCGCCAACCAGATCGAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATGGAAGCCAGGTCGTCGGAGATATTGTCTGGGATTTTGAAGGCGTTGAACGCCGGGATCACCAGATATTCGGCGAAGCAGCCCGGACGGTTAACGCCGACGCCGATCGTGTTGCGGCACAGATGGGGACGACCGCCACg  >  1:73838/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                    GCCGTGATTCATGGTGTCAAGCATGGTACGAAACGCTGGCGGCGCACCGGACATTTCCAGACCGACATCAAAACCTTCGGTCATGCCTAACTCCGCCATCACGTCATTGAGATTTTCTTTGGCGACGTTAACCGCACGGGTGATACCCATTTTACGCGCCAGCTCAAGGCGGTATTCGTTAACATCAGTGATCACCACATTGCGTGCACCAACGTGTTTCGCCACCGCCGCTGCCATAATACCAATCGGGCCTGCACCAGAAACCAGCACATCTTCGCCCACCAGATCAAACGACAGCGCGGTATGCACGGCGTTACCGAAGGGATCAAAAATTGCGGCTAAGTCATCGGAAATATTGTCGGGGATTTTGAAGGCGTTGAATGCCGGGATCACCAGATATTCGGCAAAGCAGCCCGGGCGATTAACACCAACGCCTATCGTGTTGCGGCACAAATGGGTACGACCACCACG  >  NC_000913/3790578‑3791048

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: ‑