BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,616,789	C→T	S48S (TCG→TCA)	ydeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,616,789	0	C	T	98.0%	697.0 / ‑0.8	200	S48S (TCG→TCA)	ydeB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (3/1); new base T (86/110); total (89/111)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 207 total aligned reads displayed.

AATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACCGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGTTAGCTAA  >  NC_000913/1616699‑1616881

aaTTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGa                                                                                            <  2:6793961/93‑1 (MQ=255)

 aTTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAg                                                                                           >  1:5503554/1‑93 (MQ=255)

 aTTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAg                                                                                           >  1:6197377/1‑93 (MQ=255)

  ttCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAgg                                                                                          <  1:3622925/93‑1 (MQ=255)

  ttCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAgg                                                                                          >  1:2385507/1‑93 (MQ=255)

    cGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGcc                                                                                        >  2:1931294/1‑93 (MQ=255)

      cAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCAt                                                                                      >  1:7603480/1‑93 (MQ=255)

       aGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATc                                                                                     <  1:1539646/93‑1 (MQ=255)

         aCCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCaa                                                                                   <  1:7258550/93‑1 (MQ=255)

         aCCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCaa                                                                                   >  2:3621627/1‑93 (MQ=255)

         aCCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCaa                                                                                   <  2:3567273/93‑1 (MQ=255)

         aCCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACCGAGGCCATCaa                                                                                   <  1:5927604/93‑1 (MQ=255)

          ccGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAAc                                                                                  >  2:5150614/1‑93 (MQ=255)

            gggAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGa                                                                                >  1:1364097/1‑93 (MQ=255)

             ggAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGAc                                                                               <  2:1238209/93‑1 (MQ=255)

              gAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACt                                                                              <  1:5530807/93‑1 (MQ=255)

              gAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACt                                                                              <  2:4389593/93‑1 (MQ=255)

              gAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACt                                                                              <  1:5165303/93‑1 (MQ=255)

              gAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACt                                                                              <  2:384696/93‑1 (MQ=255)

               aaTCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTg                                                                             >  2:7203968/1‑93 (MQ=255)

               aaTCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTg                                                                             <  2:7377225/93‑1 (MQ=255)

                aTCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGa                                                                            <  2:7648058/93‑1 (MQ=255)

                 tCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGAc                                                                           >  1:7526568/1‑93 (MQ=255)

                 tCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGAc                                                                           >  2:2077234/1‑93 (MQ=255)

                   gAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGa                                                                         <  1:7863969/93‑1 (MQ=255)

                      aTGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGATGATTa                                                                      <  2:4777076/93‑1 (MQ=255)

                       tGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTAc                                                                     <  1:5396436/93‑1 (MQ=255)

                         ccAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGt                                                                   <  2:460481/93‑1 (MQ=255)

                          cAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGtt                                                                  <  2:4015545/93‑1 (MQ=255)

                           aaaTGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTc                                                                 >  2:1161569/1‑93 (MQ=255)

                            aaTTTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCg                                                                >  2:6380078/1‑93 (MQ=255)

                            aaTGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCg                                                                <  1:5718453/93‑1 (MQ=255)

                             aTGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCgg                                                               <  1:2084159/93‑1 (MQ=255)

                               gTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCa                                                             <  1:7102033/93‑1 (MQ=255)

                                  tCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGATGATTACGTTCGGCACTg                                                          >  1:3374908/1‑93 (MQ=255)

                                  tCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTg                                                          <  2:3461538/93‑1 (MQ=255)

                                  tCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTg                                                          >  1:3799996/1‑93 (MQ=255)

                                    cATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGtt                                                        <  1:2847571/93‑1 (MQ=255)

                                    cATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGtc                                                        >  2:1398223/1‑92 (MQ=255)

                                    cATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTCCGGCACTGta                                                        >  1:5995779/1‑92 (MQ=255)

                                     aTCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTc                                                       <  1:4235914/93‑1 (MQ=255)

                                     aTCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTc                                                       >  1:3896982/1‑93 (MQ=255)

                                      tCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCa                                                      <  2:4730033/93‑1 (MQ=255)

                                      tCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCa                                                      >  2:2740322/1‑93 (MQ=255)

                                      tCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCa                                                      <  1:1241168/93‑1 (MQ=255)

                                      tCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCa                                                      <  1:4132744/93‑1 (MQ=255)

                                      tCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTAAGTTCCGCACTGGTCa                                                      >  2:5023063/1‑93 (MQ=255)

                                          cAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGtt                                                  <  1:1408640/93‑1 (MQ=255)

                                           aGTTGCCCGGCGTAATACGCCACAATCAAGATGGCAAATACGCCCACTGAGGCCATCTACGACTGACGATTACGGTCGGCACTGTTCCTGTtg                                                 >  2:6973913/1‑93 (MQ=255)

                                            gTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCGCTGTTCATGTtgc                                                >  2:5739659/1‑93 (MQ=255)

                                            gTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTtgc                                                >  1:4639371/1‑93 (MQ=255)

                                             ttGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTtgcc                                               <  1:6692161/93‑1 (MQ=255)

                                             ttGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTtgcc                                               <  1:1474497/93‑1 (MQ=255)

                                             ttGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTtgcc                                               <  1:2042024/93‑1 (MQ=255)

                                               gCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTtgcctg                                             >  2:852001/1‑93 (MQ=255)

                                                  cGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCa                                          >  2:808235/1‑93 (MQ=255)

                                                  cGGCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCa                                          <  2:5714903/93‑1 (MQ=255)

                                                   ggCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAg                                         <  1:6909462/93‑1 (MQ=255)

                                                   ggCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAg                                         <  2:5992442/93‑1 (MQ=255)

                                                    gCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAgg                                        >  1:505928/1‑93 (MQ=255)

                                                    gCGTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAgg                                        >  2:5659329/1‑93 (MQ=255)

                                                      gTAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGcc                                      >  2:7081213/1‑93 (MQ=255)

                                                       tAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCa                                     >  2:135686/1‑93 (MQ=255)

                                                       tAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCa                                     >  2:7553707/1‑93 (MQ=255)

                                                       tAATACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCa                                     >  2:3738071/1‑93 (MQ=255)

                                                         aTACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAg                                   >  2:6284962/1‑93 (MQ=255)

                                                         aTACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAg                                   >  1:1659725/1‑93 (MQ=255)

                                                          tACGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGa                                  >  2:1375607/1‑93 (MQ=255)

                                                           aCGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGaa                                 <  2:413781/93‑1 (MQ=255)

                                                            cGCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGaaa                                <  2:3410157/93‑1 (MQ=255)

                                                             gCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAAt                               <  1:1042485/93‑1 (MQ=255)

                                                             gCCACCATCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACCATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAAt                               >  2:706485/1‑93 (MQ=255)

                                                              ccaccaTCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATa                              <  2:1497472/93‑1 (MQ=255)

                                                               caccaTCATGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATaa                             <  1:4041385/93‑1 (MQ=255)

                                                                    tcatGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGa                        >  2:3102077/1‑93 (MQ=255)

                                                                    tcatGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGa                        >  1:4393443/1‑93 (MQ=255)

                                                                     catGATGGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGAc                       <  2:4662005/93‑1 (MQ=255)

                                                                       tgatgGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAg                     >  1:1705286/1‑93 (MQ=255)

                                                                       tgatgGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAg                     >  2:6613280/1‑93 (MQ=255)

                                                                       tgatgGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAg                     >  2:2134011/1‑93 (MQ=255)

                                                                       tgatgGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAg                     <  1:5234111/93‑1 (MQ=255)

                                                                          tgGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAgttg                  >  2:332339/1‑93 (MQ=255)

                                                                          tgGCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAgttg                  <  1:4227323/93‑1 (MQ=255)

                                                                           ggCAAATACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAgttgt                 >  1:4273720/1‑93 (MQ=255)

                                                                              aaaTACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTaa              <  1:785865/93‑1 (MQ=255)

                                                                              aaaTACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTaa              >  1:962288/1‑93 (MQ=255)

                                                                              aaaTACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTaa              <  1:3290003/93‑1 (MQ=255)

                                                                                 tACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATgg           <  1:1738264/93‑1 (MQ=255)

                                                                                 tACGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATgg           >  2:2689107/1‑93 (MQ=255)

                                                                                   cGTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGt         >  2:7864489/1‑93 (MQ=255)

                                                                                   cGTACACTGAGGCCACCAACGACTGACGATTACGTTCGGCGCTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGt         <  1:565650/93‑1 (MQ=255)

                                                                                    gTACACTGAGGCCATCAACGGCTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGtt        <  1:2721320/93‑1 (MQ=255)

                                                                                    gTACACTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGtt        >  1:236385/1‑93 (MQ=255)

                                                                                      acacTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGTTAg      <  1:6309841/93‑1 (MQ=255)

                                                                                      acacTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGTTAg      <  1:3095904/93‑1 (MQ=255)

                                                                                       cacTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGTTAGc     <  2:2150221/93‑1 (MQ=255)

                                                                                       cacTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGTTAGc     >  1:5541529/1‑93 (MQ=255)

                                                                                       cacTGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGTTAGc     <  1:7666985/93‑1 (MQ=255)

                                                                                          cGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGTTAGCTaa  <  2:5433881/93‑1 (MQ=255)

                                                                                          cGAGGCCAACAAAGACTGATGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGTTAGCTaa  <  1:2660045/93‑1 (MQ=255)

AATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATACGTACACCGAGGCCATCAACGACTGACGATTACGTTCGGCACTGTTCATGTTGCCTGCCAGGCCAAGAAATAACGCGACAGTTGTTAATGGGTTAGCTAA  >  NC_000913/1616699‑1616881

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