Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,073,610	G→A	N134N (AAC→AAT)	rutA ←	pyrimidine monooxygenase RutA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,073,610	0	G	A	100.0%	164.7 / NA	48	N134N (AAC→AAT)	rutA	pyrimidine monooxygenase RutA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/26);  total (22/26)

GTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGG  >  NC_000913/1073315‑1073894                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                             gtgtATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGa                                                                                                                                                                                                                                                                                          >  2:611446/1‑300 (MQ=255)     aTTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGcc                                                                                                                                                                                                                                                                                      >  1:225647/1‑300 (MQ=255)     aTTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGcc                                                                                                                                                                                                                                                                                      >  2:609246/1‑300 (MQ=255)                    gAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTTAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCaga                                                                                                                                                                                                                                                                       <  1:379821/300‑1 (MQ=255)                                  aaTCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCAAAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATgg                                                                                                                                                                                                                                                         <  1:488088/300‑1 (MQ=255)                                         aTACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGc                                                                                                                                                                                                                                                  >  1:439875/1‑300 (MQ=255)                                           aCCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCa                                                                                                                                                                                                                                                <  1:629839/300‑1 (MQ=255)                                              gggCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATAc                                                                                                                                                                                                                                             <  2:474701/300‑1 (MQ=255)                                               ggCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATAc                                                                                                                                                                                                                                             >  1:260957/1‑299 (MQ=255)                                                                            gCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAgg                                                                                                                                                                                                               <  2:599079/300‑1 (MQ=255)                                                                            gCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCAGCGCCAGGCCAGATACCCATCTGCTCATACTCGCGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAgg                                                                                                                                                                                                               <  2:622850/300‑1 (MQ=255)                                                                            gCCCGGAGCAGATCACTTTCATGGGGACACTCGTTTGCGTACTCACGCGACAATCATTCATGTTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACGGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCAGCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAgg                                                                                                                                                                                                               <  2:622852/300‑1 (MQ=255)                                                                                   gcAGATCACTTTCATGTNGACACTCGGTTGCGGACTNACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTNACGTTAAGGTGgcagc                                                                                                                                                                                                        <  1:602517/300‑1 (MQ=255)                                                                                     aGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGTTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACACCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGcagcag                                                                                                                                                                                                      <  2:77574/300‑1 (MQ=255)                                                                                                attgggACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTTCCCCACAGGTCGCGCAGCACCTGAACATATTCGTTTAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGTCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATc                                                                                                                                                                                           <  1:383563/297‑1 (MQ=255)                                                                                                      acacTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAGATCGCCTGTAAAATCGCCTTTCCCCGTGCCCCACAGGGCGCGCAGCACCTGAACATATGCGTTGAAATAGTCGTAACGACGGGAGAAATATGCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACACCCGCCATAGATGGAGTCAATGGTTGCGGCCATACGGGCGACGATTGCCGTAGGTAACGtcaatcc                                                                                                                                                                                                              <  2:398606/275‑7 (MQ=255)                                                                                                                       cACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGgccgc                                                                                                                                                                    >  1:264106/1‑300 (MQ=255)                                                                                                                         cgcgACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCa                                                                                                                                                                  >  1:620272/1‑300 (MQ=255)                                                                                                                             aCAATCATTCATGGTGAAAAAATCACCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGACAGAGGTGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCAAGGCCGCCAGcc                                                                                                                                                               >  2:368229/1‑299 (MQ=255)                                                                                                                                  cattcatGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGACTGCGCGAGGTCACGGCCGCCAGCCCCGcac                                                                                                                                                         >  2:93794/1‑298 (MQ=255)                                                                                                                                                        tttAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAATATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGTCCAGATACCCATCTGCTCATACTCGTGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGATTCAAGGt                                                                                                                                    <  2:647003/299‑1 (MQ=255)                                                                                                                                                        tttAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGt                                                                                                                                    <  1:293243/299‑1 (MQ=255)                                                                                                                                                        tttAAAATCGCGTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCTGTTAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGTCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCATAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGTAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGt                                                                                                                                    <  2:385588/299‑1 (MQ=255)                                                                                                                                                                           gTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCAGCAAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGtt                                                                                                                >  2:301535/1‑300 (MQ=255)                                                                                                                                                                            tGCCCCACAGGTCGCGCAGCACCTGAACATATTCGTTGAGATAGTCGTAACGACGGTAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCNAAACGCCCGCCATAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCACCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGttt                                                                                                               <  2:229693/300‑1 (MQ=255)                                                                                                                                                                                             agcaCCTGAACATATTCGGTGAGATAGTCGTAACGACGGTAGAAATAGTCATCGCCATGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCATTCACGAGATTGACGCCAAAACGCCCGCCATAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGc                                                                                              <  2:338987/300‑1 (MQ=255)                                                                                                                                                                                              gcaCCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCa                                                                                             >  1:111029/1‑300 (MQ=255)                                                                                                                                                                                              gcaCCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCAAGCa                                                                                             >  1:111079/1‑300 (MQ=255)                                                                                                                                                                                               caCCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCAATCAAGGTGAACGACTCAAGGTTGTGATCCAAGAACTCAGTTTTGCAGCCAAAGCCACGCa                                                                                             >  2:263955/1‑299 (MQ=255)                                                                                                                                                                                                 ccTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGtt                                                                                          >  2:652598/1‑300 (MQ=255)                                                                                                                                                                                                            tCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGAc                                                                               <  1:607120/300‑1 (MQ=255)                                                                                                                                                                                                                             acgacgGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCATAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGTAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGcgaaatcgaaagg                                                              <  1:363041/300‑3 (MQ=255)                                                                                                                                                                                                                                                 cGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCANACGGGCGACGATTGCCGGAGGTAACGTTAAGGNGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGca                                          >  2:529032/1‑300 (MQ=255)                                                                                                                                                                                                                                                       gccagATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAAt                                      >  1:165958/1‑298 (MQ=255)                                                                                                                                                                                                                                                          agATACACATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCAAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCtt                                 >  2:394950/1‑300 (MQ=255)                                                                                                                                                                                                                                                            aTACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCACCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCAAGAACTCAGTTTTGCCGCCAAAGCCAAGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTa                               >  1:529304/1‑300 (MQ=255)                                                                                                                                                                                                                                                                  aTCTGCTCATACTCGTGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCATCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTc                            <  2:600933/297‑1 (MQ=255)                                                                                                                                                                                                                                                                  aTCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCATCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTc                            <  2:601160/297‑1 (MQ=255)                                                                                                                                                                                                                                                                  aTCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTc                            <  1:484379/297‑1 (MQ=255)                                                                                                                                                                                                                                                                    cTGCTCATACTCGGGTTTTTGCAAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCAAGGCCGCCAGCCACNCCATCAAGGTGAACGACTCAAGGTTGTGATCCAAGAACTCCGTTTTGCCGCAAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGACATGGTCGTGCTCCGCTTTTTGCACAATGGCTTTagtc                            >  2:233859/1‑292 (MQ=255)                                                                                                                                                                                                                                                                      gCTCATACTCGNGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAgg                                                                                                                                                                                                               <  2:545613/114‑1 (MQ=255)                                                                                                                                                                                                                                                                      gCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAgg                                                                                                                                                                                                               >  1:545613/1‑114 (MQ=255)                                                                                                                                                                                                                                                                      gCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCaa                     >  1:406306/1‑300 (MQ=255)                                                                                                                                                                                                                                                                          aTACTCTGGCTTTTGCCAGCCAGTCACGAGATGGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAgg                  <  1:609237/299‑1 (MQ=255)                                                                                                                                                                                                                                                                          aTACTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAgg                  <  1:609246/299‑1 (MQ=255)                                                                                                                                                                                                                                                                           tccTCGGGCTTTTGCCAGCCAGTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCAGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTc                >  1:20501/3‑300 (MQ=255)                                                                                                                                                                                                                                                                                           agccagTCACGAGATTGACGCCAAAACGCCCGCCATAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTgcgg  <  1:549898/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                           agccagTCACGAGATTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTgcgg  <  1:549166/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                             GTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCACTTTCATGGGGACACTCGGTTGCGGACTCACGCGACAATCATTCATGGTGAAAAAATCGCCTTTAAAATCGCTTTTCCCCGTGCCCCACAGGTCGCGCAGCACCTGAACATATTCGGTGAGATAGTCGTAACGACGGGAGAAATAGTCATCGCCAGGCCAGATACCCATCTGCTCATACTCGGGCTTTTGCCAGCCAGTCACGAGGTTGACGCCAAAACGCCCGCCAGAGATGGAGTCGATGGTTGCGGCCATACGGGCGACGATTGCCGGAGGTAACGTTAAGGTGGCAGCAGTGGCGTAAATCTGAATGCGCGAGGTCACGGCCGCCAGCCCCGCCATCAAGGTGAACGACTCAAGGTTGTGATCCCAGAACTCAGTTTTGCCGCCAAAGCCACGCAGTTTGATCATCGACAGGGCGAAATCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGG  >  NC_000913/1073315‑1073894

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