BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,860,862	T→C	A123A (GCA→GCG)	b4553 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,860,862	0	T	C	100.0%	689.0 / NA	219	A123A (GCA→GCG)	b4553	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (118/101); total (118/101)

Only 100 of 228 total aligned reads displayed.

AGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCTGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGAAAC  >  NC_000913/3860773‑3860950

aGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGcc                                                                                       >  2:2889458/1‑93 (MQ=255)

 gcgcCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCa                                                                                      >  1:4203278/1‑93 (MQ=255)

 gcgcCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCCCTACCCGccc                                                                                      >  2:169741/1‑92 (MQ=255)

  ctccGTCACAATGCCGGAAATGATAATCGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCCCCAc                                                                                     <  2:4457111/91‑1 (MQ=255)

     cGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGGTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACcgc                                                                                  >  1:4612150/1‑93 (MQ=255)

       tCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGt                                                                                >  1:2229853/1‑93 (MQ=255)

         acaATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTcc                                                                              >  1:4285016/1‑93 (MQ=255)

         acaATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTcc                                                                              >  1:3209489/1‑93 (MQ=255)

         acaATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTcc                                                                              >  1:3024897/1‑93 (MQ=255)

          caATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCt                                                                             <  1:1105940/93‑1 (MQ=255)

             tGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAg                                                                          <  1:3953304/93‑1 (MQ=255)

               ccGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTg                                                                        >  2:1007881/1‑93 (MQ=255)

               ccGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTg                                                                        >  2:3979537/1‑93 (MQ=255)

                cGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGa                                                                       >  1:1087466/1‑93 (MQ=255)

                  gAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGaaa                                                                     <  1:3068601/93‑1 (MQ=255)

                   aaaTGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAAt                                                                    >  1:4648396/1‑93 (MQ=255)

                   aaaTGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAAt                                                                    <  2:1261203/93‑1 (MQ=255)

                   aaaTGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAAt                                                                    >  2:2200991/1‑93 (MQ=255)

                    aaTGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATc                                                                   <  2:541090/93‑1 (MQ=255)

                     aTGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCg                                                                  >  2:2006441/1‑93 (MQ=255)

                     aTGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCg                                                                  >  2:855520/1‑93 (MQ=255)

                     aTGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCg                                                                  >  2:4071454/1‑93 (MQ=255)

                      tGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGa                                                                 <  2:3876207/93‑1 (MQ=255)

                       gATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGag                                                                >  1:2093280/1‑92 (MQ=255)

                         tAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGc                                                              >  1:2171100/1‑93 (MQ=255)

                         tAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGc                                                              >  1:1368268/1‑93 (MQ=255)

                           aTTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCg                                                            >  2:2960870/1‑93 (MQ=255)

                            ttGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGa                                                           >  1:557793/1‑93 (MQ=255)

                            ttGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGa                                                           >  2:149807/1‑93 (MQ=255)

                            atGCGCCGATAATTCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGGTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGa                                                           <  2:3151987/92‑1 (MQ=255)

                             tGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGaa                                                          <  1:826912/93‑1 (MQ=255)

                             tGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGaa                                                          >  1:1999875/1‑93 (MQ=255)

                                gcCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTa                                                       <  2:3310999/93‑1 (MQ=255)

                                  cGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGc                                                     <  1:795632/93‑1 (MQ=255)

                                  cGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGc                                                     >  2:3392296/1‑93 (MQ=255)

                                   gATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCt                                                    >  1:1678100/1‑93 (MQ=255)

                                   gATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCt                                                    <  1:3228063/93‑1 (MQ=255)

                                    aTAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCtt                                                   >  2:970985/1‑93 (MQ=255)

                                      aaTGCTGGTATCGAGGGGTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTcc                                                 >  2:2481637/1‑93 (MQ=255)

                                       aTGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTccc                                                <  1:3181250/93‑1 (MQ=255)

                                        tGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTcccc                                               <  2:2712544/93‑1 (MQ=255)

                                          cTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAg                                             <  1:3747941/93‑1 (MQ=255)

                                           tGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAgt                                            >  1:3909046/1‑92 (MQ=255)

                                           tGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAgg                                            <  2:2781215/93‑1 (MQ=255)

                                           tGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAgg                                            >  1:430998/1‑93 (MQ=255)

                                            ggTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGt                                           >  1:1756870/1‑93 (MQ=255)

                                            ggTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGt                                           <  1:2807423/93‑1 (MQ=255)

                                              tATCGAGGGGTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCa                                         >  2:1124275/1‑93 (MQ=255)

                                                tCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCAtt                                       >  1:2057053/1‑93 (MQ=255)

                                                tCGAGGGGTTTAATACCGGCCATCATTGACAGAACGCTACTCCCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTTGTTTCCCCAGGTCAtt                                       >  2:1595137/1‑93 (MQ=255)

                                                  gAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTcc                                     <  2:3816718/93‑1 (MQ=255)

                                                   aGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTccc                                    >  1:3371879/1‑93 (MQ=255)

                                                   aGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTccc                                    >  2:3249917/1‑93 (MQ=255)

                                                    gggTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCa                                   >  2:4563684/1‑93 (MQ=255)

                                                     ggTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCAt                                  >  1:1035651/1‑93 (MQ=255)

                                                      gTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATc                                 >  1:4069605/1‑93 (MQ=255)

                                                      gTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATc                                 <  1:3631979/93‑1 (MQ=255)

                                                      gTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATc                                 >  2:1178942/1‑93 (MQ=255)

                                                         ttAATCCCGGCTATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATcgcg                              <  2:508579/93‑1 (MQ=255)

                                                          tAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGt                             >  2:153214/1‑93 (MQ=255)

                                                          tAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGt                             <  1:1565482/93‑1 (MQ=255)

                                                            aTCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtg                           >  2:706066/1‑93 (MQ=255)

                                                             tCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtga                          <  1:2362916/93‑1 (MQ=255)

                                                               ccGGCCACTATTCTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtgatg                        <  1:339369/93‑1 (MQ=255)

                                                                cGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtgatga                       <  2:127904/93‑1 (MQ=255)

                                                                cGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtgatga                       >  2:2525156/1‑93 (MQ=255)

                                                                 ggCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGaa                      <  2:235967/93‑1 (MQ=255)

                                                                   ccATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAAt                    <  1:1657312/93‑1 (MQ=255)

                                                                    catcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATa                   <  2:1415462/93‑1 (MQ=255)

                                                                    catcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATa                   >  1:895898/1‑93 (MQ=255)

                                                                     atcatTGTCAGACCGCTACCCTCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAg                  <  1:4528657/93‑1 (MQ=255)

                                                                     atcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAg                  >  2:1307577/1‑93 (MQ=255)

                                                                     atcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAg                  >  1:1434218/1‑93 (MQ=255)

                                                                     atcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAg                  <  2:4275606/93‑1 (MQ=255)

                                                                      tcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGt                 <  2:2940359/93‑1 (MQ=255)

                                                                      tcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGt                 >  1:3776626/1‑93 (MQ=255)

                                                                      tcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGt                 <  1:2915181/93‑1 (MQ=255)

                                                                      tcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGGCATTCCCATCGCGTTGATGAAATAGt                 >  1:2572615/1‑93 (MQ=255)

                                                                       catTGTCAGGCCGCTCCCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGtt                <  1:1011016/93‑1 (MQ=255)

                                                                       catTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGtt                <  1:1238375/93‑1 (MQ=255)

                                                                       catTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGtt                >  1:1543847/1‑93 (MQ=255)

                                                                         ttGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTcc              >  2:3641985/1‑93 (MQ=255)

                                                                          tGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCa             <  1:3004692/93‑1 (MQ=255)

                                                                          tGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCa             <  1:4532564/93‑1 (MQ=255)

                                                                          tGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCa             <  1:3510406/93‑1 (MQ=255)

                                                                          tGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCa             <  1:246559/93‑1 (MQ=255)

                                                                          tGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCa             <  2:260992/93‑1 (MQ=255)

                                                                            tCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAgg           >  2:2110037/1‑93 (MQ=255)

                                                                            tCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAgg           <  2:2998508/93‑1 (MQ=255)

                                                                            tCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAgg           <  1:4084064/93‑1 (MQ=255)

                                                                             cAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGt          <  1:1171662/93‑1 (MQ=255)

                                                                               gACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCGGGCCATTCCCATCGCGTTGATGAAAAAGTTCCAGGtag        >  2:4618927/1‑91 (MQ=255)

                                                                               gACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGGTGATGAAATAGTTCAAGGTcc        >  2:3228297/1‑92 (MQ=255)

                                                                                aCCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAg       <  1:3280246/93‑1 (MQ=255)

                                                                                 ccGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAgg      >  2:4007614/1‑93 (MQ=255)

                                                                                  cGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGa     <  1:4676729/93‑1 (MQ=255)

                                                                                   gCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGaa    >  2:2768144/1‑93 (MQ=255)

                                                                                     tACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGAAAc  >  1:108063/1‑93 (MQ=255)

                                                                                     tACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGAAAc  <  1:1823659/93‑1 (MQ=255)

                                                                                     caCCCGCCACCTCGTCCTGAGTAAACTCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGAAAc  <  2:3663577/92‑1 (MQ=255)

AGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCTGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGAAAC  >  NC_000913/3860773‑3860950

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