Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	674,259	G→A	C175C (TGC→TGT)	leuS ←	leucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	674,259	0	G	A	95.0%	170.1 / 0.5	61	C175C (TGC→TGT)	leuS	leucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (29/28);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

AGGTACAACCCATAAAGGTGTCCGGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGA‑TCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAGCAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGT  >  NC_000913/674012‑674503                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     aggtaCAACCCATAAAGGTGTCCGGGCGGGTAGTGTAACCGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGA‑TCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCcagcaacag                                                                                                                                                                                                                                                   <  2:201981/251‑1 (MQ=255)         cccATAAAGGTGTCCGGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGA‑TCTCCACGCCTTCGGAACGACCGATCCATTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATACAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGAGCCAGCACCAGCCGtaca                                                                                                                                                                                                                                            >  2:605874/1‑247 (MQ=255)          ccATAAAGGTGTCCGGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGCAGGTGA‑TCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGAt                                                                                                                                                                                                                                           >  1:297184/1‑250 (MQ=255)          ccATAAAGGTGTCCGGGCGGGAAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGA‑TCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGAt                                                                                                                                                                                                                                           >  2:441859/1‑250 (MQ=255)                aGGTGTCCGGCCGGGGAGTGTAAACGGTCAGCGTTTTTTCATATTCGTGTACGTTGAAGGTGA‑TCTCCACGCCGTCTGAACGGCCTATCCAGTTACGCTGCCTGGTTTTAAGGGTGTCTGGCCAGTGATTCAGTTTATCCAGATCGTTGAGCAGCTCGTTAGCGTAAGCAGGGATTTTGATAATCCACTTCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCTTTAGCAAGAGCCGTCGATAACgtgt                                                                                                                                                                                                                                    <  2:583032/251‑1 (MQ=255)                  gtgtCCGGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGA‑TCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCCGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTgttc                                                                                                                                                                                                                                  >  1:654974/1‑251 (MQ=255)                  gtgtCCGGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAATCGTTAACGTTGAAGGTGA‑TCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGCTCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGACCGCAACAGCCGTCGATAACTTgttc                                                                                                                                                                                                                                  >  1:195756/1‑251 (MQ=255)                                   gtAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGA‑TCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCt                                                                                                                                                                                                                 >  1:551964/1‑251 (MQ=255)                                      aaCGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGGGA‑TCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTggg                                                                                                                                                                                                              >  1:520210/1‑250 (MQ=255)                                                gtgtTGTCATAGTCGTTAACGTTGAGGGTGA‑TCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGTTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCa                                                                                                                                                                                                    <  1:431991/251‑1 (MQ=255)                                                     gtCATAGTCGTTAACGTTGAAGGTGA‑TCTCCAGGCCTTCGGAACGGCCGATCCAGGTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTGTATCCAGATCGTTGGGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCGCTGCGGGGTCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGt                                                                                                                                                                                               <  1:614612/251‑1 (MQ=255)                                                      tCATAGTCGTTAACGTTTAAGGTGA‑TCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAGCAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGtt                                                                                                                                                                                              >  1:379259/1‑251 (MQ=255)                                                           tgctctcaCGTTGAAGGTGAGTCACC‑CGACTTCGGACCGCCCGATCCAGTTACGTGGAATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATTCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGAGCTCTTTACGTCCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTGCGGGCACCAGTTGAcc                                                                                                                                                                                          <  2:520210/243‑1 (MQ=255)                                                                                      cGCCTTCGGAACGACCGTTCCATTTACGCTGGATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGGCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTTCGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCt                                                                                                                                                               <  2:643925/251‑1 (MQ=255)                                                                                             ggAACGACCGATCCAGTTACGCTGCGTGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTtat                                                                                                                                                        >  1:652914/1‑251 (MQ=255)                                                                                             ggAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTtat                                                                                                                                                        <  1:267664/251‑1 (MQ=255)                                                                                             ggAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGAGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCCGTGATTTTGATAAAACACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCCCCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTtat                                                                                                                                                        >  1:636276/1‑251 (MQ=255)                                                                                              gAACGTCCGATCCAGTTACGCTGCATGGGTTTAACGGTGTCTGGCCAGTGAGCCAGTTTATCCAGATCGTTGGGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTTCGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACTGTCTGGTCGTTCGGTCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCGTTTTtata                                                                                                                                                       <  2:297184/251‑1 (MQ=255)                                                                                              gAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTtata                                                                                                                                                       <  1:436367/251‑1 (MQ=255)                                                                                               aaCGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATAc                                                                                                                                                      >  1:179795/1‑251 (MQ=255)                                                                                                  gACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGc                                                                                                                                                   <  1:294238/251‑1 (MQ=255)                                                                                                        tCCAGTTACGCTGCATGGTTTTAACGGTGTCTAGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTACGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGAAACAGCAGTAGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGCCCGCAGAAGTCTTCTTATATACAAGCCCTTTTTTATACAGCTCGGTg                                                                                                                                             >  2:680877/1‑251 (MQ=255)                                                                                                         ccAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTg                                                                                                                                             <  1:384574/250‑1 (MQ=255)                                                                                                             ttACGCTGCATGCTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTgaaga                                                                                                                                         >  1:158278/1‑250 (MQ=255)                                                                                                                 gCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTc                                                                                                                                    >  1:542999/1‑251 (MQ=255)                                                                                                                        ggTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTcc                                                                                                                              <  1:458781/250‑1 (MQ=255)                                                                                                                                cGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGgt                                                                                                                       >  1:664352/1‑249 (MQ=255)                                                                                                                                 ggggtcgtgccaatgactcAGTTTATCCAGATCGTTGTGCATCTCGTAAGCGTAAGCAGGGATTTTGATAAACCACTGCTGGATCATTTTACGTTCATCTTTGGTTGCGCAGCGCCAGCAACAGCCGTCTATAACTTGTTCGTTCGCCAGTACGGTCTGGTTGTTCGGGCTCCAGTTGGCCGCAGAAGTCTTCTTATATACCAGTCCTTTTTTATACAGCGCGGTGAAGAATTTCTGTTCCCAGCGgtagt                                                                                                                    <  2:610424/233‑1 (MQ=255)                                                                                                                                   tgtCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTAt                                                                                                                  >  2:46017/1‑251 (MQ=255)                                                                                                                                    gttcggCCAGTGGTCCAGTTTATCCAGATCGTGGGGCAGCTCGTCAGCGTAAGCTGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCACCTTTGGTATCGCAGCGCCAGCAACAGCCGGCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGTCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTAtt                                                                                                                 <  2:24096/247‑1 (MQ=255)                                                                                                                                      cgtgCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTcc                                                                                                               <  1:466679/248‑1 (MQ=255)                                                                                                                                      cTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTcc                                                                                                               >  1:590946/1‑251 (MQ=255)                                                                                                                                         gctaGTGATCCAGTTTTTCCAGATCGTTGAGCAGCTCGTCAGCTTAAGCAGTTATTTTGATAAACCATTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCgg                                                                                                             <  2:615462/247‑1 (MQ=255)                                                                                                                                           cAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGt                                                                                                          >  1:166103/1‑251 (MQ=255)                                                                                                                                                                  gTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAATCCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTcc                                                                                   <  2:295080/251‑1 (MQ=255)                                                                                                                                                                   ttGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTcc                                                                                   >  2:122408/1‑250 (MQ=255)                                                                                                                                                                     gAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGt                                                                                >  1:148149/1‑251 (MQ=255)                                                                                                                                                                      agcagcTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTc                                                                               >  2:394339/1‑251 (MQ=255)                                                                                                                                                                        caccTCGTCAGCATAGACAGTGGTTTTGATAAACAACTGCGGGATCTCTTTACGTTAAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGATCGCGCACCAGTTGACCGCAGAAGTCTTCTTATATACCCTGCCTTGTTTATACAGCGCGGTGCAGACTTTCTGTTCCCAACTGTAGTCTTCCGGCGAACAGGTTGCTAGATCGCGGCTCCAGTCAt                                                                             >  2:245598/4‑251 (MQ=255)                                                                                                                                                                             gttctgCGTAAGCAATGATTTTGACAAACCATTGCAGGATCTCTTTACGCACAAATTTGGTACTGAAGCGCCCGCAACAGCCGTCGCTAACTTGTTCGTTCGCAAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCACCTTATATACCAGGCCTTTTTTATCCAGCTCGGTGAACAATTTCTGTTCCCAACGGTATTATTCCGGCATACAGGTTGCCAGATCGCGGCTCCAGTCCTAACCa                                                                        <  2:653327/246‑1 (MQ=255)                                                                                                                                                                                 aGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGc                                                                    >  1:234076/1‑251 (MQ=255)                                                                                                                                                                                           gTGATTTTGATAAACCACGGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCAttt                                                           <  2:590946/250‑1 (MQ=255)                                                                                                                                                                                            tGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCACAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTCCAGGTTGCCAGCTCGCGGCTCCAGTCATCACCAAAGCCCCGCATTTTg                                                         >  2:170860/1‑251 (MQ=255)                                                                                                                                                                                             gATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCACAGCCCAGCATTTTgt                                                        >  1:660811/1‑250 (MQ=255)                                                                                                                                                                                                 ttGATAAACCACTGCGGGATCTCTTTACGTTCACCTTGGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCt                                                     <  2:306374/250‑1 (MQ=255)                                                                                                                                                                                                     tAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGttt                                                >  1:147262/1‑251 (MQ=255)                                                                                                                                                                                                      aaaCCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTTTTCCCAACgg                                                                                                                        <  2:289225/178‑1 (MQ=255)                                                                                                                                                                                                      aaaCCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACgg                                                                                                                        >  1:289225/1‑178 (MQ=255)                                                                                                                                                                                                      aaaCCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAACTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGtttt                                               <  1:46017/251‑1 (MQ=255)                                                                                                                                                                                                        aCCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTc                                             <  1:680877/251‑1 (MQ=255)                                                                                                                                                                                                        aCCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTc                                             >  2:520427/1‑251 (MQ=255)                                                                                                                                                                                                        aCCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCACAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTc                                             >  1:664108/1‑251 (MQ=255)                                                                                                                                                                                                         ccACTGCGGGATTTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCa                                            <  1:87003/251‑1 (MQ=255)                                                                                                                                                                                                              gCGGGTTCTTTTTACTTTCAACTTGGGTATGGCAGTGCCAGCAACGGCCGTCGATAACTTGTTCGTTCGCCAGTTCGGTTTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCAtata                                        <  2:474225/250‑1 (MQ=255)                                                                                                                                                                                                                           cctttCAACTTTGGTATCACAGCGCCTGCAACAGCCGTCGTTATCTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCTGGCACCCGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGACGAATTTCTGTGCCCAACGGTAGTATTCCGGCGTACAGGTTGCTATCGCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACTCGATGTTGTCg                           <  2:626582/247‑1 (MQ=255)                                                                                                                                                                                                                           aCGTTCAACTTTGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCg                           <  1:441859/250‑1 (MQ=255)                                                                                                                                                                                                                                 aaCTTTGGTATCGCAGCGCCAGCAGCAGCCGTCTATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCt                                                     <  2:104479/218‑1 (MQ=255)                                                                                                                                                                                                                                 aaCTTTGGTATCGCAGCGCCAGCAGCAGCCGTCTATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCt                                                     >  1:104479/1‑218 (MQ=255)                                                                                                                                                                                                                                     ttGGTATCGCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACg                 <  1:260219/250‑1 (MQ=255)                                                                                                                                                                                                                                             gCAGCGCCAGCAACAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCgg          >  1:114021/1‑249 (MQ=255)                                                                                                                                                                                                                                               aGCGCCAGCAACAGCCGTGGATATCTTTTTCGTTCGCCATTACGGTCTGGTCGGTCGGGCCCCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCCTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCTCGGTGCCGGAg        <  2:5304/249‑1 (MQ=255)                                                                                                                                                                                                                                                   ccAGCAACAGCAGTCGATACCTTGTTCGTTCGCCAGTACGGTCTGGCCGTTCGGGCCCCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGCATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGAGGCGgtgt  >  2:395949/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     AGGTACAACCCATAAAGGTGTCCGGGCGGGTAGTGTAAACGGTCAGCGTGTTGTCATAGTCGTTAACGTTGAAGGTGA‑TCTCCACGCCTTCGGAACGACCGATCCAGTTACGCTGCATGGTTTTAACGGTGTCTGGCCAGTGATCCAGTTTATCCAGATCGTTGAGCAGCTCGTCAGCGTAAGCAGTGATTTTGATAAACCACTGCGGGATCTCTTTACGTTCAACTTTGGTATCGCAGCGCCAGCAGCAGCCGTCGATAACTTGTTCGTTCGCCAGTACGGTCTGGTCGTTCGGGCACCAGTTGACCGCAGAAGTCTTCTTATATACCAGGCCTTTTTTATACAGCTCGGTGAAGAATTTCTGTTCCCAACGGTAGTATTCCGGCGTACAGGTTGCCAGCTCGCGGCTCCAGTCATAACCAAAGCCCAGCATTTTGAGCTGGTTTTTCATATACGCGATGTTGTCGTACGTCCACGGTGCCGGAGCGGTGT  >  NC_000913/674012‑674503

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