Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,196,232	C→T	T370T (ACC→ACT)	icd →	isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,196,232	0	C	T	100.0%	131.7 / NA	45	T370T (ACC→ACT)	icd	isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/31);  total (14/31)

GAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCA  >  NC_000913/1195942‑1196448                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                         gagaTCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGCAGCCACCCACGGTACTGCGCCGAAATATGCCGGTTAGGACAAAGTAAATCCTGGCCCTATTCTTCTCTCCGCTGAGATGATTCTGCGCCATATGGGTTGGACTGAAgggg                                                                                                                                                                                                                   >  2:837436/1‑298 (MQ=255)                  cctgagtgcgcaTGCATTCCTGCAACTGATCCTGCTGCGTCCGGCCTAATATCATGGTATGTCCTGTATGAACCTGAACGGTGACTACATCTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGt                                                                                                                                                                                               <  2:84058/290‑1 (MQ=255)                  cGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGt                                                                                                                                                                                               <  2:778840/301‑1 (MQ=255)                  cGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGACGTTATCGCCTGTATTAACCTGAACGGTGACTACATGTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCTCTGGTGCAAACATCGGTGACGAATTCGCCCTGTTTGATGCCACCCTCGGCACTGCGCCGAAATATCCCGGTCAGGTCAAAGTAAATCCTGGCTCTTTTGTTCTCTCCGCTGTGAGTTTGCTTCCCCAAATGGGTCGGACTGAAGCGGCTGAccgtacttggaaaggg                                                                                                                                                                                               >  1:456994/1‑286 (MQ=255)                        tGCTGATTCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCGGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGaa                                                                                                                                                                                         <  1:137669/301‑1 (MQ=255)                        tGCTGATGCATTCCTGCAACAGCTCCTGCTGCGTCCGGCTGAATATTATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGaa                                                                                                                                                                                         <  2:457007/301‑1 (MQ=255)                         gCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGCCCAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCCCCATATGGGGTGGACTTAAGCGGCTGCCCTGATTGGTAAAGGTATgggaa                                                                                                                                                                                        >  1:349742/1‑298 (MQ=255)                         gCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCTCCATATGGGGTGGACTGAAGCGGCTGACCTGATTGGTAAAGGTATGGaaa                                                                                                                                                                                        >  1:305948/1‑300 (MQ=255)                             aTGCATTCCTGCAACAGGTCCTGGTTCGTCCGGCTGAATATGATGTGATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcgc                                                                                                                                                                                    <  2:433962/301‑1 (MQ=255)                             aTGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGAATGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcg                                                                                                                                                                                     >  2:874912/1‑300 (MQ=255)                                     cTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaatg                                                                                                                                                                            >  2:396414/1‑299 (MQ=255)                                            tgagccgggggtgtcacgCTGAATATGAGGGTATGGCCTGTATGAACCTGCACGGTGAATACATTTCTGAAAGCC‑GGCACGGGCAGGTTGGGGGGAACGGTTTCGCCCCTGGTGCAAAAATCGGTGACGAATGCGCCCTGATTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAGATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAGAGGTATGGAAGGCGCAATCAATGCCAAGAc                                                                                                                                                                     <  1:270277/284‑1 (MQ=25)                                              aTCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTTCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCg                                                                                                                                                                   >  2:687194/1‑301 (MQ=255)                                                   gctgcGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAc                                                                                                                                                               <  2:57109/300‑1 (MQ=255)                                                    ttgcGTCCGGCTGAATATGATGTTTACGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCG‑GGGCAGGTTGGCGGTATCGGTATCGCCCCTGGTTCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAActt                                                                                                                                                             <  1:687168/300‑3 (MQ=38)                                                                       atgTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGGCTgat                                                                                                                                          >  2:103098/1‑301 (MQ=255)                                                                        tgtgctCGCCAGTATGACGCTGAACGGTTGTCACATGTCGGACGACCTGGCA‑GCGCAGGTTGGCGGGATCGGTATCGCCCCTGGTGCAAAAATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatg                                                                                                                                         <  1:291961/296‑1 (MQ=25)                                                                        tgTTTTAGCCCGTATGAACCTGAACGGTGACTACATTGCTTACGCCCGGGCA‑GCGCAGGGTGGCGGGATCGGGAGCGCCCCTGGTTCAAACATCGGTGACGAAGTCGCCCTGTTTGAAGACACCCCCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTAGTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatg                                                                                                                                         <  1:730673/301‑1 (MQ=25)                                                                         gTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGCTGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAACGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACGTATGACGTcagacttttgatgg                                                                                                                                        >  1:845976/1‑288 (MQ=38)                                                                          ttATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatgg                                                                                                                                        <  2:297141/300‑1 (MQ=255)                                                                          ttATCGCCTGTATGAACCTGAACGGTGACTACATTTCGGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatgg                                                                                                                                        <  2:196608/300‑1 (MQ=255)                                                                               gCCTGTATGAACCTGAACGGTAACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGCGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGCTGTTGCTGCGCCATATGGGTTGGTCTGCAGCGGCTGACCTGCTTGTTAAAGGTAGGGAACGCGCAATCAATGACACGCCCGTAACTTATGACTTCGCACGCCTGAGGGAAGgcg                                                                                                                                  >  2:99320/1‑301 (MQ=255)                                                                                   gTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTaa                                                                                                                              >  2:361161/1‑301 (MQ=255)                                                                                     tcaacactaaccggtggtcacgtttCGGACGCCCTGGCA‑GCGGAGGTTGGCGGTACTGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCGCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCtgagc                                                                                                                            <  1:78165/279‑5 (MQ=21)                                                                                     aTGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTaagc                                                                                                                            <  2:765321/301‑3 (MQ=37)                                                                                     aTGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTaagc                                                                                                                            <  2:445294/301‑3 (MQ=37)                                                                                          ccTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCATGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGATATGGAAGGCGCAATCAATGCCAAGAACGTAACTTATGAATTCGCAAGGCTGCTGGAAGGCGCTAAGccgcta                                                                                                                       >  1:743553/1‑300 (MQ=37)                                                                                                      ttccATGTCTGACGCGCTGGCA‑GAGCAGGTTGGCGGTATCGGGATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCTAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCagag                                                                                                           <  1:655987/298‑1 (MQ=21)                                                                                                                 aCGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGaag                                                                                                 <  2:477317/300‑3 (MQ=25)                                                                                                                                        cggtatcggtatCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtc                                                                           <  2:762990/300‑2 (MQ=21)                                                                                                                                                gtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTTTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtctctccatg                                                                   <  2:826949/300‑10 (MQ=18)                                                                                                                                                 tatcgCCCCTGGTGCAAACATCGGTGACGAATGCGACCTGTTTGAAGCCACCCACGGTACTGCGCCGAAAGATGCCGGGCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGAGGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtctctccatgtg                                                                 <  1:396403/301‑12 (MQ=17)                                                                                                                                                       cccTGGTGCAAACAGCCGGGACGAATCAGCCCGTGTTGAAGCCACCCCCGTTACTGCGCCGAAATATGGCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCGGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtctctccatgtgttaaat                                                           <  1:103093/301‑18 (MQ=14)                                                                                                                                                       cccTGGTGAAAACATAGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAATTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtctctccatgtgttaaat                                                           <  2:806189/301‑18 (MQ=17)                                                                                                                                                       cccTGGGGCAAAAATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtctctccatgtgttaaat                                                           <  2:846008/301‑18 (MQ=17)                                                                                                                                                                                            aaaGCCACCCACGGTAATGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCGCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttattt                      <  1:361151/300‑5 (MQ=11)                                                                                                                                                                                                       cGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatatt           >  2:265037/1‑286 (MQ=14)                                                                                                                                                                                                         ttACTTCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatattat         <  2:197130/300‑18 (MQ=14)                                                                                                                                                                                                          tACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTTTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatattat         <  2:349752/300‑18 (MQ=14)                                                                                                                                                                                                           tcTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGCCTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatattatta       <  1:265027/300‑20 (MQ=14)                                                                                                                                                                                                             tGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatattattaat     <  2:305958/301‑22 (MQ=14)                                                                                                                                                                                                              gcgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtc                                                                           >  2:630980/1‑229 (MQ=14)                                                                                                                                                                                                              gcgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAAtc                                                                           <  1:630960/230‑2 (MQ=14)                                                                                                                                                                                                              gcgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatattattaat     <  2:467450/300‑22 (MQ=14)                                                                                                                                                                                                                 ccGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTAATCTCTCCATGTGTTAAATATTGAAACGGGCGTATAACACGCCCGTTGttttatttatgtggatattattaatagc  <  2:901291/300‑25 (MQ=17)                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                         GAGATCGTCATTAAAGACGTGATTGCTGATGCATTCCTGCAACAGATCCTGCTGCGTCCGGCTGAATATGATGTTATCGCCTGTATGAACCTGAACGGTGACTACATTTCTGACGCCCTGGCA‑GCGCAGGTTGGCGGTATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCA  >  NC_000913/1195942‑1196448

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑