BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,161,391	A→G	T90A (ACC→GCC)	fabR →	transcriptional repressor of fabA and fabB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,391	0	A	G	100.0%	411.2 / NA	113	T90A (ACC→GCC)	fabR	transcriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (53/60); total (53/60)

Only 100 of 113 total aligned reads displayed.

AGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTAT  >  NC_000913/4161243‑4161537

aGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGcc                                                                                                                                                  >  4:274293/1‑151 (MQ=255)

   gggCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCg                                                                                                                                               <  2:1321297/151‑1 (MQ=255)

   gggCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCg                                                                                                                                               <  3:297333/151‑1 (MQ=255)

    ggCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCgg                                                                                                                                              <  4:150002/151‑1 (MQ=255)

      cATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGtc                                                                                                                                            <  2:1065323/151‑1 (MQ=255)

       aTTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGtct                                                                                                                                           <  2:532630/151‑1 (MQ=255)

             cccACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCcaca                                                                                                                                      <  2:1089597/150‑1 (MQ=255)

             cccACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACAt                                                                                                                                     <  3:774518/151‑1 (MQ=255)

             cccACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACAt                                                                                                                                     >  3:919656/1‑151 (MQ=255)

                aCCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACAttt                                                                                                                                   >  2:380958/1‑150 (MQ=255)

                aCCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTa                                                                                                                                  <  3:516583/151‑1 (MQ=255)

                aCCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTa                                                                                                                                  <  3:704310/151‑1 (MQ=255)

                     tttttATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATgg                                                                                                                               >  2:896669/1‑149 (MQ=255)

                      ttttATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGa                                                                                                                              <  3:720595/149‑1 (MQ=255)

                       tttATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGtt                                                                                                                           >  2:605663/1‑151 (MQ=255)

                       tttATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGtt                                                                                                                           >  2:246560/1‑151 (MQ=255)

                        ttATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGTGTGTGATCCGCGCCTCGGCCTCCACATTTATGGAGTTc                                                                                                                          >  3:612138/1‑151 (MQ=255)

                         tATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTc                                                                                                                          >  1:448832/1‑150 (MQ=255)

                          aTCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCAt                                                                                                                        >  2:1010425/1‑151 (MQ=255)

                            cGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCg                                                                                                                      >  3:1190229/1‑151 (MQ=255)

                            cGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCg                                                                                                                      >  1:160348/1‑151 (MQ=255)

                            cGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCg                                                                                                                      <  4:228311/151‑1 (MQ=255)

                               cATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGta                                                                                                                   <  3:1255666/151‑1 (MQ=255)

                                aTTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCTTGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGta                                                                                                                   >  2:315586/1‑150 (MQ=255)

                                 tttCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGtaa                                                                                                                  <  4:112393/150‑1 (MQ=255)

                                  ttCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGtaa                                                                                                                  >  3:32900/1‑149 (MQ=255)

                                    ccGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTaataat                                                                                                              <  3:1135846/151‑1 (MQ=255)

                                    ccGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTaataat                                                                                                              <  3:246541/151‑1 (MQ=255)

                                        gACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCt                                                                                                           >  4:222779/1‑150 (MQ=255)

                                        gACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTa                                                                                                          <  4:257553/151‑1 (MQ=255)

                                           gTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACg                                                                                                       >  3:1163929/1‑151 (MQ=255)

                                              gACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCt                                                                                                    <  2:1267039/151‑1 (MQ=255)

                                                cGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTc                                                                                                  <  3:344926/151‑1 (MQ=255)

                                                cGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTc                                                                                                  >  3:702037/1‑151 (MQ=255)

                                                  aaCTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCg                                                                                                <  1:61269/151‑1 (MQ=255)

                                                   aCTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCAGGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTGCCgg                                                                                               <  2:802863/151‑1 (MQ=255)

                                                      gggTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGtt                                                                                             <  1:397642/150‑1 (MQ=255)

                                                       ggTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttat                                                                                           >  3:278511/1‑151 (MQ=255)

                                                       ggTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGttat                                                                                           <  1:176993/151‑1 (MQ=255)

                                                          cTGCCCATGGTTGATGAGAGCGGTTTAAGGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATACGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTtattat                                                                                        >  3:1050369/1‑151 (MQ=255)

                                                            gACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTtattat                                                                                        <  1:302021/149‑1 (MQ=255)

                                                            gACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTg                                                                                      <  2:352020/151‑1 (MQ=255)

                                                             aCCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGc                                                                                     >  1:156195/1‑151 (MQ=255)

                                                                aTGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCggg                                                                                  >  2:526740/1‑151 (MQ=255)

                                                                aTGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCggg                                                                                  >  2:91421/1‑151 (MQ=255)

                                                                  ggTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGaa                                                                                <  2:955073/151‑1 (MQ=255)

                                                                  ggTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGa                                                                                 <  3:500329/150‑1 (MQ=255)

                                                                   gTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAAc                                                                               >  2:954376/1‑151 (MQ=255)

                                                                    ttGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACg                                                                              >  4:454767/1‑151 (MQ=255)

                                                                      gatgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGc                                                                             <  3:222784/150‑1 (MQ=255)

                                                                       atgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGc                                                                             >  2:53095/1‑149 (MQ=255)

                                                                       atgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTc                                                                           <  2:343376/151‑1 (MQ=255)

                                                                        tgaGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTcc                                                                          >  2:1174511/1‑151 (MQ=255)

                                                                         gagaGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCg                                                                         <  3:86879/151‑1 (MQ=255)

                                                                          agagCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCgg                                                                        <  3:438799/151‑1 (MQ=255)

                                                                          agagCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCgg                                                                        >  2:147537/1‑151 (MQ=255)

                                                                          agagCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCgg                                                                        <  2:786705/151‑1 (MQ=255)

                                                                              cGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCAcc                                                                    >  3:1338938/1‑151 (MQ=255)

                                                                                gTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTTCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCCTTCGGTTATTATTGCGGGGACGCTTCGGCAcc                                                                    >  2:1278436/1‑149 (MQ=255)

                                                                                gTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTc                                                                  <  3:605690/151‑1 (MQ=255)

                                                                                 tttAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTcc                                                                 >  4:260752/1‑151 (MQ=255)

                                                                                  ttAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCg                                                                >  3:42860/1‑151 (MQ=255)

                                                                                    aaTGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGGGATCCGCGCCTCGGTTTCCACATTTAGGGAGTTCATCGGTAATAATCCTAACGCCGTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCt                                                              <  4:19800/151‑1 (MQ=255)

                                                                                          aCGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGttt                                                        >  2:976896/1‑151 (MQ=255)

                                                                                           cGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTTATCCGCGCCTCGGTTTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGtt                                                         <  2:914458/149‑1 (MQ=255)

                                                                                            gCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCg                                                      <  1:424094/151‑1 (MQ=255)

                                                                                                 cTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTgccg                                                 >  2:423553/1‑151 (MQ=255)

                                                                                                  tCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTgccgc                                                >  1:167393/1‑151 (MQ=255)

                                                                                                      gcgcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGtt                                            <  4:167391/151‑1 (MQ=255)

                                                                                                        gcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTgc                                          <  2:585932/151‑1 (MQ=255)

                                                                                                        gcCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTg                                           <  2:612111/150‑1 (MQ=255)

                                                                                                          cAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTgcgc                                        >  3:588002/1‑151 (MQ=255)

                                                                                                          cAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTgcgc                                        >  3:390333/1‑151 (MQ=255)

                                                                                                          cAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTgcgc                                        >  1:206274/1‑151 (MQ=255)

                                                                                                           aGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGcgcg                                       >  3:1056973/1‑151 (MQ=255)

                                                                                                            ggCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGt                                      <  3:687922/151‑1 (MQ=255)

                                                                                                             gcgcGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTg                                     <  3:985397/151‑1 (MQ=255)

                                                                                                              cgcgTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTg                                     >  3:139443/1‑150 (MQ=255)

                                                                                                                 gTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGaaa                                  <  2:32893/150‑1 (MQ=255)

                                                                                                                  tCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAAtt                                >  3:443252/1‑151 (MQ=255)

                                                                                                                   cAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAAtt                                <  4:66284/150‑1 (MQ=255)

                                                                                                                   cAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTc                               <  2:1241628/151‑1 (MQ=255)

                                                                                                                    aGCGTATCGCCAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTc                               >  2:339491/1‑150 (MQ=255)

                                                                                                                              cAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCAtt                    >  3:1253086/1‑151 (MQ=255)

                                                                                                                              cAAAGGCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCAtt                    <  3:1174567/151‑1 (MQ=255)

                                                                                                                                  ggCGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCgg                >  2:538824/1‑151 (MQ=255)

                                                                                                                                    cGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGaa              <  1:61299/151‑1 (MQ=255)

                                                                                                                                    cGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGaa              <  3:1010468/151‑1 (MQ=255)

                                                                                                                                    cGGGAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGa               >  3:1315489/1‑150 (MQ=255)

                                                                                                                                      ggAGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACt            <  2:278484/151‑1 (MQ=255)

                                                                                                                                        aGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTg          <  1:202928/151‑1 (MQ=255)

                                                                                                                                        aGTGTGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTg          <  2:1190172/151‑1 (MQ=255)

                                                                                                                                         gtgtGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGc         >  4:152993/1‑151 (MQ=255)

                                                                                                                                         gtgtGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGc         <  3:315616/151‑1 (MQ=255)

                                                                                                                                         gtgtGATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGc         >  3:1349434/1‑151 (MQ=255)

                                                                                                                                             gATCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGAc     >  3:1225481/1‑151 (MQ=255)

                                                                                                                                               tCCGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACt    >  2:776032/1‑150 (MQ=255)

                                                                                                                                                ccGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGa               >  2:1112210/1‑138 (MQ=255)

                                                                                                                                                ccGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGa               <  3:1112263/138‑1 (MQ=255)

                                                                                                                                                ccGCGCCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTAt  <  3:380994/151‑1 (MQ=255)

AGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCGTTGCGCGTGAAATTCAGCACTTCATTGCGGAACTTGCGGACTAT  >  NC_000913/4161243‑4161537

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