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%	148.9 / NA	44	D666D (GAC→GAT)	torZ	trimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/23);  total (21/23)

CTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGC‑ATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTC‑TCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAGTG  >  NC_000913/1954763‑1955250                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   cTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCtcatcg                                                                                                                                                                                                                                                 >  2:177632/1‑251 (MQ=255)   tCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATcgg                                                                                                                                                                                                                                                <  1:254291/250‑1 (MQ=255)           gCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCtt                                                                                                                                                                                                                                       >  2:117673/1‑251 (MQ=255)                 cACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCAtt                                                                                                                                                                                                                                 <  2:64286/251‑1 (MQ=255)                  aCGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATtc                                                                                                                                                                                                                                <  1:29456/251‑1 (MQ=255)                                       gtcaCCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTgg                                                                                                                                                                                                           <  2:458319/251‑1 (MQ=255)                                                   ccTGTCAAAATCTGACCACGTTTGTTCCAGACCCGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCGTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCa                                                                                                                                                                                               <  2:403177/251‑1 (MQ=255)                                                    cTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTTCATTTATTACGAGCCAGCCAGGGTGGGTGTGCCGGGCa                                                                                                                                                                                               >  2:501826/1‑250 (MQ=255)                                                        cAAAATCTGACCACGTTTGTTCCAGACACGCACCATATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGGGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCt                                                                                                                                                                                          <  1:233804/251‑1 (MQ=255)                                                         aaaaTCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGTTTCACGATCGGCGATCGCATATTTTTTACGCAGTGCCGCATAGTTAAGCTGACTGTGTAATCGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCTGCGGTACCGTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCtt                                                                                                                                                                                         <  1:518341/251‑1 (MQ=255)                                                          aaaTCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCCGGTTGGGTGTGCCGGGCAATCCtt                                                                                                                                                                                         >  1:515477/1‑250 (MQ=255)                                                                gACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCa                                                                                                                                                                                  <  1:276424/251‑1 (MQ=255)                                                                  ccACGTTTGTTCCAGACACGCACCAGATCGCAATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCACTCCTTATAGCCaa                                                                                                                                                                                 >  1:137802/1‑250 (MQ=255)                                                                    aCGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAAtt                                                                                                                                                                              <  1:117673/251‑1 (MQ=255)                                                                          gTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCa                                                                                                                                                                        <  2:174020/251‑1 (MQ=255)                                                                              cAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGt                                                                                                                                                                    <  1:432801/251‑1 (MQ=255)                                                                                gACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGCACCTTTACATTCCTCAGGAGGCAGCCAGGTTGGGAGTGCCGGGCAATCCTTATAGCCCAATTTTTCCAACGttt                                                                                                                                                                  >  1:406969/1‑251 (MQ=255)                                                                                 acacGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCCGGTTGGGGGCGCCGGGCCATCCTTATAGCCAAATTTTTCCCACGtttt                                                                                                                                                                 >  2:466435/1‑251 (MQ=255)                                                                                 acacGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGtttt                                                                                                                                                                 <  1:501826/251‑1 (MQ=255)                                                                                 acacGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTAACTTTCCATTCATCAGGAGCCAGCCAGGTTGGGGGTGCCGGGCACTCCTTATTGCCCAATTTTTCCAACGtttt                                                                                                                                                                 >  2:405227/1‑251 (MQ=255)                                                                                            aTCGCCATTCGCAATACCAAAGCGTGCAGCATCTGCGGTGTGAATAGTGAATGGTTCACGATAGGCGAGCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTACCGGATGTGCGTTCAGAAGCTGTAACTGCTTATCATCGGCGGGAACTTTCGTGTCGTCAAGAGGCAGCCAGGGTGGGTGTGGCGGGCAATCCTTATATCGAAATATTTTCAACGTTTTTGAATAGAttt                                                                                                                                                      >  2:96606/1‑251 (MQ=255)                                                                                                  aTTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTTCTCAGGAGGCAGCCCGGGTGGGGGTGCCGGGCAATCCTTTTAGCCAAATTTTTCCAACGTTTTTGATTAGATTTTAAttt                                                                                                                                                >  2:240254/1‑251 (MQ=255)                                                                                                      gCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGcc                                                                                                                                            <  2:317666/251‑1 (MQ=255)                                                                                                       cAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCg                                                                                                                                           <  1:533076/251‑1 (MQ=255)                                                                                                             cAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTAACTTTTCATTCATCAGGAGACAGCCAGGTTGGGGGTGACGGGCAAACCTTATAACCAAAATTTTCCAAAGGTTTTGAATAGATTTCAATTTTGCCGCTTggg                                                                                                                                     >  1:59038/1‑250 (MQ=255)                                                                                                                aGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAAGCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTa                                                                                                                                  <  1:54006/251‑1 (MQ=255)                                                                                                                   gTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCCTTCATCCGGAGCCCACCAGGGTGGGGGTGGCGGGGCATCCTTTTTACACAATTTTTCCCACGTTTTTGGATAGATTTTATTTTTGACGCTTGGctgaccc                                                                                                                               >  2:307058/1‑245 (MQ=255)                                                                                                                                       aTAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGGACCTTTCCATTCCTACGGAGCCAGCCAGGTTGGGTGGGCCGGGCAATTCTTATCACCAAATTTTTCCAACGTTTTTGAATAGATTTTAATTTTTTCGCTTTGCGTACCCCACGC‑ATTTTTTACCTGATc                                                                                                          >  2:326145/1‑251 (MQ=255)                                                                                                                                         aGTGATTGGTTCGCGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTTAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGC‑ATTTTTCTCCGGATCg                                                                                                         <  2:554378/250‑1 (MQ=255)                                                                                                                                           tGATTGGTTCACGATCGTCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGTATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGCACCTTTTCAATTATGAGGGACCAGCCAGGGTGGGTGTGCCGGGGAATCCCTATTACAAAATTTTTCCCACGTTTTTGTATAGATTTTCATTATGCCGCGTGGCGTAccccgagagttttttctcggatcggtg                                                                                                       >  2:443763/1‑227 (MQ=255)                                                                                                                                                      cGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGGGGGCCGGGCAATCCTTATAGCCCAATTTTTCCCACGGTTTTGAATAGATTTCAATTTTTCCGCTTGGCGTACCCCGCGC‑ATTTTTCACCGGATCGGCGCGGGAATCAcc                                                                                           >  1:387042/1‑251 (MQ=255)                                                                                                                                                      cGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTAATCAGGAGCCAGCCAGGGTGGGGGGGCCGGGCAATACTTATAGACCAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTTGCGTACCCAGCGC‑ATTTTTCCCCGGATCGGCGCGGAAATCAcc                                                                                           >  1:551976/1‑251 (MQ=255)                                                                                                                                                              gATCGCATATTTTTTACGCAGTTCCGCATCGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCCGGTTGGGGGGGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGC‑ATTTTTCACCGGATCGGCGCGGAAATCACCATaacga                                                                                    >  1:209250/1‑250 (MQ=255)                                                                                                                                                               aTCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCCGGGTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTTCCGCTTGGCGTAACCAGCGC‑ATTTTTCACCGGATCGGCGCGGAAATCACCATaacgaa                                                                                   >  1:241661/1‑250 (MQ=255)                                                                                                                                                               aTCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGC‑ATTTTTCACCGGATCGGCGCGGAAATCACCATaacgaac                                                                                  <  1:447613/251‑1 (MQ=255)                                                                                                                                                                                  gTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGC‑ATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTC‑TCg                                                              <  2:141023/251‑1 (MQ=255)                                                                                                                                                                                        cATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCAGAGGTACCCTTCCCTTCATCACGGGCCCGCCAGGGTGGGTGTGCCCGGCAATCCTTATAACCAAATTTTTCTAACGTTTTTGAATAGATTTCAATTTTTCCCCTTTGAGTACCCAGCGC‑ATTTTTCACCGGATCGGCGCGGGAATCACCATCACCAACGCACTGTTCGTTCTTC‑TCGCTGTgg                                                        >  1:253311/1‑251 (MQ=255)                                                                                                                                                                                                                ggTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTCCCATTCCATTTATTAGGAGCCCCACAGGTTGGGGGGGTCGGGCAATCCTTATTACCAAAATTTTCCCACGGTTTTGAATCGATTTTCATTTTTCCGCCTGGCGTACCCCGCGCCATTTTT‑ACCGGATCAGCGCGGGAATCACCCTAACGCACGAACTTTTTGTTTTTCCTCCTTGTGGCG‑ATTTTGATCAGATTTTTTtgt                                >  2:298006/1‑250 (MQ=255)                                                                                                                                                                                                                 gtgtGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGC‑ATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTC‑TCGCTGTGGCGCATTTCGATCAGTTTATTTtgct                               <  1:177632/251‑1 (MQ=255)                                                                                                                                                                                                                      ccGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTTCCCAGCGC‑ATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTC‑TCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAg                          <  2:137802/251‑1 (MQ=255)                                                                                                                                                                                                                       cGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTTATCAAGAGCCAGCCAGGGTGGGTGTGGCGGGCAATCCTTATATCCAAATTTTTTCAACGTTTTTGAATAGATTTCAATTTTGACGCTTGGCCCACCCAGCGC‑ATTTTTTAACGGATCTGCGCGGAAATCACCATAACGAACGTACTGTTCGCTCTTC‑TCGCTGTGGCGCATTTCGTGCAGTTTATTTtgctgac                            >  1:403867/1‑246 (MQ=255)                                                                                                                                                                                                                                     gAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGC‑ATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTC‑TCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACAtt           <  1:466435/251‑1 (MQ=255)                                                                                                                                                                                                                                     gAAGCTGTAACTGCTTCTCATCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATACTTATAGCCAAATTTTTCCAACGGTTTTGAATAGATTTCAATTTTTCCGCTTGGCGTACCCAGCGC‑CTTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTTTTCGTTCTTC‑TCGCTGTGGCGCATTTCGGTCAGGTTATTTTTCTTCCAGAAGGCATTAAAAAtt           >  1:136986/1‑251 (MQ=255)                                                                                                                                                                                                                                     gAAGCTGTAACTGCTTCTCAGCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGTGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGC‑ATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTC‑TCGCTGTGGCGCATGTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACAtt           <  1:309229/251‑1 (MQ=255)                                                                                                                                                                                                                                              aCTGCTTCTCATCGGCGGTCCCTTACCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGC‑ATTTGTCACCGGATCGGCGCGGAAATCTCCATAACGAACGTACTGTTCGTTCTTC‑TCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCAGTAAACATTGGCATAGTg  <  2:304070/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   CTTCATGCACGCATACCACGCCTTTTTTGATCCCGTCAGTCACCACCGCGCCTGTCAAAATCTGACCACGTTTGTTCCAGACACGCACCAGATCGCCATTCGCAATACCAAAGCGTGCAGCATCTTCGGTGTGAATAGTGATTGGTTCACGATCGGCGATCGCATATTTTTTACGCAGTTCCGCATAGTTAAGCTGACTGTGTAAACGGTGTGCCGGATGTGCGGTCAGAAGCTGTAACTGCTTCTCGTCGGCGGTACCTTTCCATTCATCAGGAGCCAGCCAGGTTGGGTGTGCCGGGCAATCCTTATAGCCAAATTTTTCCAACGTTTTTGAATAGATTTCAATTTTGCCGCTTGGCGTACCCAGCGC‑ATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAACGTACTGTTCGTTCTTC‑TCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAGTG  >  NC_000913/1954763‑1955250

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