BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,452,089	T→C	Q316R (CAG→CGG)	yheD ←	YheD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,452,089	0	T	C	82.4%	616.0 / 126.0	262	Q316R (CAG→CGG)	yheD	YheD
Reads supporting (aligned to +/- strand): ref base T (25/21); new base C (139/77); total (164/98)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 266 total aligned reads displayed.

TTAAGCGATGCGATCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATCACCACATGAG  >  NC_000913/3451999‑3452179

ttAAGCGATGCGATCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGa                                                                                          >  2:5207100/1‑93 (MQ=255)

   aGCGATGCGATCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCg                                                                                       <  1:5288908/93‑1 (MQ=255)

   aGCGATGCGATCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCg                                                                                       >  1:5510530/1‑93 (MQ=255)

   aGCGATGCGATCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCg                                                                                       <  2:1454504/93‑1 (MQ=255)

   aGCGATGCGATCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCg                                                                                       <  2:2710009/93‑1 (MQ=255)

    gcgatgcgatCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGc                                                                                      <  2:6398841/93‑1 (MQ=255)

      gatgcgatCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCAc                                                                                    >  2:1984956/1‑93 (MQ=255)

       atgcgatCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCAcc                                                                                   >  1:1935884/1‑93 (MQ=255)

       atgcgatCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCAcc                                                                                   <  2:6207308/93‑1 (MQ=255)

       atgcgatCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCAcc                                                                                   <  1:2663791/93‑1 (MQ=255)

        tgcgatCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCa                                                                                  <  2:4881817/93‑1 (MQ=255)

         gcgatCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAg                                                                                 <  2:6045428/93‑1 (MQ=255)

          cgatCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAgg                                                                                >  2:2346907/1‑93 (MQ=255)

          cgatCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAgg                                                                                >  2:5822732/1‑93 (MQ=255)

           gatCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGc                                                                               >  1:1472545/1‑93 (MQ=255)

             tCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCg                                                                             >  1:3452474/1‑93 (MQ=255)

               cAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGaa                                                                           >  1:4898033/1‑93 (MQ=255)

               cAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGaa                                                                           >  2:3986136/1‑93 (MQ=255)

               cAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGaa                                                                           >  1:5902590/1‑93 (MQ=255)

               cAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGaa                                                                           >  2:5309718/1‑93 (MQ=255)

                aGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAAt                                                                          <  1:5409510/93‑1 (MQ=255)

                 gggCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATc                                                                         >  1:2596922/1‑93 (MQ=255)

                 gggCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATc                                                                         <  2:6495287/93‑1 (MQ=255)

                  ggCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCa                                                                        <  1:423641/93‑1 (MQ=255)

                    cAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGc                                                                      >  2:54196/1‑93 (MQ=255)

                     aaaTTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCt                                                                     >  2:3081514/1‑93 (MQ=255)

                     aaaTTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCt                                                                     >  2:2616381/1‑93 (MQ=255)

                      aaTTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTg                                                                    <  2:2280096/93‑1 (MQ=255)

                      aaTTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTg                                                                    >  1:4993147/1‑93 (MQ=255)

                       aTTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTgg                                                                   <  1:1961312/93‑1 (MQ=255)

                        ttCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGCTggg                                                                  >  1:1256646/1‑93 (MQ=255)

                        ttCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTggg                                                                  >  1:788801/1‑93 (MQ=255)

                           tGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGAcc                                                               <  2:3495451/93‑1 (MQ=255)

                            gATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCt                                                              >  2:6545420/1‑93 (MQ=255)

                             aTCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTc                                                             >  2:4407462/1‑93 (MQ=255)

                              tCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCg                                                            >  1:1797373/1‑93 (MQ=255)

                              tCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCg                                                            >  2:1208494/1‑93 (MQ=255)

                              tCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCg                                                            <  1:2309907/93‑1 (MQ=255)

                               cTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGt                                                           >  2:3338241/1‑93 (MQ=255)

                                tACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTa                                                          >  2:4277743/1‑93 (MQ=255)

                                tACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTa                                                          >  2:2493677/1‑93 (MQ=255)

                                  ccAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTAcc                                                        >  1:6000522/1‑93 (MQ=255)

                                   cAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTAccc                                                       >  2:6348423/1‑93 (MQ=255)

                                    aaTGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCCTAcccc                                                      >  1:3782615/1‑93 (MQ=255)

                                     aTGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTAccccc                                                     <  1:324319/93‑1 (MQ=255)

                                     aTGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTAccccc                                                     <  1:74021/93‑1 (MQ=255)

                                      tGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCa                                                    <  2:4707228/93‑1 (MQ=255)

                                      tGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCa                                                    <  1:2965165/93‑1 (MQ=255)

                                      tGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCa                                                    >  2:1390068/1‑93 (MQ=255)

                                       gTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCAt                                                   <  1:3319000/93‑1 (MQ=255)

                                        tttGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATa                                                  >  1:351417/1‑93 (MQ=255)

                                              aagaagCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGt                                            >  1:44922/1‑93 (MQ=255)

                                                gaagCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGtat                                          <  1:6376364/93‑1 (MQ=255)

                                                 aagCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGACCGAATCAGCTGTGACCTCGTACCCCCATACGCCGcata                                         >  1:1184506/1‑93 (MQ=255)

                                                  agCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATaa                                        >  1:1045770/1‑93 (MQ=255)

                                                  agCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGACGTATaa                                        >  1:2202540/1‑93 (MQ=255)

                                                   gCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAac                                       >  1:2263711/1‑93 (MQ=255)

                                                    cATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAaca                                      <  2:426463/93‑1 (MQ=255)

                                                     aTTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAacac                                     <  1:2874753/93‑1 (MQ=255)

                                                      ttCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACg                                    >  2:5023979/1‑93 (MQ=255)

                                                      ttCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTTGGACCTCGTACCCCCATACGCCGTATAACACg                                    >  2:477754/1‑93 (MQ=255)

                                                      ttCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACg                                    >  2:5686928/1‑93 (MQ=255)

                                                       tCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGc                                   >  2:3138826/1‑93 (MQ=255)

                                                         ccACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCtt                                 >  1:272911/1‑93 (MQ=255)

                                                         ccACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCtt                                 >  2:897571/1‑93 (MQ=255)

                                                          cACTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTa                                <  2:2813391/93‑1 (MQ=255)

                                                           aCTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAg                               >  1:554507/1‑93 (MQ=255)

                                                           aCTTTTACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAg                               >  2:3599737/1‑93 (MQ=255)

                                                             ttttACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAgcg                             >  1:2710998/1‑93 (MQ=255)

                                                               ttACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGcc                           <  2:3568114/93‑1 (MQ=255)

                                                               ttACATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGcc                           <  1:567930/93‑1 (MQ=255)

                                                                 aCATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTc                         <  2:4171884/93‑1 (MQ=255)

                                                                 aCATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTc                         >  1:1414046/1‑93 (MQ=255)

                                                                 aCATGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTc                         >  2:1846054/1‑93 (MQ=255)

                                                                  cATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCa                        >  2:4450022/1‑93 (MQ=255)

                                                                   aTGCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTCGCGCCCCAc                       >  2:1626848/1‑93 (MQ=255)

                                                                     gCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTa                     >  2:2554843/1‑93 (MQ=255)

                                                                     gCCAGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTa                     <  1:4822408/93‑1 (MQ=255)

                                                                        aGCCCCGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATcg                  >  1:577485/1‑93 (MQ=255)

                                                                         gCCCCGCTGTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATcgc                 >  1:6499909/1‑93 (MQ=255)

                                                                          ccccGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATcgct                >  1:4280620/1‑92 (MQ=255)

                                                                           cccGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCTTATAACACGCTTAGCGCCTCACTATCgcgc               >  2:4072983/1‑93 (MQ=255)

                                                                           cccGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCgcgc               <  2:4021436/93‑1 (MQ=255)

                                                                           cccGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCgcgc               >  1:2003336/1‑93 (MQ=255)

                                                                            ccGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCa              <  2:3969861/93‑1 (MQ=255)

                                                                            ccGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCa              >  2:379565/1‑93 (MQ=255)

                                                                            ccGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCa              >  1:1490859/1‑93 (MQ=255)

                                                                             cGCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCAt             >  1:363245/1‑93 (MQ=255)

                                                                              gCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATc            >  1:3059425/1‑93 (MQ=255)

                                                                              gCTCTTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATc            <  2:4496220/93‑1 (MQ=255)

                                                                               ctctTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATca           >  2:2981590/1‑93 (MQ=255)

                                                                               ctctTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATca           >  2:6111162/1‑93 (MQ=255)

                                                                                 ctTACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATcacc         >  1:65809/1‑93 (MQ=255)

                                                                                  ttACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATcacca        <  1:907945/93‑1 (MQ=255)

                                                                                   tACTGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATcaccac       >  2:5902731/1‑93 (MQ=255)

                                                                                     cTGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATCACCACAt     >  2:769932/1‑93 (MQ=255)

                                                                                      tGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATCACCACATg    <  1:3812305/93‑1 (MQ=255)

                                                                                      tGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATCACCACATg    <  2:2680941/93‑1 (MQ=255)

                                                                                      tGCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATCACCACATg    >  1:467288/1‑93 (MQ=255)

                                                                                       gCCCGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATCACCACATGa   <  1:1338366/93‑1 (MQ=255)

TTAAGCGATGCGATCCAGGGCAAATTCTGATCTACCAATGTTTGTAAAGAAGCATTCCCACTTTTACATGCCAGCCCCGCTCTTACTGCCTGATCGCACCAGGCCGAATCAGCTGGGACCTCGTACCCCCATACGCCGTATAACACGCTTAGCGCCTCACTATCGCGCATCACCACATGAG  >  NC_000913/3451999‑3452179

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: `‑`