BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,295,019	T→C	D254G (GAT→GGT)	yraQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,295,019	0	T	C	82.3%	492.6 / 104.1	215	D254G (GAT→GGT)	yraQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (20/18); new base C (89/88); total (109/106)
Fisher's exact test for biased strand distribution p-value = 8.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 223 total aligned reads displayed.

CAATCGGAATTTCTGCTGCCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGAAAAGCGTCCATA  >  NC_000913/3294927‑3295109

cAATCGGAATTTCTGCTGCCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCAt                                                                                            <  1:5988603/93‑1 (MQ=255)

  aTCGGAATTTCTGCTGCCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCg                                                                                          <  2:2958925/93‑1 (MQ=255)

  aTCGGAATTTCTGCTGCCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCg                                                                                          >  1:3396346/1‑93 (MQ=255)

  aTCGGAATTTCTGCTGCCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCg                                                                                          <  2:4801546/93‑1 (MQ=255)

      gAATTTCTGCTGCCGTTGGAATGACAAACAAGCATCCTGCTAGCGCCATCGCCACCACCCACATCAGTCTGTTATCGACAGCACCACCGGCAt                                                                                      >  1:3628859/1‑93 (MQ=255)

      gAATTTCTGCTGCCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCAt                                                                                      >  1:5621816/1‑93 (MQ=255)

         tttCTGCTGCCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATggg                                                                                   <  1:1775323/93‑1 (MQ=255)

           tCTGCTGCCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATggggg                                                                                 <  1:105206/93‑1 (MQ=255)

             tgctgcCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGaa                                                                               <  2:192422/93‑1 (MQ=255)

             tgctgcCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGaa                                                                               >  2:1555031/1‑93 (MQ=255)

              gctgcCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAAt                                                                              >  2:1622885/1‑93 (MQ=255)

               ctgcCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATa                                                                             >  1:2210247/1‑93 (MQ=255)

               ctgcCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATa                                                                             >  2:2697531/1‑93 (MQ=255)

                 gcCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGCATAAc                                                                           >  2:6189636/1‑93 (MQ=255)

                  ccGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAAcc                                                                          >  2:2281593/1‑93 (MQ=255)

                  ccGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAAcc                                                                          <  1:2914229/93‑1 (MQ=255)

                   cGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAAccg                                                                         >  1:6316693/1‑92 (MQ=255)

                    gTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAAACAg                                                                        >  2:3941055/1‑93 (MQ=255)

                    gTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAg                                                                        <  2:2593988/93‑1 (MQ=255)

                    gTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAg                                                                        <  1:6574239/93‑1 (MQ=255)

                     tGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGa                                                                       >  1:1714438/1‑93 (MQ=255)

                      gggAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAACCAGAc                                                                      <  1:2157180/93‑1 (MQ=255)

                      gggAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGAc                                                                      >  2:719641/1‑93 (MQ=255)

                       ggAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGAcg                                                                     >  1:6603705/1‑93 (MQ=255)

                        gAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGAcgc                                                                    <  2:3530180/93‑1 (MQ=255)

                         aaTGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGAcgcg                                                                   >  2:4347117/1‑93 (MQ=255)

                         aaTGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGAcgcg                                                                   <  2:696120/93‑1 (MQ=255)

                          aTGACAAACAAGCATCCTGCTTCCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGa                                                                  <  1:5267227/93‑1 (MQ=255)

                          aTGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGa                                                                  <  2:595038/93‑1 (MQ=255)

                          aTGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGa                                                                  <  1:1181881/93‑1 (MQ=255)

                           tGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAg                                                                 <  1:1239297/93‑1 (MQ=255)

                           tGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAg                                                                 <  1:3778439/93‑1 (MQ=255)

                            tacaaacaaGCATCGTGCTACCGCCATCGCCACCACCCACATCAGGCTGGTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGc                                                                <  2:1399319/92‑1 (MQ=255)

                             acaaacaaGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCg                                                               >  1:2021289/1‑93 (MQ=255)

                               aaacaaGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGc                                                             >  1:4785741/1‑93 (MQ=255)

                                aacaaGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAACCAGACGCGAGCGGCa                                                            >  1:5823816/1‑93 (MQ=255)

                                aacaaGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCa                                                            >  2:3149636/1‑93 (MQ=255)

                                   aaGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCAccc                                                         <  2:1661742/93‑1 (MQ=255)

                                        tCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGATAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAacac                                                    >  1:452852/1‑93 (MQ=255)

                                        tCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAacac                                                    >  2:5308130/1‑93 (MQ=255)

                                         ccTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACAcc                                                   <  2:1787979/93‑1 (MQ=255)

                                          cTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCa                                                  <  2:3699530/93‑1 (MQ=255)

                                          cTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGCCGCGAGCGGCACCCAACACCa                                                  <  2:3671758/93‑1 (MQ=255)

                                          cTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCa                                                  <  1:1603218/93‑1 (MQ=255)

                                          cTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCa                                                  >  2:2834539/1‑93 (MQ=255)

                                            gCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGt                                                >  2:1072359/1‑93 (MQ=255)

                                            gCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGt                                                >  1:1578527/1‑93 (MQ=255)

                                             cTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTa                                               <  1:675734/93‑1 (MQ=255)

                                               aCCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACt                                             >  2:5205412/1‑93 (MQ=255)

                                               aCCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACt                                             >  2:1224172/1‑93 (MQ=255)

                                                ccGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTg                                            >  2:4425991/1‑93 (MQ=255)

                                                ccGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTg                                            <  1:749255/93‑1 (MQ=255)

                                                 cGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGc                                           >  1:2993010/1‑93 (MQ=255)

                                                 cGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGc                                           <  2:4733807/93‑1 (MQ=255)

                                                 cGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGc                                           >  1:35604/1‑93 (MQ=255)

                                                 cGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGc                                           <  1:1469076/93‑1 (MQ=255)

                                                  cccATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGAAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCa                                          <  1:3695679/92‑1 (MQ=255)

                                                   ccATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCaa                                         >  1:455309/1‑93 (MQ=255)

                                                   ccATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCaa                                         <  2:5843520/93‑1 (MQ=255)

                                                    cATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAg                                        <  1:6100020/93‑1 (MQ=255)

                                                    cATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAg                                        >  1:2357105/1‑93 (MQ=255)

                                                    cATCGCCACCACCCACATCAGGCTGTTATCGACAGCACAACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAg                                        >  1:2857110/1‑93 (MQ=255)

                                                     aTCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAgg                                       <  1:5075767/93‑1 (MQ=255)

                                                      tCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGa                                      <  2:507487/93‑1 (MQ=255)

                                                      tCGCCACCACCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGa                                      <  2:2237019/93‑1 (MQ=255)

                                                        gCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATg                                    >  1:1084939/1‑93 (MQ=255)

                                                         cccccacccacATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGt                                   <  2:4046380/90‑1 (MQ=255)

                                                         ccaccaCCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGt                                   <  1:189716/93‑1 (MQ=255)

                                                         ccaccaCCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGt                                   >  1:5923378/1‑93 (MQ=255)

                                                          caccacCCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTa                                  >  2:4171652/1‑93 (MQ=255)

                                                           accaccCACATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTaa                                 >  1:5397272/1‑93 (MQ=255)

                                                             cacccacATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAAc                               >  2:4477674/1‑93 (MQ=255)

                                                             cacccacATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAAc                               <  2:5596029/93‑1 (MQ=255)

                                                               cccacATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCg                             >  2:1996681/1‑93 (MQ=255)

                                                               cccacATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCg                             >  1:944699/1‑93 (MQ=255)

                                                               cacacaTCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAACCCg                             <  1:5724013/91‑1 (MQ=255)

                                                                ccacATCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCgg                            >  2:4384865/1‑93 (MQ=255)

                                                                  acaTCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGa                          >  2:3368891/1‑93 (MQ=255)

                                                                  acaTCAGACTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACCCCAGGACTGCAAGGATGTAAACCgggg                          >  1:3841842/1‑92 (MQ=255)

                                                                   caTCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGAt                         <  2:6483864/93‑1 (MQ=255)

                                                                   caTCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGAt                         <  1:3455306/93‑1 (MQ=255)

                                                                   caTCAGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGAt                         <  1:5376170/93‑1 (MQ=255)

                                                                     tCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCg                       <  1:2745912/93‑1 (MQ=255)

                                                                       aGGCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTa                     >  1:5033505/1‑93 (MQ=255)

                                                                        ggCTGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTAc                    <  1:3283434/93‑1 (MQ=255)

                                                                           tGTTATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACAAGTACTGCAAGGATGTACACCGGGATCGTATTcc                 >  1:2282339/1‑93 (MQ=255)

                                                                            gTTATCGACAGCACCATCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCa                >  1:4338295/1‑93 (MQ=255)

                                                                              tATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGa              >  2:5769892/1‑93 (MQ=255)

                                                                              tATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGa              <  1:3677193/93‑1 (MQ=255)

                                                                              tATCGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGa              >  2:770733/1‑93 (MQ=255)

                                                                                 cGACAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGaaaa           >  2:2433461/1‑93 (MQ=255)

                                                                                   aCAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGAAAAGc         >  1:6286624/1‑93 (MQ=255)

                                                                                    cAGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGAAAAGCg        <  1:4950472/93‑1 (MQ=255)

                                                                                     aGCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGAAAAGCGt       >  1:3097873/1‑93 (MQ=255)

                                                                                      gCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGAAAAGCGTc      <  2:3876980/93‑1 (MQ=255)

                                                                                      gCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGAAAAGCGTc      >  2:4941574/1‑93 (MQ=255)

                                                                                      gCACCACCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAAAACCAGTACTGCAAGGCTGTAAACCGGGATCGTACTCCAGAAAAGCGTc      >  1:1012253/1‑93 (MQ=255)

                                                                                       caccacCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGAAAAGCGTcc     >  1:3574731/1‑93 (MQ=255)

                                                                                        aCCATCGGCATGGGGGAATACCCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGAAAAGCGTCCa    <  1:4704378/93‑1 (MQ=255)

                                                                                         ccACCGGCATAGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGAAAAGCGTCCAt   >  2:2930661/1‑93 (MQ=255)

CAATCGGAATTTCTGCTGCCGTGGGAATGACAAACAAGCATCCTGCTACCGCCATCGCCACCACCCACATCAGGCTGTTATCGACAGCACCATCGGCATGGGGGAATAACCAGACGCGAGCGGCACCCAACACCAGTACTGCAAGGATGTAAACCGGGATCGTACTCCAGAAAAGCGTCCATA  >  NC_000913/3294927‑3295109

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