Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,955,010	G→A	D666D (GAC→GAT)	torZ ←	trimethylamine N‑oxide reductase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,955,010	0	G	A	100.0%	155.8 / NA	46	D666D (GAC→GAT)	torZ	trimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/25);  total (21/25)

CATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGA‑ACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAG  >  NC_000913/1954766‑1955248                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 cATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCACAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATcggcg                                                                                                                                                                                                                                            >  2:371182/1‑251 (MQ=255)     cACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTAc                                                                                                                                                                                                                                        >  1:479882/1‑251 (MQ=255)          aTACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTc                                                                                                                                                                                                                                   <  2:71550/251‑1 (MQ=255)            aCCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACAACGTTTGTTCCAGCCACGCACCACATCGCCATTCGCAATACCAAAGCGAGCAGCATCTTCGGCGTGAATAGCGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCCGATCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCCTTccc                                                                                                                                                                                                                                 >  1:84827/1‑250 (MQ=255)                     tttttGATCCCGTCAGTCACCCCCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAg                                                                                                                                                                                                                         <  2:565955/250‑1 (MQ=255)                     tttttGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAgg                                                                                                                                                                                                                        >  1:129030/1‑251 (MQ=255)                         tGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGAACACGTTTGTTCCAGACACGCACCAGATCGCCATTAGCAATACCAAAGCGTGCAGCCTCTTCGGTGTGAATAGTGATTGGTTAACGATCGGAGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGACGGCTGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGATCCTTTcccatcactaggagcc                                                                                                                                                                                                                    >  2:688898/1‑237 (MQ=255)                          gATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGagc‑c                                                                                                                                                                                                                   <  2:11538/250‑1 (MQ=255)                           aTCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGagc‑cag                                                                                                                                                                                                                 <  2:588130/251‑1 (MQ=255)                              ccGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGTTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCATGAgc‑cagcca                                                                                                                                                                                                              <  2:707786/251‑1 (MQ=255)                                   agtcaCCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTg                                                                                                                                                                                                         <  2:128300/251‑1 (MQ=255)                                    gtcaCCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCGTCGGTGTGAATAGTGATTGGTTCACGATCCGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTTTGTAAACGGTTTGCCGGATGTGCGGTCAGAAGCTGTAACGGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CATCCAGGTTgg                                                                                                                                                                                                        <  1:688898/251‑1 (MQ=255)                                        ccaccGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGGTGGgtgt                                                                                                                                                                                                    <  2:383042/251‑1 (MQ=255)                                                              cccccGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGAGCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCaa                                                                                                                                                                              <  2:192561/247‑1 (MQ=255)                                                                   gTTTGTTCCAGCCACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCATTCCCTTTTTCTGGAGC‑CCGCCCGGGTGGGTGTGCCGGGGCATCCTTATTGCCAAAattt                                                                                                                                                                         >  2:311249/1‑251 (MQ=255)                                                                            aGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CATCCAGGTTGGGGGTGCCGGGCCGTCCCTTTAACCACATTTTTCCCACGtt                                                                                                                                                                >  2:67138/1‑251 (MQ=255)                                                                                acGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTg                                                                                                                                                            <  2:557967/251‑1 (MQ=255)                                                                                  gCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGaa                                                                                                                                                          <  1:199925/251‑1 (MQ=255)                                                                                            gCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCa                                                                                                                                                 >  1:8951/1‑250 (MQ=255)                                                                                             ccATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTTCATTCATCAGGTGC‑CCGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAACATTTTCCAACGTTTTTGAATAGATTTCAAt                                                                                                                                               >  1:138637/1‑251 (MQ=255)                                                                                               aTTCGCAATACCAAAGCGTGCAGCATCGTCTGTGGGAATAGTTATTGGTTCCCCATCGGCGATCTCATATTTTTTACGCAGTTCCTCATAGTTAAGCTGACTGTGTAAACTGTGTGCCTGATGTGCGGTCAGAAGCTTTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGTTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATATATTGCAAttt                                                                                                                                             <  2:348327/251‑1 (MQ=255)                                                                                                  cGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCCTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGc                                                                                                                                          <  2:612658/251‑1 (MQ=255)                                                                                                    cAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGGGTGCCGGGCAATCCTTCTAGCCAAATTTTTCCACCGTTTTTGAATAGATTTCTATTTTTCCg                                                                                                                                        >  2:579841/1‑251 (MQ=255)                                                                                                    cAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGATCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCg                                                                                                                                        <  1:215336/251‑1 (MQ=255)                                                                                                        aCCAAAGCGTGCAGCATCGTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTg                                                                                                                                    <  2:347268/251‑1 (MQ=255)                                                                                                            aaGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTTCATTCATCAGGAGC‑CAGCCAGGTTGGGGGTGCCGGGCAATCCTTATAGCCCAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTTCCGCTTTGCGt                                                                                                                                >  2:409968/1‑251 (MQ=255)                                                                                                                  gcagcaTCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCACA‑CCAGGCTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGGTTTTGAATAGATTTTAATTTTTCCGCTTGGCGTaacc                                                                                                                            >  1:479496/1‑246 (MQ=255)                                                                                                                            cGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCAttttt                                                                                                                 <  2:479882/250‑1 (MQ=255)                                                                                                                                            tGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGG‑CAGCCAGGGTGGGGGTGCCGGGGAATCCTTATAGACAAATTTTTCTAAAGGTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGGTAGGCGCgg                                                                                                 >  1:356901/1‑250 (MQ=255)                                                                                                                                                   cGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTCACTGCTTCTCATCGGCGGTACCTTTTCATTACTCCGGAGG‑CCGGCACGTTGGGCGAGACCGGCAATCCTTCTAGCCAAATTTTTCCAACGATTTTGAATAGATTTCAATTTTTCCGCTTTGCGTACCCCACGCCTTTTTCACCCGCTCGGCGCGGAACTCAcc                                                                                         >  1:110263/1‑251 (MQ=255)                                                                                                                                                         gCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATaacg                                                                                   <  2:126338/251‑1 (MQ=255)                                                                                                                                                           gATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATaacga‑a                                                                                <  1:67138/251‑1 (MQ=255)                                                                                                                                                             tCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGGTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTTAATAGATTTCAATTTTTCCGCTTGGCGTACCCAGCGCCTTTTTCTCCGGATCGGCGCGGGAATCAACAtcaccaacc                                                                               >  1:369147/1‑243 (MQ=255)                                                                                                                                                               gCATATTTTTTACGCAGTGCCGCATAGTGAAGCTGACGGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCAGCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGTCAAATTGTTCCAACGTTGTTGAATAGATTGCAATGGTGCCGCTTGGCGTTCCCAGCGTATTTGTCACCGGATCGGCGCGGAATTCACCagatcttacgta                                                                             <  2:691255/251‑12 (MQ=255)                                                                                                                                                                cATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CCGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTTCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCCTAACGA‑ACGTAc                                                                           >  1:48599/1‑251 (MQ=255)                                                                                                                                                                       tttACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGA‑ACGTACTgttcg                                                                     <  1:409968/250‑1 (MQ=255)                                                                                                                                                                        ttACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGA‑ACGTACTgttcgt                                                                    <  2:422649/250‑1 (MQ=255)                                                                                                                                                                             cAGTTCCGCATAGTGAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTGTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGA‑ACGTACTGTTCGTtcttct                                                              <  2:447959/251‑1 (MQ=255)                                                                                                                                                                                  ccGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGGGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCCGCGCCTTTTTCACCGGATCGGCGCGGAAATCACCATAACGA‑ACGTACTGTTCGTTCTTCTCGCtg                                                         >  1:310406/1‑251 (MQ=255)                                                                                                                                                                                              gCTGACTGTGTAAACGGTGTGCCGAATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCCTTCATCAGGAAC‑CCGCCAGGTTGGGTGTGCCGGGCAATCCCTATAGCCAAATTTTTTCAACGTTTTTGAATAGATTTTAATTTTGCCGCTTTGCCGACCCAGCGCCTTTTTCACCCGATCGGCGCGGAAATCACCATAACGA‑ACGTACTTTTGCTTCTTCTCCCTGTGGCCCATTTCg                                             >  2:295728/1‑251 (MQ=255)                                                                                                                                                                                                  aCTGTGTATACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGAGGGACTTTGCCATTCATCAAGAGG‑CAGCCAGGTTGGGTGTACCGGGCCCTCCTTTTAGCCACATTTTTCCAACGGTTTTGAATTGGTTTCAATTTTATCGATCTGCGTACCCAGCGCATTTTTCCCCGGGTCCGCGCTGAAGTCACCCTAAAGC‑ACGTAATGTTCGTTCTTCTTGCCAGGGCGCCTTTCGATcc                                         >  2:328148/1‑250 (MQ=255)                                                                                                                                                                                                  aCTGTGTAAACGGGTTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTTCCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTCTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGA‑ACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCa                                         <  1:494496/251‑1 (MQ=255)                                                                                                                                                                                                       gtAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGA‑ACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGttt                                     <  2:129030/250‑1 (MQ=255)                                                                                                                                                                                                                   ccGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGCACCCTTTCATTCATCAGGGGC‑CCGCCAGGGTGGGTGTGCCGGGCAATCCTTTTAACCAAATTTTTCCAACGTTTTTGAATAGATTTCACTTTTTCCGCTTGGCGGACCCACCGGCTTTTTCACCGGGGCGGAGCGGGAAACACCCTAACTACACG‑ACTGTTTGTTCTTCTCGCAGTTGGGCCTTTCGCTCACTTTAtttttcttcccg                        >  2:238222/1‑243 (MQ=255)                                                                                                                                                                                                                     ggATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGAGCCGGGCAATCCTTATAACCAAATTTTTCCCAAGTTTTTGAATAGATTTTAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAAAGA‑ACGTACTGTTCGTTCTTCTCGCTGTGGCGCCTTTCGTTCAGTTTTTTTTTCTGcccga                       >  1:251519/1‑247 (MQ=255)                                                                                                                                                                                                                                          aaCTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGA‑ACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAg  <  2:70253/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 CATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGC‑CAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGCATTTTTCACCGGATCGGCGCGGAAATCACCATAACGA‑ACGTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAG  >  NC_000913/1954766‑1955248

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