BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,091,676	A→G	Q53R (CAG→CGG)	yqgF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,091,676	0	A	G	99.5%	673.7 / ‑4.8	208	Q53R (CAG→CGG)	yqgF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (85/122); total (86/122)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

Only 100 of 227 total aligned reads displayed.

TTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑AGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTA  >  NC_000913/3091586‑3091765

    cGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGa                                                                                     <  2:3479332/93‑1 (MQ=255)

    cGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGa                                                                                     <  1:1361911/93‑1 (MQ=255)

      gCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGGCGGACg                                                                                   >  2:1235072/1‑93 (MQ=255)

      gCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACg                                                                                   <  1:1422320/93‑1 (MQ=255)

       cACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGGCGa                                                                                  >  2:4471702/1‑93 (MQ=255)

         ccgcccgcCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGaaa                                                                                <  2:2148780/93‑1 (MQ=255)

         ccGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGaaa                                                                                <  1:5601844/93‑1 (MQ=255)

         ccGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGaaa                                                                                <  2:2270622/93‑1 (MQ=255)

          cGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAat                                                                               <  1:2739169/93‑1 (MQ=255)

           gCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGGCCGGCGAAata                                                                              >  1:283228/1‑92 (MQ=255)

           gCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGACGGGCAAAatc                                                                              >  1:5560822/1‑93 (MQ=255)

           ccTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAatc                                                                              <  2:6129512/92‑1 (MQ=255)

             tCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGGCGGACGAAatcat                                                                            >  2:4252966/1‑93 (MQ=255)

             tCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAatcat                                                                            <  1:6325596/93‑1 (MQ=255)

              cGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAatcatc                                                                           <  1:4958637/93‑1 (MQ=255)

               gCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGGCGGACGAAATCAtcg                                                                          >  1:6544330/1‑93 (MQ=255)

                cccTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCAtcgt                                                                         <  1:2731811/93‑1 (MQ=255)

                cccTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCAtcgt                                                                         <  2:2277889/93‑1 (MQ=255)

                 ccTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCAtcgtc                                                                        <  2:5355765/93‑1 (MQ=255)

                   tttGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCgg                                                                      <  2:4066350/93‑1 (MQ=255)

                    ttGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGt                                                                     <  2:6562781/93‑1 (MQ=255)

                     tGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGtt                                                                    <  2:1377463/93‑1 (MQ=255)

                     tGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGtt                                                                    <  1:4783742/93‑1 (MQ=255)

                     tGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGtt                                                                    <  1:6515105/93‑1 (MQ=255)

                     tGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGtt                                                                    <  1:1228481/93‑1 (MQ=255)

                     tGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGtt                                                                    <  2:4978333/93‑1 (MQ=255)

                     tGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGtt                                                                    <  2:3704399/93‑1 (MQ=255)

                     tGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGtt                                                                    <  2:3748220/93‑1 (MQ=255)

                       ccTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGGCGGACGAAAACATCGTCGGGTTg                                                                  >  1:2804141/1‑93 (MQ=255)

                       ccTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGGCCGGCGAAAACATCGGCGGGTTg                                                                  >  2:3642577/1‑93 (MQ=255)

                       ccTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTg                                                                  <  2:1651967/93‑1 (MQ=255)

                        cTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGc                                                                 <  1:5977461/93‑1 (MQ=255)

                         tGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAGAGAGTGGC‑GGCCGGACGAAATAATCGTCGGTTTGcc                                                                >  1:2156733/1‑93 (MQ=255)

                         tGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGGCGGGAGAAATCATCGTCGGGTTggc                                                                >  2:2439258/1‑91 (MQ=255)

                         tGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGGAATTCTCCTCGGGTTGcc                                                                >  1:3304161/1‑93 (MQ=255)

                         tGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGcc                                                                <  2:3194270/93‑1 (MQ=255)

                          gCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGGCGGGCGAAAACATCGTCGGGTTGGCg                                                               >  1:4255671/1‑93 (MQ=255)

                           cAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGc                                                              <  1:4988306/93‑1 (MQ=255)

                           cAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGc                                                              <  1:4476492/93‑1 (MQ=255)

                            aaTTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCCGGACG‑AATCATTGTAGGGTTGTCGCt                                                             >  1:4308552/1‑93 (MQ=255)

                            aaTTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAAACATCGTCGGGTTGTCGCt                                                             >  2:5738937/1‑93 (MQ=255)

                               tAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGaa                                                          <  2:781447/93‑1 (MQ=255)

                                aaaGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGGCGGGCGAAATCATCGTCGGGTTGGCGCTGAat                                                         >  1:3967347/1‑93 (MQ=255)

                                  aGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGGCGAAAACATCGTCGGGTTGCCGCTGAatat                                                       >  2:2419425/1‑93 (MQ=255)

                                   gCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATg                                                      <  2:2603311/93‑1 (MQ=255)

                                   gCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATg                                                      <  2:3321174/93‑1 (MQ=255)

                                    cacaGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATgg                                                     <  2:4798869/93‑1 (MQ=255)

                                     acaGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGGCGAAATCATCGTCGGGTTGTCCCTGAATATGGa                                                    >  2:3401473/1‑93 (MQ=255)

                                     acaGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGa                                                    <  2:4280557/93‑1 (MQ=255)

                                      caGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGGCGAAAACATCGTCGGGTTGGCGCTGAATATgggg                                                   >  1:255403/1‑91 (MQ=255)

                                        ggACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGGTTGCCGCTGGATATgggagg                                                 >  1:2303277/1‑89 (MQ=255)

                                         gACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGGCGAAATCATCCTCGGTTTGTCGCTGAATATGggcggc                                                >  2:3862500/1‑93 (MQ=255)

                                           cGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCAc                                              <  2:3595371/93‑1 (MQ=255)

                                             gTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCg                                            <  1:4046675/93‑1 (MQ=255)

                                              tACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGGCCACAACATCGTTGGCTTGCCGCCGAATATGGGCGGCAcccg                                           >  1:1673499/1‑91 (MQ=255)

                                              tACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGa                                           <  1:1307551/93‑1 (MQ=255)

                                               aCGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATTATCGTCGGTTTGCCCCCCAATATGGAAGGGACCGac                                          >  2:3006222/1‑92 (MQ=255)

                                               aCGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAg                                          <  2:2732094/93‑1 (MQ=255)

                                                cGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGc                                         <  1:4699093/93‑1 (MQ=255)

                                                   cGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAAc                                      <  2:72382/93‑1 (MQ=255)

                                                      aCTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGGCGAAATCATCGTCGGGTTGCCGCTGAATATGGACGGCACCCAGCAACCAt                                   >  1:328141/1‑93 (MQ=255)

                                                       cTNGAANATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCAtt                                  <  2:1719129/93‑1 (MQ=255)

                                                       cTGTAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCAtt                                  <  2:5106718/93‑1 (MQ=255)

                                                       cTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGGCGAAATCATCGTCGGGTTTCCGCTGAATATGGACGGGACCCAGCAACCAtt                                  >  1:1919284/1‑93 (MQ=255)

                                                       cTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCCTCGGGTTGCCGCTGAATATGGACGGCACCGCGCAAACact                                  >  2:3514096/1‑91 (MQ=255)

                                                        tGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGGTTGCCGCTGAATATGGACGGGAACGAGCAAACAttt                                 >  2:5443302/1‑92 (MQ=255)

                                                        tGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCCGGCGAAATCATCGTCGGGTTGTCGCTGAATATGGGCGGGAACGCGCAACCATTg                                 >  1:218473/1‑93 (MQ=255)

                                                         ggAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGGCGGGCGAAAACATCGTCCGGTTGCCGCTGGATATTGACGGGGACGAGGAAACATTGa                                >  1:3221999/1‑93 (MQ=255)

                                                         ggAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGa                                <  2:1478901/93‑1 (MQ=255)

                                                         ggAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGGTTGCCGCTGTATATGGACGGCACCGAGGAAACATTgc                                >  2:5241325/1‑92 (MQ=255)

                                                          gAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGAc                               <  2:6333060/93‑1 (MQ=255)

                                                          gAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGAc                               <  1:1755342/93‑1 (MQ=255)

                                                          gAACATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGAc                               <  2:6686380/93‑1 (MQ=255)

                                                            aCATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTg                             <  1:3482504/93‑1 (MQ=255)

                                                            aCATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCAACGTCGGGTTGCCGCTGAATATGGACGGCACCGGGGAACCACTggatg                             >  2:4420389/1‑89 (MQ=255)

                                                             cATTATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGGCGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGGACCGCGGAACCATTTACttc                            >  1:3063207/1‑91 (MQ=255)

                                                              aTTATCGAGCGTTTACTGAAAGAGTGGC‑GGGCGGACGAAATCATCGTCGGTTTGCCGCTGAATATTGAACGCACCGAGCAAACATTGACTggc                           >  1:777772/1‑91 (MQ=255)

                                                               ttATCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAAACATCGTCGGTTTGCCGCCGGATATGGACGGCACCGAGCAACCATTGACTGCCa                          >  2:6431481/1‑93 (MQ=255)

                                                                 aTCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGGCGAAATCATCGTCGGTTTGCCGCTGAATATTGGCGGCCCCGAGCAACCATTTACTGCCaga                        >  2:4917854/1‑93 (MQ=255)

                                                                 aTCGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCaga                        <  2:1999616/93‑1 (MQ=255)

                                                                   cGAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAgc                      <  2:1430776/93‑1 (MQ=255)

                                                                    gAGCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAgcg                     <  1:4273621/93‑1 (MQ=255)

                                                                      gCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGcgcg                   <  2:1560759/93‑1 (MQ=255)

                                                                      gCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGcgcg                   <  2:4228508/93‑1 (MQ=255)

                                                                      gCGTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGGTTGCCGCTGAATATGGACGGCACCGAGCAAACATTGACTGCCAGAGcgcg                   >  1:1002527/1‑93 (MQ=255)

                                                                        gTTTACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTa                 <  1:5010341/93‑1 (MQ=255)

                                                                           tACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCAGCGAGCAACCATGGACTGCCAGAGCGCGTAAAt              <  1:6436447/93‑1 (MQ=255)

                                                                           tACTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAAt              <  2:6148580/93‑1 (MQ=255)

                                                                            aCTGAAAGAGTGGC‑GGCCGGGCGAAATCATCGTCGGGTTGCCGCTGAATATGGGCGGGACCGAGCAACCATTGACTGCCAGAGCGCGGAAAtt             >  2:2123454/1‑93 (MQ=255)

                                                                             cTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAAttt            <  1:1096721/93‑1 (MQ=255)

                                                                             cTGAAAGAGTGGC‑GGCCGGACGAAATCATCGTCGGGTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCCGAGCGCGTAAAttt            >  2:4499861/1‑93 (MQ=255)

                                                                              tGAAAGAGTGGC‑GGCCGGACGAAATCATCCTCGGTTTGCCGCTGATTATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTg           >  1:2364099/1‑93 (MQ=255)

                                                                                aaaGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGcc         <  1:2046723/93‑1 (MQ=255)

                                                                                 aaGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTTGATATGGACGGCACCGAGCAACCATTGACTTGCAGAGCGCGTAAAATTTCCa        >  2:1799735/1‑93 (MQ=255)

                                                                                 aaGAGTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTTACTGCCAGAGCGCGTAAATTTGCCa        >  2:5936634/1‑93 (MQ=255)

                                                                                  agagTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGGAACCATTTACTGCCAGAGCGCGTAAATTTGCCaa       >  2:4658106/1‑93 (MQ=255)

                                                                                  agagTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCaa       >  2:3219010/1‑93 (MQ=255)

                                                                                  agagTGGC‑GGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCaa       <  1:38862/93‑1 (MQ=255)

                                                                                       ggc‑ggcCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTa  <  2:439170/93‑1 (MQ=255)

                                                                                       ggc‑ggcCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTa  >  1:4172438/1‑93 (MQ=255)

TTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGC‑AGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTA  >  NC_000913/3091586‑3091765

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