Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,894,839	T→C	L12P (CTC→CCC)	pabB →	aminodeoxychorismate synthase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,894,839	0	T	C	100.0%	147.1 / NA	41	L12P (CTC→CCC)	pabB	aminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/14);  total (27/14)

CGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACTCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGC  >  NC_000913/1894607‑1895073                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              cGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCCGGAcgccgc                                                                                                                                                                                                                            >  2:290503/1‑251 (MQ=255)  gCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCATGG‑CGTCAGGACGCCGCt                                                                                                                                                                                                                           >  1:376644/1‑251 (MQ=255)      ctctTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAAt                                                                                                                                                                                                                       >  1:227511/1‑251 (MQ=255)                              cTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCAcc                                                                                                                                                                                               >  1:714042/1‑251 (MQ=255)                                ggCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCGTAAGCCACCt                                                                                                                                                                                              <  2:501156/250‑1 (MQ=255)                                           aTCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGAt                                                                                                                                                                                  >  1:589012/1‑251 (MQ=255)                                                       aTCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAACCCTGCAAACATTGTTAACTCCTGCTCAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGATGCCGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGACCGCCGCTGAATTTTATTTCTCCCGCTTACGCCACCTGCCGTGGGCGATGCTTTTACACTc                                                                                                                                                                      >  2:362464/1‑251 (MQ=255)                                                       aTCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTc                                                                                                                                                                      <  2:227511/251‑1 (MQ=255)                                                                                             ttCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGTTATtg                                                                                                                                >  2:389227/1‑251 (MQ=255)                                                                                               cATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATtgtg                                                                                                                              >  1:505434/1‑251 (MQ=255)                                                                                                aTGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTgg                                                                                                                             >  1:406260/1‑251 (MQ=255)                                                                                                 ggggggagggaAGGTAAACCTGCAAACATTGTTATCTCCTTCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTCCTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCCCCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGTTATTGTggg                                                                                                                            <  2:170268/242‑2 (MQ=255)                                                                                                      gagtgaAGGTACACCTGCAAACATTGTTAACTCCTGCTCAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAATCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCgccg                                                                                                                      >  2:90972/1‑251 (MQ=255)                                                                                                      gagtgaAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCgccg                                                                                                                      <  1:301450/251‑1 (MQ=255)                                                                                                          gaAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCgccgagcc                                                                                                                  >  1:513119/1‑251 (MQ=255)                                                                                                             ggTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGCCTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTCCTTTACCCTGG‑CGTCCGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGAt                                                                                                               >  1:513637/1‑251 (MQ=255)                                                                                                                aaaCCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGAttt                                                                                                             <  1:195768/250‑1 (MQ=255)                                                                                                                  aCCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCAGAGCCGATTT‑Gcc                                                                                                         >  1:758592/1‑250 (MQ=255)                                                                                                                   ccTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTCCCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCAc                                                                                                        <  1:102081/251‑1 (MQ=255)                                                                                                                   ccTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCAc                                                                                                        <  1:628680/251‑1 (MQ=255)                                                                                                                       cAAACATTGTTAAATCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGG‑GATTACTTTACCCTGG‑CGTCAGGTCGCCGCTGAATTTTATTTCTCCCGCTTACGCCAACTGCCGTGGGCGATGCTTTTACACTCCGGCTATGACGATCATCACAATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCTGTTT‑GCACttt                                                                                                     >  2:781113/1‑250 (MQ=255)                                                                                                                        caaCATTGATAACTCCTGCTAAATTGTTGGCGCTAATTCTTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCATTTACACTCCGGCTATGCCGCTCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTa                                                                                                    >  2:72038/2‑250 (MQ=255)                                                                                                                        aaaCATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGACGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGACTATGCCGATCATCCGAATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTa                                                                                                    >  2:469965/1‑250 (MQ=255)                                                                                                                            aTTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCAc                                                                                               <  1:690049/251‑1 (MQ=255)                                                                                                                             ttggacaccccTGATAAATGGTTGGCGCTAATTATTGCATTCTACCCGTCCCAGAGCGATGAGAGTCA‑G‑ACTTTATGAAGACGTTATCTCCTTCGTT‑GCTTACTTT‑CCCTGGGCGTCAGGACGCCGATGATTTTTATTCCGCCCGCTAAAGCCACCTGCCGTGTGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATGGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCAc                                                                                               <  2:578098/241‑1 (MQ=255)                                                                                                                               tgtAACCCCTGCTAAATGGTTGGCGATACTTATTCCATTCTACCGGTCACATAGCCATTAGAGTAA‑GGCATG‑ATGAAGACTTTATTCCCCGCTTT‑GTTTCCTTTACCCTGG‑CGTGAGGACGCCGCTGATTTTTATTTCTCCCGCTTAAGCCACCTGTCGTGGGCGTTGCTTTTACACTCCGGCTATGCCGATCACCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCATTTTAACCACtt                                                                                             <  2:758592/248‑1 (MQ=255)                                                                                                                                    cTCCTGCTAACTTGTTGACGCTAATTATTTCATGCTACCCGTCCCATAGCCAGTAGAGGCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCAGCTTAAGACACCTGCCGTAGGCGATGCTTTTAAACTACGGCTATGCCGATCATCAGTATAGCCGCTTTGATATTGTGG‑TCGCCGCGCCTATTT‑GCACTTTAACCACTTTCGGt                                                                                        >  2:514251/1‑250 (MQ=255)                                                                                                                                      ccTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCAATTTCGGTaa                                                                                      >  2:239304/1‑250 (MQ=255)                                                                                                                                            aaaTTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCACTTTCGGTAAAGAAAcc                                                                               >  2:9134/1‑251 (MQ=255)                                                                                                                                                      cgcTAATTATTTCATGCTACCCGGCACCTAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCAGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGCATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGATCCGATTT‑GCCCTTTACCCACTTTCGGTCAAGAAACCGTTGTTAGTg                                                                     >  2:37950/1‑251 (MQ=255)                                                                                                                                                                  catactcccaagCACATAGACAGTAGAGTCA‑GGACTG‑ATCACATCCTTATCTCCCGCTGT‑GATTACTTTCCCCTGG‑CGTCAGGACGCCGATGAATTTTATTTCTCCCGCTTACGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATAATCCAAATAACCGCTTTGATATTGTGG‑TCGC‑AAACCAATTTTGCACTTTAAACACTTTCGGTAAAGAAACCTTTGTTAGTAAAAACAAAAAAc                                                         >  2:91023/12‑251 (MQ=255)                                                                                                                                                                         cccGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAAc                                                  >  1:437069/1‑251 (MQ=255)                                                                                                                                                                          ccGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑GGTCAGGACGCAGCTGAAGTTTATTTCTCCCGCTTAAGCCACCTGCGGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCACCCGGATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCTCTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACg                                                 <  1:362464/251‑1 (MQ=255)                                                                                                                                                                           cGGCACATGCCCATTAGGGTCA‑GGACGA‑ATGAAAACGTTATCCCCGGCTTT‑GATTACTTTCCCCTGG‑CGTCGGGACGCCGCTGATTTTTTTTTCTCCCCCTTAAGCCCCCTGCCGGGGGCGATGTTTTTACACTCCGGCTATGCCAATCATCCGTATAGCCGCTTTGATAT‑GTGGGTCGCCGAGCCGATTT‑GCACTTTAACCACTTTCGGTAA‑GAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGAc                                               <  1:99942/251‑1 (MQ=255)                                                                                                                                                                              cacaTAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCACTTTCGGTAAAGAAACCgttgtt                                                                         <  1:285448/223‑1 (MQ=255)                                                                                                                                                                              cacaTAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCACTTTCGGTAAAGAAACCgttgtt                                                                         >  2:285448/1‑223 (MQ=255)                                                                                                                                                                               acaTTGCCAGTAGAG‑CACGGACTT‑ATGACGACGTTCT‑TCCCCGTGTTGATTTCTTTACCCTGG‑CGTCAGGACGCCGTTGAATTTTATTTCTCCCGCTTAAGCCCCCTGCCGTGGGCGATGCTTTTACACGCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCAc                                            <  2:585122/251‑1 (MQ=255)                                                                                                                                                                                       aGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACACCGACCACTGATAAcc                                    >  2:402242/1‑251 (MQ=255)                                                                                                                                                                                                a‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAgg                           >  1:678525/1‑251 (MQ=255)                                                                                                                                                                                                                  ttATCTCCCGCTGT‑GATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCt           >  1:299662/1‑251 (MQ=255)                                                                                                                                                                                                                              gt‑gATTACTTTACCCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTCGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCgcgc  >  2:669185/1‑248 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              CGCAGCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCA‑GGACTG‑ATGAAGACGTTATCTCCCGCTGT‑GATTACTTTACTCTGG‑CGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGG‑TCGCCGAGCCGATTT‑GCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGC  >  NC_000913/1894607‑1895073

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