Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,188,087	T→C	T110A (ACG→GCG)	yehA ←	putative fimbrial adhesin YehA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,188,087	0	T	C	100.0%	180.4 / NA	51	T110A (ACG→GCG)	yehA	putative fimbrial adhesin YehA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/28);  total (23/28)

TATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAA  >  NC_000913/2187849‑2188321                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           tATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAAcc                                                                                                                                                                                                                                >  1:611356/1‑251 (MQ=255)    gACAATGTCAGATTAAAGGAATAGTCACCTGCATTTTGTGCCAATGCAGTTCGGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGAGGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACt                                                                                                                                                                                                                             <  1:556694/251‑1 (MQ=255)          gTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGa                                                                                                                                                                                                                        >  1:845371/1‑250 (MQ=255)              gATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTg                                                                                                                                                                                                                   >  1:732709/1‑251 (MQ=255)               aTTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTAGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTACATTACTGGCCACAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTgg                                                                                                                                                                                                                  >  2:866271/1‑251 (MQ=255)                  aaaGGAATAGACACCTGCATTTTGAGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTc                                                                                                                                                                                                                <  1:204423/250‑1 (MQ=255)                  aaaGGAATAGAAACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACAGGCCGCAGGTGTCCCCGCACCTGAATAACCACTGTTGAAACTGGTCt                                                                                                                                                                                                               >  2:123121/1‑251 (MQ=255)                          aaacacCTGCCTTATGTGCCAATGCATTTCTTTAGCTATTTACACATGCATCCCATGTCGTTCCACTAAATAAACAAAGTAGAAATGGTACTGCGGCACTCCCGCGTGTCGATGATGACGCACGGGATAATGACTCAACCCC‑ATTAATCATCCCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCAACGCAACTGAATAACCACTGGTGAAACTGGTCTGattatt                                                                                                                                                                                                        <  2:28912/248‑1 (MQ=255)                               ccTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCa                                                                                                                                                                                                                               >  1:591175/1‑222 (MQ=255)                                             cAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGtttt                                                                                                                                                                                    <  2:487780/251‑1 (MQ=255)                                                       cTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGGGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGCGCTGTTGTTTTTATTCGTccc                                                                                                                                                                          >  2:450686/1‑250 (MQ=255)                                                                 ttACACATGCACCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCCCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATg                                                                                                                                                                 <  1:510686/250‑1 (MQ=255)                                                                   acacATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGACACTCCGGCGTGTCGATGATGAAGCACTGGATAATGTCGCAAAACC‑ATTAATAATCACTGAACCGCCGCCCCCCGCAGAGACATTCACTTTTGCATTACTGGACGCAGGTGTCACCGCTACTGAATAACCACTGGTGGAACTGGTCTGATTATTTGACAACGTGTTGTTGATTTAGTTGGTCCATTCAACATGTa                                                                                                                                                               >  1:874532/1‑250 (MQ=255)                                                                         gCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCTCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAaatttaa                                                                                                                                                        <  1:855891/251‑1 (MQ=255)                                                                          cATCCCATGACTTGCCACCTAATAAACACAGTAGACATTGTACTGAGGCACTCACGCGTGTCGCTGATGAAGCCCTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACAGCCGCCCCCCGCAGAAAAATTCACTTTTGCATTACTGGCCGCCGGTGTCACCGCAACTGAATAACAACTGGTGAAACTGATCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAATATGTAaatttaat                                                                                                                                                       >  2:376893/1‑251 (MQ=255)                                                                                  ttcTTGCCACCTTATAAAACCAGAAAAAATTTTACTGCGGCCCACCCTCGGCTCGATGATGAAGCACTGGATAATGTC‑CAAAACCCATTAATCATCACTGAACCCCCGCCCCCCGCAGAAACATTCACTTTTGCATTTCTTGCCGCAGGTGTCACCGCACCTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTGTTATTGGTCCATTCAAAATGTAAATTTAATTTTATgg                                                                                                                                                <  2:678184/248‑1 (MQ=255)                                                                                    cTTGCCCCCTAATAAACACAGTAGAAATGTTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCAtt                                                                                                                                                                        <  2:554923/224‑1 (MQ=255)                                                                                    cTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCAtt                                                                                                                                                                        >  1:554923/1‑224 (MQ=255)                                                                                      tGCCTCCTAATAAGCACAGTAGAAATTTTTGCGCGTCACTACCCCGTGTAGTTGATGAATCACTGGATAATGTCGCAACACT‑ATTAATCCTCACTGAACAGCCTCCCACCGCAGAAACATTCTCTTTTGCATTACTTGCCGCATGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATTTTGTTTTTATTGGTCTAGTCAAAAGGTTAATTTAATGTTATGGTATc                                                                                                                                            <  2:106689/250‑1 (MQ=255)                                                                                        ccACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAACTTTAATTTTATGGTATCGtt                                                                                                                                         >  1:625626/1‑251 (MQ=255)                                                                                             taataaACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATaa                                                                                                                                     >  1:159932/1‑250 (MQ=255)                                                                                                        gTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACGGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACaata                                                                                                                         <  1:28254/251‑1 (MQ=255)                                                                                                            aaaTTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAAGTTAATTTTATGGTATCGTTATAAAGTCCAACAATAatat                                                                                                                     <  1:320683/251‑1 (MQ=255)                                                                                                            aaaTTGTACGGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCTCAACCCC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAataataa                                                                                                                     <  2:626569/251‑2 (MQ=255)                                                                                                              aTTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATcc                                                                                                                   <  1:473443/251‑1 (MQ=255)                                                                                                              aTTGTACTGCGGCACTCCCGCGTGTCGAAGATGAAGCACTGGATAATGACGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTTGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATcc                                                                                                                   <  2:412253/251‑1 (MQ=255)                                                                                                                       cgacaccccccagtgtCTATGATGAAGCACGGGATAAGGTAGCAACCCC‑ATTAATCATCCCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGTTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATGCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGcgc                                                                                                          <  2:403438/239‑3 (MQ=255)                                                                                                                                        gatgatGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTTTAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTAtt                                                                                         <  1:752707/251‑1 (MQ=255)                                                                                                                                        gatgatGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTAtt                                                                                         >  1:650756/1‑251 (MQ=255)                                                                                                                                         atgatgAAGCACTGGATAATGTCGCAACACA‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTAttt                                                                                        >  1:258509/1‑251 (MQ=255)                                                                                                                                             tgaAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACggg                                                                                                             >  1:819884/1‑226 (MQ=255)                                                                                                                                                      ggATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCAt                                                                            <  1:565172/250‑1 (MQ=255)                                                                                                                                                       gATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCTCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAATTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATGTGTACTTACGCATGa                                                                          <  2:737843/251‑1 (MQ=255)                                                                                                                                                                ccAACACC‑ATTAATCATCCCGGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCAtt                                                                 <  2:611356/250‑1 (MQ=255)                                                                                                                                                                 cAACACC‑ATTAATGATCACTGAACCGCCGCCACCCGCAGAAACATTCACTTTTGCATTACTCGCCGCAGATGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATCTGACAACGTGATGTTGTTTTTATTGGTCAATTAAAAATGTAAATTTAATTTTATGGTATCGTTATAATGTCCAACAATAATATCCGACGGGCTAATAGGGGTGATAGTATTTGTACTTACGTATGACGatacatt                                                                 >  2:800534/1‑244 (MQ=255)                                                                                                                                                                      cc‑ATTAATCATCCCTTAACCGCCGCCCCCCCCAGAAACATTCACTTTTGCATTCCTGGCTGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAAtt                                                           <  2:781427/251‑1 (MQ=255)                                                                                                                                                                         aTTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTAt                                                         <  1:105145/251‑1 (MQ=255)                                                                                                                                                                         aTTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTAt                                                         >  1:101352/1‑251 (MQ=255)                                                                                                                                                                               catcaCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCa                                                    <  1:661277/250‑1 (MQ=255)                                                                                                                                                                                   aCTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGt                                               >  1:676103/1‑251 (MQ=255)                                                                                                                                                                                           gccgccCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACTGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCAttt                                       <  1:510579/251‑1 (MQ=255)                                                                                                                                                                                                      cAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGCATCATTAAAATTATGCTCATATg                                                >  2:314034/1‑231 (MQ=255)                                                                                                                                                                                                      cAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATg                                                <  1:314034/231‑1 (MQ=255)                                                                                                                                                                                                         aaaCATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGt                         >  1:621713/1‑251 (MQ=255)                                                                                                                                                                                                                    tttGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATcc                <  1:603323/249‑1 (MQ=255)                                                                                                                                                                                                                     ttGCATTACTGGCCGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATcccc              <  1:450686/250‑1 (MQ=255)                                                                                                                                                                                                                              tGGCCGCAGGGGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAAcc    >  1:791880/1‑251 (MQ=255)                                                                                                                                                                                                                                gCCGCAGGTGTCACCGCAACTGAATAACCACTGGTTAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCaa  <  1:690298/251‑1 (MQ=255)                                                                                                                                                                                                                                 ccGCAGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCa   >  1:44151/1‑249 (MQ=255)                                                                                                                                                                                                                                     aGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCAtt                                                                 >  2:398967/1‑183 (MQ=255)                                                                                                                                                                                                                                     aGGTGTCACCGCAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCAtt                                                                 <  1:398967/183‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                           TATGACAATGTCAGATTAAAGGAATAGACACCTGCATTTTGTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACC‑ATTAATCATCACTGAACCGCCGCCCCCCGCAGAAACATTCACTTTTGCATTACTGGCCGCAGGTGTCACCGTAACTGAATAACCACTGGTGAAACTGGTCTGATTATTTGACAACGTGATGTTGTTTTTATTGGTCCATTCAAAATGTAAATTTAATTTTATGGTATCGTTATAAAGTCCAACAATAATATCCGACGGGCTCATAGGGGTGATAGTATTTGTACTTACGCATGACGAATCATTAAAATTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAA  >  NC_000913/2187849‑2188321

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