BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,031,335	T→C	W59R (TGG→CGG)	pepQ →	Xaa‑Pro dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,031,335	0	T	C	100.0%	1223.1 / NA	295	W59R (TGG→CGG)	pepQ	Xaa‑Pro dipeptidase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (156/139); total (156/139)

Only 100 of 298 total aligned reads displayed.

TCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAA  >  NC_000913/4031237‑4031432

 cAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCggg                                                                                                <  2:6640610/101‑1 (MQ=255)

 cAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCggg                                                                                                >  2:6640984/1‑101 (MQ=255)

   aGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTg                                                                                              >  1:6640987/1‑101 (MQ=255)

   aGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTg                                                                                              <  2:6640567/101‑1 (MQ=255)

   aGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTg                                                                                              <  2:6640663/101‑1 (MQ=255)

   aGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTg                                                                                              <  2:6640903/101‑1 (MQ=255)

   aGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTg                                                                                              >  2:6640989/1‑101 (MQ=255)

   aGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTg                                                                                              <  1:6640262/101‑1 (MQ=255)

      tGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCg                                                                                           >  2:6640994/1‑101 (MQ=255)

      tGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCg                                                                                           <  2:6640564/101‑1 (MQ=255)

        gATGCGTTACTTATTCCCTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGt                                                                                         <  2:6640516/101‑1 (MQ=255)

          tGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTaa                                                                                       >  1:6640999/1‑101 (MQ=255)

            cGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACt                                                                                     >  2:6641002/1‑101 (MQ=255)

             gTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTc                                                                                    >  1:6641004/1‑101 (MQ=255)

               tACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAg                                                                                  >  1:6641010/1‑101 (MQ=255)

               tACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAg                                                                                  <  1:6640619/101‑1 (MQ=255)

                aCTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAgg                                                                                 <  2:6640422/101‑1 (MQ=255)

                 cTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGt                                                                                >  2:6641014/1‑101 (MQ=255)

                 cTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGt                                                                                <  2:6640717/101‑1 (MQ=255)

                  ttATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTg                                                                               <  1:6640977/101‑1 (MQ=255)

                   tATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGc                                                                              >  1:6641017/1‑101 (MQ=255)

                     ttCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCa                                                                            <  2:6640457/101‑1 (MQ=255)

                      tCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCaa                                                                           >  2:6641022/1‑101 (MQ=255)

                      tCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCaa                                                                           <  1:6640531/101‑1 (MQ=255)

                      tCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCaa                                                                           <  1:6640650/101‑1 (MQ=255)

                      tCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCaa                                                                           >  1:6641021/1‑101 (MQ=255)

                        aCTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAAc                                                                         <  1:6640602/101‑1 (MQ=255)

                         cTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAAct                                                                        >  2:6641026/1‑101 (MQ=255)

                           ccGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAActgc                                                                      >  1:6641031/1‑101 (MQ=255)

                            cGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAActgct                                                                     <  1:6640742/101‑1 (MQ=255)

                            cGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAActgct                                                                     >  2:6641032/1‑101 (MQ=255)

                                gAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACtgctggtt                                                                 >  1:6641037/1‑101 (MQ=255)

                                gAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACtgctggtt                                                                 <  2:6640376/101‑1 (MQ=255)

                                 aGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACtgctggttg                                                                >  1:6641040/1‑101 (MQ=255)

                                 aGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACtgctggttg                                                                <  2:6640633/101‑1 (MQ=255)

                                 aGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACtgctggttg                                                                >  2:6641042/1‑101 (MQ=255)

                                   cTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACtgctggttgct                                                              >  1:6641045/1‑101 (MQ=255)

                                    tGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACtgctggttgctg                                                             <  1:6640378/101‑1 (MQ=255)

                                    tGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACtgctggttgctg                                                             >  2:6641048/1‑101 (MQ=255)

                                      ttCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACtgctggttgctggt                                                           <  2:6640436/101‑1 (MQ=255)

                                       tCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGc                                                                              >  1:6641053/1‑81 (MQ=255)

                                        cAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCtggtgg                                                         <  2:6640547/101‑1 (MQ=255)

                                         aaTGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGa                                                        >  2:6641054/1‑101 (MQ=255)

                                            gTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATgg                                                     <  2:6640356/101‑1 (MQ=255)

                                             tttttCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGc                                                    >  1:6641060/1‑101 (MQ=255)

                                              ttttCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCg                                                   >  2:6641061/1‑101 (MQ=255)

                                               tttCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGt                                                  >  2:6641064/1‑101 (MQ=255)

                                               tttCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGt                                                  <  2:6640761/101‑1 (MQ=255)

                                                ttCTCGACGATCATCCCTAACCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTg                                                 <  1:6640697/101‑1 (MQ=255)

                                                  ctcGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGaa                                               <  1:6640487/101‑1 (MQ=255)

                                                    cgacgaTCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACa                                             >  1:6641069/1‑101 (MQ=255)

                                                    cgacgaTCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACa                                             <  2:6640577/101‑1 (MQ=255)

                                                     gacgaTCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACaa                                            >  2:6641073/1‑101 (MQ=255)

                                                      acgaTCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAg                                           >  1:6641075/1‑101 (MQ=255)

                                                       cgaTCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGc                                          <  2:6640700/101‑1 (MQ=255)

                                                       cgaTCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGc                                          >  2:6641078/1‑101 (MQ=255)

                                                        gaTCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGcc                                         >  1:6641081/1‑101 (MQ=255)

                                                         atcatcCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCg                                        >  2:6641085/1‑101 (MQ=255)

                                                          tcatcCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGa                                       <  2:6640477/101‑1 (MQ=255)

                                                            atcCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGaaa                                     <  1:6640556/101‑1 (MQ=255)

                                                            atcCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGaaa                                     >  1:6641088/1‑101 (MQ=255)

                                                             tcCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAAc                                    <  1:6640704/101‑1 (MQ=255)

                                                             tcCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAAc                                    >  2:6641090/1‑101 (MQ=255)

                                                              cccTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAAc                                    <  2:6640627/100‑1 (MQ=255)

                                                              cccTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACt                                   <  2:6640496/101‑1 (MQ=255)

                                                              cccTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACt                                   >  1:6641093/1‑101 (MQ=255)

                                                               ccTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACtg                                  <  2:6640527/101‑1 (MQ=255)

                                                                cTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACtgt                                 >  1:6641098/1‑101 (MQ=255)

                                                                cTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACtgt                                 >  2:6641097/1‑101 (MQ=255)

                                                                 tATCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACtgtg                                >  1:6641101/1‑101 (MQ=255)

                                                                  aTCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTgg                               >  2:6641104/1‑101 (MQ=255)

                                                                   tCCGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGt                              >  2:6641108/1‑101 (MQ=255)

                                                                     cGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTc                            >  1:6641111/1‑101 (MQ=255)

                                                                     cGTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTc                            <  1:6640515/101‑1 (MQ=255)

                                                                      gTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCt                           >  2:6641115/1‑101 (MQ=255)

                                                                      gTTTAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCt                           >  1:6641116/1‑101 (MQ=255)

                                                                       tttAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTa                          <  1:6640859/101‑1 (MQ=255)

                                                                        ttAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTAt                         <  1:6640764/101‑1 (MQ=255)

                                                                         tAAAGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATc                        >  1:6641120/1‑101 (MQ=255)

                                                                          aaaGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCt                       <  1:6640592/101‑1 (MQ=255)

                                                                          aaaGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCt                       <  2:6640582/101‑1 (MQ=255)

                                                                            aGTGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGc                     >  2:6641123/1‑101 (MQ=255)

                                                                              tGAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCg                   <  1:6640808/101‑1 (MQ=255)

                                                                               gAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCgg                  >  1:4632701/1‑101 (MQ=255)

                                                                               gAACCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCgg                  >  2:6641127/1‑101 (MQ=255)

                                                                                aaCCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGt                 >  1:6641128/1‑101 (MQ=255)

                                                                                 aCCCGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGtt                >  2:6641132/1‑101 (MQ=255)

                                                                                   ccGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGa              <  2:6640829/101‑1 (MQ=255)

                                                                                   ccGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGa              >  1:6641136/1‑101 (MQ=255)

                                                                                    cGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGAt             >  2:6641137/1‑101 (MQ=255)

                                                                                    cGCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGAt             <  1:6640550/101‑1 (MQ=255)

                                                                                     gCAATTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTAAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTg               <  2:6640643/98‑1 (MQ=255)

                                                                                       aaTTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTAc          >  1:6641142/1‑101 (MQ=255)

                                                                                        aTTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGAtt            >  2:6641144/1‑98 (MQ=255)

                                                                                        aTTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGAt             <  1:6640434/97‑1 (MQ=255)

                                                                                        aTTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACt         <  1:6640689/101‑1 (MQ=255)

                                                                                        aTTCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACt         <  1:6640598/101‑1 (MQ=255)

                                                                                          tCAAAGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTgg       <  2:6640675/101‑1 (MQ=255)

                                                                                            aaaGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGa                                       <  1:6640391/67‑1 (MQ=255)

                                                                                              aGCGCGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGcaca   >  2:6641148/1‑101 (MQ=255)

TCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCCCTATCCGTTTAAAGTGAACCCGCAATTCAAAGCGTGGGTGCCGGTAACTCAGGTGCCAAACTGCTGGTTGCTGGTGGATGGCGTGAACAAGCCGAAACTGTGGTTCTATCTGCCGGTTGATTACTGGCACAA  >  NC_000913/4031237‑4031432

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`