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%	167.1 / NA	50	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/30);  total (20/30)

TATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652086‑1652573                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    tATGGGGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCag                                                                                                                                                                                                                                                 <  2:460917/250‑1 (MQ=255)    ggTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCagagag                                                                                                                                                                                                                                             <  2:192484/251‑1 (MQ=255)    ggTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATAGACGATCTCCGGAGTCAGATGACATTAAAGA‑GACTGAGCAAAAAGCATATCATTGATGTTCACCTCaggg                                                                                                                                                                                                                                               >  1:528173/1‑247 (MQ=255)               ggttaaaCTGGAGTTCGCCCTGGTGAGATGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGAAGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGc‑agaac                                                                                                                                                                                                                                 <  2:672151/245‑1 (MQ=255)               ggttATACTGGAGTTCGCCCTGTTGAGCTGTGGGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGACGATCAGGAGAAATTTAGCAAACGACCGCTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAAGCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GGCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGc‑agaa                                                                                                                                                                                                                                  <  2:507117/248‑1 (MQ=255)                tttATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAACCAGCAATCGAAGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGGGCCAAGAGCATATCATTGAAGTTCACCTCCGAGAGGATGGc‑agaac                                                                                                                                                                                                                                 >  1:752382/1‑250 (MQ=255)                     aCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCGCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGaaaa                                                                                                                                                                                                                           <  1:577471/251‑1 (MQ=255)                           gTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGAGCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACaaaa                                                                                                                                                                                                                     <  1:491198/251‑1 (MQ=255)                              cGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTGCACGGTACCAAAAACACAGGCGGGAACCAAGCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAAt                                                                                                                                                                                                                   <  2:133491/250‑1 (MQ=255)                              cGCCCTGGTGAGCTGTGCGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCCGATGACATTAACGC‑GACTGAGCAAAAAGCATATCATTGATGTTCACCTCCGAGAGTATACC‑AGAACAAAAAAAACAAAATg                                                                                                                                                                                                                  >  1:764885/1‑251 (MQ=255)                               gCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGc                                                                                                                                                                                                                 <  1:724024/251‑1 (MQ=255)                                cccTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAAGGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAGCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACGGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCa                                                                                                                                                                                                                <  2:728271/251‑1 (MQ=255)                                             tgtgCACTTGGATGGGAGGACATAGATATGAAAAATGGAACAATGATGAGCAGGAGAAATTTAGCAAAAGACCGGTTCACGGTACCAAAAACACAGGTGGGAACCAATCGGGTCATTCATCTTATTAATCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCGTTGTTTTTc                                                                                                                                                                                                    <  2:986078/250‑1 (MQ=255)                                               tgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCCTTTATGTTCACCCCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAc                                                                                                                                                                                                 >  1:463525/1‑251 (MQ=255)                                                gCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGA‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAc                                                                                                                                                                                                 <  1:556096/250‑1 (MQ=255)                                                          gggAGGAAATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCGTTGTTTTTCAACCTGAAGTGTc                                                                                                                                                                                       <  2:423402/250‑1 (MQ=255)                                                             aggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCg                                                                                                                                                                                   <  1:618406/251‑1 (MQ=255)                                                                                  gAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAc                                                                                                                                                              <  1:821114/251‑1 (MQ=255)                                                                                        tgatgaTCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCCATAGGGGCATTCATCTTATTAAGCCAGCCAACGACACTCTCCGGAGTCAGATGGCATTATCGA‑AGCTGCGCCAAGAACATTTCAATGGTTTTTACCCCAAAGAGTATGGCCAGA‑CAAAAAAACAAATATGCAACTTTGTTTTTTAAACTAGAGTGTCATCGAGAATAAAAAAatattgttaaccattt                                                                                                                                                        >  2:328744/1‑239 (MQ=255)                                                                                             atCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAAAAAAAATGCCCCTTTTTTTTTTAACCTGAAGGGTCAgaggggggaaaaaaatattgggaa                                                                                                                                                              >  1:903999/1‑217 (MQ=255)                                                                                             aacAGGAGAAATTTAGCAAAAGACCGTTTCACCGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAc                                                                                                                                                              <  2:903999/238‑1 (MQ=255)                                                                                               cAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGT‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGtt                                                                                                                                                  <  2:457016/250‑1 (MQ=255)                                                                                                  gagaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCCATCGGGGCATTTATCTTATTAAACCAACAATCGACGCTCTCCCGAATCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGCGTATTGC‑AGAAAAGAAAAACAAAAATGCCCATTTTTTTTTCCAACTAAAATGTCAGCGCGAGAAAAAAAATAAGGTGACCATTTTAcccgtgact                                                                                                                                              >  1:516798/1‑244 (MQ=255)                                                                                                      aaTTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGTAACAATCGGGGCGCTCAACTTTATCAGACACCAATATAACAGCTCCGGAGTTAGACGGAACTAAAGAAGA‑TGCGCGAAGAGAATATAATTTATGTTCACCTCCCAGAGGATAGA‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCAGAATTGTCAGAGAGAGTTAAAAATTATgagtaacatattacctgtgtatcaa                                                                                                                                           >  2:516049/1‑226 (MQ=255)                                                                                                              aaaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCCTCTTATTAAGCCAGCAATCGACCCTCTCCGGGGTCAGATGACATTACCGG‑GACTGCGCAAAAAGCATATCATTTATGTTCACCTCACAGAGGATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAACTGAAGTGGCAGCGAGAGAAAAAAATTAAGGTGACCATTTTACCGTTTACTCAATAAGGc                                                                                                                                     >  1:979260/1‑248 (MQ=255)                                                                                                                     cGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGa                                                                                                                           <  2:528173/251‑1 (MQ=255)                                                                                                                     cGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAAAAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGGGACCATTTTACCGTTGACTCAATAAGGCAGATGAggg                                                                                                                            >  1:378856/1‑250 (MQ=255)                                                                                                                                aCCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATa                                                                                                                 <  1:782575/250‑1 (MQ=255)                                                                                                                                  ccaaaaCACATGCGGGAAACAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATaaaa                                                                                                              <  2:667930/249‑1 (MQ=255)                                                                                                                                      aaCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTg                                                                                                          <  1:481203/251‑1 (MQ=255)                                                                                                                                                ggAACCAATCGGGTCATTCATCTTATTAACCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGCCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCg                                                                                                <  1:931486/251‑1 (MQ=255)                                                                                                                                                                   attttattaAGCCAGCAATCGACGCTCCCCGTAGACAGATGACATTAACGA‑GACCTAGGAAAGAGCATATCAGTGATGTTCACCTCCGAGCGTATTGA‑AGAACAGAAAAACATAAATGCATCTTTGTTTTTCAACCTTAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATATGGCAGATGTTGGACGCAGCGATGAAACGTGCAGGACGCTGTCATCGAAAATCACATCAg                                                                             <  2:839338/248‑1 (MQ=255)                                                                                                                                                                       tattaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCg                                                                          >  1:489985/1‑250 (MQ=255)                                                                                                                                                                          taggCCAGCAATCGACGCTCTCAGGAGTCAGATGACATTAACGA‑GACTAAGCCAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACATTTGTTTTTCAAACTGAAGTGTCAGAGAGAGTCAAAAATTATGGGGACCATTTTACCGTTGCCTCAATAAGGCAGATGTGGGACGCAGCGATAAAGCCTGCCGGCCTCCGCCATCGCAAATCATATCAGTCGAGAc                                                                      >  2:107299/4‑251 (MQ=255)                                                                                                                                                                             gCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGCGACAt                                                                    >  2:596146/1‑250 (MQ=255)                                                                                                                                                                                    aTCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGc                                                            <  1:631974/251‑1 (MQ=255)                                                                                                                                                                                      tgACGCTATCCGGAGTCAGGTGACATTAACGA‑GACTGAGCAAAGAACATATCATTGCTGTTCACCTCAGAGAGTATGGC‑AGAACAGCAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGAAAAAAGTTCTGGAGACCATTTTACCGTTGACTCAATAAGGCGGATGTGGGACGCAGCGTAAAAAAGTGCCGGACTCAGCCATCGAAAATCATATCAGTCGAAACATACTTATgact                                                          >  2:113457/2‑248 (MQ=255)                                                                                                                                                                                       gACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctg                                                         >  1:721899/1‑251 (MQ=255)                                                                                                                                                                                       gACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctg                                                         >  1:553384/1‑251 (MQ=255)                                                                                                                                                                                               ccGGAGTCAGATGTCATTAACGA‑GACTATGCAAAGATCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCGTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACAATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGc                                                 <  2:684278/251‑1 (MQ=255)                                                                                                                                                                                                 gtaggccgATGACATTAACGA‑GACTGAGCACAGAGTATATCATTGATGTTACCCTCCGATAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGGGTCAGCTAGAGTAAAAAATTATGGTGATCACTTTAACGGTGGCTCAATAAGGCAGATGTGGTACTCAGCGATAAACCGTGCCGGACTCTGCCATCGAAAATCATATCAGTCGAGACCCACTGATGTTTGCTGGTCGCTg                                               <  2:433264/244‑1 (MQ=255)                                                                                                                                                                                                       aGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCt                                         <  1:570829/251‑1 (MQ=255)                                                                                                                                                                                                                    ga‑gaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCgg                            >  1:6591/1‑251 (MQ=255)                                                                                                                                                                                                                          tGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCCTCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCAGGCAttt                       >  2:275464/1‑251 (MQ=255)                                                                                                                                                                                                                           gAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGACAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGAGCTAACACGGCAttt                       >  2:821219/1‑250 (MQ=255)                                                                                                                                                                                                                            aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTtat                     >  1:698071/1‑251 (MQ=255)                                                                                                                                                                                                                                  gagCATATCAATGATGTTCACCTCAGATAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATGTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCaaa               <  2:353153/251‑1 (MQ=255)                                                                                                                                                                                                                                    gcagatCACTGATGTTTACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCGTTGTTTTTCAACCGGAAGTGGCAGAGAGATTAAAAAAATATGGTGATCATTTTACCGTTGACTCAATAATGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTGATGCCTGCTTGTCGCTGACAGCTGTTGCTAACCCGGCATTTATAGCaaaac             <  2:855381/247‑3 (MQ=255)                                                                                                                                                                                                                                    gCATATCATTGATGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAcc             <  1:596146/251‑1 (MQ=255)                                                                                                                                                                                                                                               aTGTTCACCTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAt  <  2:463525/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    TATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGA‑GACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGC‑AGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652086‑1652573

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