BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,215,450	T→C	Q423R (CAG→CGG)	aer ←	Aer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,215,450	0	T	C	99.5%	591.7 / NA	189	Q423R (CAG→CGG)	aer	Aer
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (80/108); total (81/108)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 209 total aligned reads displayed.

GTTGAATGGCTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCTGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGC  >  NC_000913/3215360‑3215539

gTTGAATGGCTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCtgccgc                                                                                         >  1:1571391/1‑93 (MQ=255)

gTTGAATGGCTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCtgccgc                                                                                         <  1:2352547/93‑1 (MQ=255)

gTTGAATGGCTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCtgccgc                                                                                         >  1:1825982/1‑93 (MQ=255)

gTTGAATGGCTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCtgccgc                                                                                         <  1:1618593/93‑1 (MQ=255)

gTTGAATGGCTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCtgccgc                                                                                         >  1:392279/1‑93 (MQ=255)

    aaTGGCTAGTCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAAc                                                                                     <  2:3880161/93‑1 (MQ=255)

      tGGCTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCg                                                                                   <  1:4311443/93‑1 (MQ=255)

        gCTAATCTGTGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGa                                                                                 <  2:526754/93‑1 (MQ=255)

        gCTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGa                                                                                 <  2:4588046/93‑1 (MQ=255)

        gCTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGa                                                                                 >  2:3326152/1‑93 (MQ=255)

        gCTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGa                                                                                 <  1:2955917/93‑1 (MQ=255)

         cTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGAt                                                                                <  1:314900/93‑1 (MQ=255)

         cTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGAt                                                                                <  1:3803357/93‑1 (MQ=255)

          tAATCTGGGCGATCAACTGGGTGACGTTTTTCCCCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGAtt                                                                               <  1:1145835/93‑1 (MQ=255)

          tAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGAtt                                                                               >  1:1861028/1‑93 (MQ=255)

          tAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGAtt                                                                               >  2:299338/1‑93 (MQ=255)

          tAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGCCCGGCGGCGTGTACCTGCCGCGAACCGGAtt                                                                               <  2:2007288/93‑1 (MQ=255)

            aTCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGc                                                                             <  1:4652346/93‑1 (MQ=255)

            aTCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGc                                                                             >  1:2072448/1‑93 (MQ=255)

             tCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCa                                                                            <  1:1713690/93‑1 (MQ=255)

              cTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCCGATTGCAc                                                                           >  2:1923693/1‑93 (MQ=255)

                gggCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCt                                                                         <  1:3596121/93‑1 (MQ=255)

                 ggCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCtt                                                                        >  2:2589368/1‑93 (MQ=255)

                   cGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTAt                                                                      <  1:469208/93‑1 (MQ=255)

                    gATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATc                                                                     <  1:58512/93‑1 (MQ=255)

                    gATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCACCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATc                                                                     <  2:1902248/93‑1 (MQ=255)

                     aTCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCa                                                                    >  1:965535/1‑93 (MQ=255)

                       cAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGc                                                                  >  2:1097703/1‑93 (MQ=255)

                        aaCTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTCCCTGCCGCGAACCGGATTGCACCTTATCAGCa                                                                 <  2:2056512/93‑1 (MQ=255)

                        aaCTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCa                                                                 >  1:1008845/1‑93 (MQ=255)

                         aCTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCAc                                                                >  1:2743108/1‑93 (MQ=255)

                         aCTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCAc                                                                >  1:1107089/1‑93 (MQ=255)

                          cTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACt                                                               <  1:3247271/93‑1 (MQ=255)

                           tGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATTAGCACTg                                                              >  2:3392891/1‑93 (MQ=255)

                           tGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTg                                                              >  1:4114586/1‑93 (MQ=255)

                            gggTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTgg                                                             >  1:2227725/1‑93 (MQ=255)

                              gTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCa                                                           <  1:4389572/93‑1 (MQ=255)

                                gACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCTGCGAACCGGATTGCACCTTATCAGTACTGCCatc                                                         >  1:4012974/1‑93 (MQ=255)

                                gACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCatc                                                         >  1:1065941/1‑93 (MQ=255)

                                 aCGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCatca                                                        >  2:3705656/1‑93 (MQ=255)

                                 aCGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCatca                                                        <  1:4141299/93‑1 (MQ=255)

                                  cGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCatcaa                                                       >  2:731820/1‑93 (MQ=255)

                                    tttttCCCCTGTGCCACAATATCTTCCATCGTCCGCCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCatcaatc                                                     <  2:201979/93‑1 (MQ=255)

                                    tttttCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCatcaatc                                                     >  2:4113169/1‑93 (MQ=255)

                                    tttttCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCatcaatc                                                     <  1:4474351/93‑1 (MQ=255)

                                    tttttCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCatcaatc                                                     >  2:3175939/1‑93 (MQ=255)

                                     ttttCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCatcaatca                                                    <  2:1393027/93‑1 (MQ=255)

                                     ttttCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCatcaatca                                                    <  1:1301140/93‑1 (MQ=255)

                                       ttCTCCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGc                                                  <  1:2799560/93‑1 (MQ=255)

                                       ttCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGc                                                  <  2:2042592/93‑1 (MQ=255)

                                        tCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCt                                                 <  1:4613329/93‑1 (MQ=255)

                                         cgccTGTGCCAAAATATCATCCATCGTCCGTCGGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCtt                                                <  2:1198317/91‑1 (MQ=255)

                                          aCCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTg                                               >  2:2429444/1‑93 (MQ=255)

                                          aCCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTg                                               >  1:2858039/1‑93 (MQ=255)

                                            cTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCg                                             >  2:4022411/1‑93 (MQ=255)

                                             tgtgCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGa                                            >  1:2283749/1‑93 (MQ=255)

                                              gtgCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGaa                                           <  2:70431/93‑1 (MQ=255)

                                                 ccACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAatat                                        <  1:1212524/93‑1 (MQ=255)

                                                  cacaATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATc                                       >  2:1450628/1‑93 (MQ=255)

                                                  cacaATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATc                                       >  2:4468936/1‑93 (MQ=255)

                                                  cacaATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATc                                       <  1:2140501/93‑1 (MQ=255)

                                                  cacaATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATc                                       >  1:4358804/1‑93 (MQ=255)

                                                   acaATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCg                                      >  2:227991/1‑93 (MQ=255)

                                                    caATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGt                                     <  2:897230/93‑1 (MQ=255)

                                                    caATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGt                                     <  1:2943255/93‑1 (MQ=255)

                                                     aaTATTTTCCATCGTCCGTCCGGCGGCTTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGtt                                    <  1:3794093/93‑1 (MQ=255)

                                                     aaTATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGtt                                    >  1:355846/1‑93 (MQ=255)

                                                     aaTATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGtt                                    <  1:2982687/93‑1 (MQ=255)

                                                      atatCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTg                                   >  1:4456508/1‑93 (MQ=255)

                                                      atatCTTCCATCGTCCGTCCGGCGGCGTCTCCCTGCCGCGAACCGGATTGCCCCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTg                                   <  2:4342257/93‑1 (MQ=255)

                                                       tatCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTgg                                  >  1:98684/1‑93 (MQ=255)

                                                       tatCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTgg                                  <  1:2029115/93‑1 (MQ=255)

                                                       tatCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTgg                                  <  1:3161809/93‑1 (MQ=255)

                                                         tCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGgca                                <  2:1529178/93‑1 (MQ=255)

                                                            tCCATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGgcagca                             <  2:3517127/93‑1 (MQ=255)

                                                             ccATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCAt                            >  2:2336620/1‑93 (MQ=255)

                                                             ccATCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTCCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCAt                            <  2:2666520/93‑1 (MQ=255)

                                                             cagtCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCAt                            <  1:2609487/90‑1 (MQ=255)

                                                               aTCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGCATATCGTTGGCAGCCTTa                          >  2:669679/1‑93 (MQ=255)

                                                                tCGTCCGTCCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAg                         <  1:423542/93‑1 (MQ=255)

                                                                ccgtccgtcCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAg                         <  2:1123285/92‑1 (MQ=255)

                                                                ccgcccgtccgGCGGCGTGTCCCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAg                         <  2:4015112/89‑1 (MQ=255)

                                                                 cgtccgtcCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAgc                        <  1:3058548/93‑1 (MQ=255)

                                                                 cgtccgtcCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAgc                        >  1:1230130/1‑93 (MQ=255)

                                                                 cgtccgtcCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAgc                        <  1:2507500/93‑1 (MQ=255)

                                                                 cgtccgtcCGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAgc                        <  2:2888199/93‑1 (MQ=255)

                                                                  gtccgtccGGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAgcg                       >  2:859827/1‑93 (MQ=255)

                                                                   tccgtccgGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAgcgc                      >  1:3383101/1‑93 (MQ=255)

                                                                   tccgtccgGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAgcgc                      >  1:3641698/1‑93 (MQ=255)

                                                                    ccgtccgGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCt                     <  2:2976405/93‑1 (MQ=255)

                                                                      gtccgGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTgg                   >  2:3809664/1‑92 (MQ=255)

                                                                      gtccgGCGGCGTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGc                   >  1:740299/1‑93 (MQ=255)

                                                                             ggcgTGTACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCggctggc            >  2:4283879/1‑93 (MQ=255)

                                                                                 tgtACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTaaa        >  1:913295/1‑93 (MQ=255)

                                                                                   tACCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATg      <  1:1166390/93‑1 (MQ=255)

                                                                                    aCCTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGa     <  1:165959/93‑1 (MQ=255)

                                                                                     ccTGCCGCGAACCGGATTGCTCCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGAc    <  1:2317202/93‑1 (MQ=255)

                                                                                     ccTGCCGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGAc    <  2:3242199/93‑1 (MQ=255)

                                                                                       tgccgcGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGc  >  2:3921271/1‑93 (MQ=255)

                                                                                       tgccgcGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGc  <  2:606540/93‑1 (MQ=255)

GTTGAATGGCTAATCTGGGCGATCAACTGGGTGACGTTTTTCACCTGTGCCACAATATCTTCCATCGTCCGTCCGGCGGCGTGTACCTGCTGCGAACCGGATTGCACCTTATCAGCACTGGCATCAATCAGCTTGCGAATATCGTTGGCAGCATTAGCGCTGCGGCTGGCTAAATGACGC  >  NC_000913/3215360‑3215539

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