Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,652,331	T→C	intergenic (+2346/+397)	ydfD → / ← ynfP	lysis protein/protein YnfP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,652,331	0	T	C	100.0%	152.2 / NA	42	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/17);  total (25/17)

GGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGC  >  NC_000913/1652089‑1652575                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     ggTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCagaga                                                                                                                                                                                                                                               <  1:378389/250‑1 (MQ=255)          cTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCCCTCCGGAGTCAGATGAAATTAAAGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCCGAGAGTATGGc                                                                                                                                                                                                                                        >  1:361933/1‑248 (MQ=255)                     ggAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaaa                                                                                                                                                                                                                            <  1:403075/249‑1 (MQ=255)                          tCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGGCATTCATCTT‑TTTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGGGACTGAGCAAA‑AGACATATCCTTGATTTTCACCTCAGAGAGTATGGCAGAACAGAaaaaaaaaaa                                                                                                                                                                                                                      >  2:144667/1‑250 (MQ=255)                           cGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATg                                                                                                                                                                                                                    <  1:153081/251‑1 (MQ=255)                                             gCACTTGGATGGGAGGACATAGATCTGAAAAATGTAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAcc                                                                                                                                                                                                  <  2:301826/251‑1 (MQ=255)                                                      tGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAAAAGAAAAACAAAAATGCACCTTTGTTTTTTAACCTGAAGTGTc                                                                                                                                                                                         >  1:337499/1‑251 (MQ=255)                                                               aTAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGGCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAG‑CATATCATTGATGTTTACCTCCGAGAGTATGGCAGAAAAGAAAAACAAAAATGCACCTTTGTTTTTCAAACTGAAGTGTCAGCGggagg                                                                                                                                                                                >  1:366571/1‑250 (MQ=255)                                                                tAGATCTGAAAAATGGAACAATGATGATCCGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGAAGCTCTCCGGAGTCAGATGACATTAAAGAGACTGAGCAAAGAG‑CATATCATTTATGTTCACCTCAGAGAGTATGGGAGAAAAGAAAAACAAAAATGCCCCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGt                                                                                                                                                                                >  2:176092/1‑250 (MQ=255)                                                                 aGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCGCCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTa                                                                                                                                                                               <  1:32277/250‑1 (MQ=255)                                                                           aaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAA‑CCCGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAG‑CATATCATTGATGTTAACCTCAGAGAGTATGGCAGAAAAGAAAAAAAAAAAAGCACCTTTGTTTTTAAACCTGAAGTGTCAGCGAGAGGAAAAAATTATgg                                                                                                                                                                    >  2:302738/1‑251 (MQ=255)                                                                              ggAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATACGGTCAATCATATT‑ATTAAA‑CAAGCAAAAGACACTATCCGGAGTCAGATGACATTAACGAGAAAGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGCAAAGGAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCATCGGGAGTAAAAAATTATgggt                                                                                                                                                                  >  1:16190/1‑248 (MQ=255)                                                                                      gatgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAtttt                                                                                                                                                          <  2:62596/250‑1 (MQ=255)                                                                                         gatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCg                                                                                                                                                      >  1:438771/1‑251 (MQ=255)                                                                                               gagaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGACCAAACCGGTCAATCATATTCTTTAAA‑CC‑GCAATCATCGCTCTCCGGAGACAGACGGCCTTAACGAGACCGAGCAAAGAG‑CATATCATTGATGTTTACCTCAGAGGGGATGGCAGGACCGACAAACAAAAATGCACCTTTGTTTTTCAAACTGAAGCGTTAGCGTGGGTAAAAAATTTTGGTGACCATTTTACCGTTGaat                                                                                                                                                >  1:402713/1‑249 (MQ=255)                                                                                                             aaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑TTTAAGACCA‑CAAACGGCGCTCTCCGGGGGCAGATGACATTAACGAGACTGAGAAAAAAG‑CCTATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGGTACCATTTTTCCGTTGCCTCAATAAGGCAgctg                                                                                                                                  >  2:437236/1‑248 (MQ=255)                                                                                                                                            gggAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTcc                                                                                                   <  2:337499/251‑1 (MQ=255)                                                                                                                                            gggAACCAATCGGGTCATTCATCTT‑AATAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTcc                                                                                                   <  1:430058/251‑1 (MQ=255)                                                                                                                                                  cAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATc                                                                                             <  2:366571/251‑1 (MQ=255)                                                                                                                                                         gTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATc                                                                                      >  1:10693/1‑251 (MQ=255)                                                                                                                                                              tcatCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCatat                                                                                  <  1:162754/250‑1 (MQ=255)                                                                                                                                                                    t‑attaAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTc                                                                             <  1:176092/249‑1 (MQ=255)                                                                                                                                                                              cAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACt                                                                   <  1:144667/251‑1 (MQ=255)                                                                                                                                                                                      gACGCTCTCCGGATTCAGATGACATTAACGAGACTGAGCAAAGCG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGATTTTCAACCTGAAGAGTCAGCGCGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGGGGGACGCAGCGATAAAACGTGCCGGGCTCCGCCATCGAAACTCATATCAGTCCAGACATACTTATGccgg                                                           >  2:266359/1‑249 (MQ=255)                                                                                                                                                                                         gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATAGCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTGTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGGCCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAGCGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgctg                                                        <  2:505656/251‑1 (MQ=255)                                                                                                                                                                                                    tCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTAAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCAGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACa                                              >  2:445785/1‑250 (MQ=255)                                                                                                                                                                                                     cAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGCTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGc                                            >  2:113554/1‑251 (MQ=255)                                                                                                                                                                                                     cAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGc                                            >  2:415397/1‑251 (MQ=255)                                                                                                                                                                                                        aTGACAATAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTgg                                         <  2:15177/251‑1 (MQ=255)                                                                                                                                                                                                                        tGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATg                                                                                                                                                                     >  2:470589/1‑111 (MQ=255)                                                                                                                                                                                                                        tGAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATg                                                                                                                                                                     <  1:470589/111‑1 (MQ=255)                                                                                                                                                                                                                         gAGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAtt                          <  1:7610/249‑1 (MQ=255)                                                                                                                                                                                                                          aGCAAAGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTtat                       >  1:339719/1‑251 (MQ=255)                                                                                                                                                                                                                              aaGAG‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGc                    >  2:300175/1‑250 (MQ=255)                                                                                                                                                                                                                               cgag‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGAAAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGAACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAACCGGCATTTATAGCa                   >  2:121899/2‑250 (MQ=255)                                                                                                                                                                                                                                 ag‑CATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGACAAACAAACATGCACCTTTGCTTTTAAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAc                >  1:166302/1‑251 (MQ=255)                                                                                                                                                                                                                                        tccttGATGTTCACCTCAGAGAGTCTGGCAGAACAGACAACCAAAAATGCACCTTTGTTTTTCAAGCTGAAAGGTCAGCGAGAGAAAAAACTTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAACGCTAAAACGTGCCGGACTCCGCAATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGGCGCTGACAGCTGGTGCTAACCCGGCCTTTATATCAAACCAGATg          >  2:163229/4‑251 (MQ=255)                                                                                                                                                                                                                                          aTTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATggg        >  2:176906/1‑251 (MQ=255)                                                                                                                                                                                                                                           ttGATGTTCACCTCAGCGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTACAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGACGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGc       >  1:127741/1‑251 (MQ=255)                                                                                                                                                                                                                                             gATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCAGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGcc      >  2:517674/1‑250 (MQ=255)                                                                                                                                                                                                                                             gATGTTCACCTCAGAGAGTATGGCAGAAAAGAAAAACAAAAAGGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAACCGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGCTACATACTTATGCCTGCTGGTCGCTGACATCTGGTGTTAACCCGGCATTTATAGCAAACCAGATGGGCCa     <  1:5836/251‑1 (MQ=255)                                                                                                                                                                                                                                                gTTCACCTCAGAGCGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGc  >  1:13165/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     GGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTT‑ATTAAG‑CCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAG‑CATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGC  >  NC_000913/1652089‑1652575

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