BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,458,752	T→C	R86R (CGT→CGC)	b1396 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,458,752	0	T	C	100.0%	818.9 / NA	254	R86R (CGT→CGC)	b1396	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (116/138); total (116/138)

Only 100 of 266 total aligned reads displayed.

ATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGTTAA  >  NC_000913/1458662‑1458841

aTTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGccc                                                                                         >  2:5776971/1‑93 (MQ=255)

aTTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGccc                                                                                         >  2:4159645/1‑93 (MQ=255)

aTTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGccc                                                                                         <  2:5526959/93‑1 (MQ=255)

aTTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTcc                                                                                         >  2:5887407/1‑93 (MQ=255)

 ttttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCa                                                                                        <  2:4710601/93‑1 (MQ=255)

 ttttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCa                                                                                        <  1:6206834/93‑1 (MQ=255)

 ttttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCa                                                                                        >  1:3352692/1‑93 (MQ=255)

 ttttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCa                                                                                        <  1:555624/93‑1 (MQ=255)

 ttttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCa                                                                                        <  1:1917843/93‑1 (MQ=255)

 ttttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCa                                                                                        <  2:1668477/93‑1 (MQ=255)

  tttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAg                                                                                       <  2:643519/93‑1 (MQ=255)

    tcgcTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGc                                                                                     <  2:4914767/90‑1 (MQ=255)

     cACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCt                                                                                    <  2:1230687/93‑1 (MQ=255)

     cACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCt                                                                                    <  2:1162714/93‑1 (MQ=255)

       cTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCttt                                                                                  >  1:1989019/1‑93 (MQ=255)

       cTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCttt                                                                                  >  1:124761/1‑93 (MQ=255)

         ggCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGc                                                                                >  1:16289/1‑93 (MQ=255)

          gCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGcc                                                                               >  2:1104726/1‑93 (MQ=255)

          gCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGcc                                                                               <  2:6110917/93‑1 (MQ=255)

          gCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGcc                                                                               >  2:531147/1‑93 (MQ=255)

            tGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCgg                                                                             >  2:1536723/1‑93 (MQ=255)

             gATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGga                                                                            >  2:4363951/1‑93 (MQ=255)

             gATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGga                                                                            <  1:4683497/93‑1 (MQ=255)

                aCTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGac                                                                         <  1:5217497/93‑1 (MQ=255)

                    ccTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCt                                                                     <  1:112022/93‑1 (MQ=255)

                    ccTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCt                                                                     <  2:6503306/93‑1 (MQ=255)

                    ccTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCt                                                                     <  2:3655333/93‑1 (MQ=255)

                     cTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCtt                                                                    >  2:3733047/1‑93 (MQ=255)

                     cTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCtt                                                                    <  2:2777327/93‑1 (MQ=255)

                      tttGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTTGCTTCTGCCTGCACGATTGATTTTTTGTGCCCAGGCTTTGCCGGAGACACCTTa                                                                   >  1:2670541/1‑93 (MQ=255)

                      tttGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTa                                                                   >  1:2299480/1‑93 (MQ=255)

                        tGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAAc                                                                 <  1:4331859/93‑1 (MQ=255)

                        tGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAAc                                                                 <  2:5908535/93‑1 (MQ=255)

                        tGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAAc                                                                 <  1:1737949/93‑1 (MQ=255)

                          ccTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTg                                                               >  2:2866899/1‑93 (MQ=255)

                           cTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGc                                                              >  1:3707517/1‑93 (MQ=255)

                           cTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGc                                                              >  1:1042863/1‑93 (MQ=255)

                            tACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCt                                                             >  1:2056133/1‑93 (MQ=255)

                            tACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCt                                                             >  2:483693/1‑93 (MQ=255)

                            tACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCt                                                             <  2:2648080/93‑1 (MQ=255)

                             aCGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTa                                                            <  2:3479255/93‑1 (MQ=255)

                                 cTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgc                                                        >  2:3020357/1‑93 (MQ=255)

                                 cTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgc                                                        >  2:5430324/1‑93 (MQ=255)

                                   gCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgcgc                                                      <  2:84978/93‑1 (MQ=255)

                                         gCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGTGACACCTTAACTGCTACTGCGCAGGTAc                                                >  2:6734444/1‑93 (MQ=255)

                                         gCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTAc                                                >  2:5242520/1‑93 (MQ=255)

                                            aGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtc                                             >  1:4843189/1‑93 (MQ=255)

                                            aGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtc                                             <  1:3955998/93‑1 (MQ=255)

                                             ggggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtca                                            <  2:2340238/93‑1 (MQ=255)

                                             ggggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtca                                            >  2:6368353/1‑93 (MQ=255)

                                             ggggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtca                                            <  2:5962277/93‑1 (MQ=255)

                                              gggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtcat                                           >  1:4208669/1‑93 (MQ=255)

                                               ggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtcatc                                          >  1:5975985/1‑93 (MQ=255)

                                                 cTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAg                                        <  1:5427275/93‑1 (MQ=255)

                                                   ggCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAggg                                      >  2:4192658/1‑93 (MQ=255)

                                                       gCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGggaag                                  >  2:795048/1‑93 (MQ=255)

                                                       gCCGTCGCTTCTGCCTGCACGATTGATTGTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaag                                  >  1:3222963/1‑93 (MQ=255)

                                                         cGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaagca                                <  1:4291458/93‑1 (MQ=255)

                                                           tCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCaaa                              <  1:4945031/93‑1 (MQ=255)

                                                            cGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAAc                             <  1:5550360/93‑1 (MQ=255)

                                                            cGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAAc                             <  2:3458508/93‑1 (MQ=255)

                                                            cGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAAc                             <  1:1566466/93‑1 (MQ=255)

                                                              cTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCg                           <  1:6633047/93‑1 (MQ=255)

                                                               ttCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCgg                          >  2:1865924/1‑93 (MQ=255)

                                                               ttCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCgg                          >  2:2511093/1‑93 (MQ=255)

                                                               ttCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGGACGTCATCAGGGCGAGCAAACCgg                          >  1:63014/1‑93 (MQ=255)

                                                               ttCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGACGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCgg                          >  1:3530206/1‑93 (MQ=255)

                                                                tCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCggc                         >  1:1405310/1‑92 (MQ=255)

                                                                tCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgt                         <  2:6052053/93‑1 (MQ=255)

                                                                tCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgt                         <  1:4428578/93‑1 (MQ=255)

                                                                  tgcctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgtgt                       <  2:1296600/93‑1 (MQ=255)

                                                                  tgcctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgtgt                       <  2:2186800/93‑1 (MQ=255)

                                                                  tgcctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgtgt                       <  1:4141157/93‑1 (MQ=255)

                                                                  tgcctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgtgt                       <  1:5838770/93‑1 (MQ=255)

                                                                    cctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGttt                     <  2:916875/93‑1 (MQ=255)

                                                                    cctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGttt                     >  2:1953687/1‑93 (MQ=255)

                                                                    cctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGttt                     <  1:4736442/93‑1 (MQ=255)

                                                                     ctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTa                    >  2:956260/1‑93 (MQ=255)

                                                                     ctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTa                    <  1:2835283/93‑1 (MQ=255)

                                                                     ctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTa                    <  1:2955975/93‑1 (MQ=255)

                                                                     ctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTa                    <  2:6192383/93‑1 (MQ=255)

                                                                     ctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTa                    <  1:5678847/93‑1 (MQ=255)

                                                                      tgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTAc                   >  1:5897666/1‑93 (MQ=255)

                                                                      tgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTAc                   <  1:76734/93‑1 (MQ=255)

                                                                       gcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACg                  >  1:241284/1‑93 (MQ=255)

                                                                       gcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACg                  <  1:6514309/93‑1 (MQ=255)

                                                                        cACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGa                 <  2:6606057/93‑1 (MQ=255)

                                                                         aCGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGAc                >  2:5630566/1‑93 (MQ=255)

                                                                          cGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACa               <  1:5246041/93‑1 (MQ=255)

                                                                           gattgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACAt              >  2:6489421/1‑93 (MQ=255)

                                                                           gattgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACAt              >  2:3871373/1‑93 (MQ=255)

                                                                            attgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATc             >  2:2477282/1‑93 (MQ=255)

                                                                               gattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGaa          >  1:332100/1‑93 (MQ=255)

                                                                               gattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGaa          >  2:1453353/1‑93 (MQ=255)

                                                                                attTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGaaa         <  1:2366324/93‑1 (MQ=255)

                                                                                 ttttttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAAt        >  2:5555789/1‑93 (MQ=255)

                                                                                 ttttttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAAt        <  2:6646847/93‑1 (MQ=255)

                                                                                  tttttGCGCTCAGGCTTTGCCGGAGTCGCCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTCCGACATCGAAAtt       <  2:927747/93‑1 (MQ=255)

                                                                                    tttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGt     <  1:3238091/93‑1 (MQ=255)

                                                                                     ttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGtt    <  1:6383281/93‑1 (MQ=255)

ATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGTTAA  >  NC_000913/1458662‑1458841

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