BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,817,156	T→C	S550G (AGT→GGT)	yicF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,817,156	0	T	C	97.5%	487.9 / 0.0	157	S550G (AGT→GGT)	yicF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/2); new base C (67/86); total (69/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 160 total aligned reads displayed.

ACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACT  >  NC_000913/3817067‑3817245

aCGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGcc                                                                                        >  1:2546409/1‑93 (MQ=255)

aCGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGcc                                                                                        >  1:1906052/1‑93 (MQ=255)

 cGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGccc                                                                                       <  2:4361682/93‑1 (MQ=255)

  gCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCa                                                                                      >  2:2922654/1‑93 (MQ=255)

  gCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCa                                                                                      >  1:1531410/1‑93 (MQ=255)

     ttAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAACCGTGTGATCTGCTGGGCAGCCAGCGAACGTCCCAGct                                                                                   >  1:4333839/1‑93 (MQ=255)

     ttAAGCTGGGATTACCGGTATTTTACTACAGCCACGAAGGACACTAAGGGTCAAAACCTGTGATCTGCTGGGCAGCCAGCAAACCGCCCAGct                                                                                   >  1:1073580/1‑93 (MQ=255)

      tAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGATCAAAACCTGTGATCTGCTAGGCAGCCAGCCAACCGCCCAGctt                                                                                  >  2:4702317/1‑93 (MQ=255)

       aaGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttc                                                                                 <  1:3653557/93‑1 (MQ=255)

       aaGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttc                                                                                 <  2:473677/93‑1 (MQ=255)

        aGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttct                                                                                <  1:3208968/93‑1 (MQ=255)

         gcggggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttctt                                                                               <  1:713009/90‑1 (MQ=255)

            gggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAt                                                                            <  1:1208759/93‑1 (MQ=255)

              gATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAttt                                                                          <  2:3367447/93‑1 (MQ=255)

              gATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAttt                                                                          <  1:3068644/93‑1 (MQ=255)

              gATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAttt                                                                          >  2:257385/1‑93 (MQ=255)

              gATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAttt                                                                          <  2:158613/93‑1 (MQ=255)

                ttACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgc                                                                        <  1:348871/93‑1 (MQ=255)

                ttACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgc                                                                        <  2:1975769/93‑1 (MQ=255)

                 tACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgcg                                                                       >  1:3354985/1‑93 (MQ=255)

                    cGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAt                                                                    >  2:3487205/1‑93 (MQ=255)

                      gTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAttt                                                                  <  1:1002900/93‑1 (MQ=255)

                      gTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAttt                                                                  <  2:2750260/93‑1 (MQ=255)

                        cTTTTTCTACAGCCACGAAGGACGCTAAGGGTCAAAACCTGTGAGCTGCTGGGCAGCCAGCCATCCGCCCAGCTTCTTGATTTGCGCATTttc                                                                <  1:626397/93‑1 (MQ=255)

                         ttttACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttcc                                                               <  1:4482766/93‑1 (MQ=255)

                         ttttACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttcc                                                               <  1:3763629/93‑1 (MQ=255)

                          tttACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttcct                                                              >  2:1353529/1‑93 (MQ=255)

                          tttACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttcct                                                              <  2:4032257/93‑1 (MQ=255)

                          tttACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTctcct                                                              <  2:2951333/93‑1 (MQ=255)

                            tactacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttccttc                                                            <  1:964062/93‑1 (MQ=255)

                            tactacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttccttc                                                            >  1:684921/1‑93 (MQ=255)

                             actacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                           >  1:329068/1‑93 (MQ=255)

                              ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttccttcct                                                          >  2:2139556/1‑92 (MQ=255)

                              ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          <  2:224563/93‑1 (MQ=255)

                              ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          >  1:3005038/1‑93 (MQ=255)

                              ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          <  1:476552/93‑1 (MQ=255)

                              ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          <  1:1229779/93‑1 (MQ=255)

                               tacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCAt                                                         >  1:2180884/1‑93 (MQ=255)

                               tacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCAt                                                         <  1:2840121/93‑1 (MQ=255)

                                tcAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCAtt                                                        <  1:3456622/92‑1 (MQ=255)

                                acAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCAtt                                                        >  2:419252/1‑93 (MQ=255)

                                acAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCAtt                                                        >  2:419253/1‑93 (MQ=255)

                                acAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCAtt                                                        >  2:4472327/1‑93 (MQ=255)

                                 cAGCCCGGAAGGACACTAAGGTTCAAAACCTGGGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTc                                                       <  1:3988179/93‑1 (MQ=255)

                                 cAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTc                                                       <  2:3222511/93‑1 (MQ=255)

                                  aGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCa                                                      >  2:1257517/1‑93 (MQ=255)

                                   gCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCaa                                                     <  1:3013231/93‑1 (MQ=255)

                                    ccACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAAt                                                    <  1:4381668/93‑1 (MQ=255)

                                     cACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATa                                                   <  1:4134753/93‑1 (MQ=255)

                                         aaGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCt                                               <  1:3301753/93‑1 (MQ=255)

                                         aaGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCt                                               <  1:3380528/93‑1 (MQ=255)

                                          aGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTg                                              <  2:659315/93‑1 (MQ=255)

                                           ggACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGt                                             >  2:3490228/1‑93 (MQ=255)

                                             acacTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCt                                           <  2:2668299/93‑1 (MQ=255)

                                              cacTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTg                                          <  1:3817848/93‑1 (MQ=255)

                                              cacTAAGGTTCAAAACCTGTGATATGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATACAATAACCTGTCTg                                          >  2:400846/1‑93 (MQ=255)

                                               acTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTgg                                         >  1:469198/1‑93 (MQ=255)

                                               acTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTgg                                         <  2:994506/93‑1 (MQ=255)

                                               acTAAGGTTCAAAACCTGTGATCTGCTGGGCAACCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTgg                                         >  2:620395/1‑93 (MQ=255)

                                                cTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgc                                        >  1:199274/1‑93 (MQ=255)

                                                cTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgc                                        <  1:246629/93‑1 (MQ=255)

                                                cTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgc                                        <  1:2930973/93‑1 (MQ=255)

                                                cTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgc                                        >  1:2117241/1‑93 (MQ=255)

                                                  aaGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgcgc                                      >  2:22383/1‑93 (MQ=255)

                                                   aGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGcgcg                                     >  1:3782296/1‑93 (MQ=255)

                                                   aGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGcgcg                                     <  1:4542119/93‑1 (MQ=255)

                                                   aGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGcgcg                                     <  2:3042003/93‑1 (MQ=255)

                                                     gTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTc                                   <  1:1615778/93‑1 (MQ=255)

                                                      ttCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTcc                                  >  2:1107187/1‑93 (MQ=255)

                                                       tCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTccc                                 <  2:3074334/93‑1 (MQ=255)

                                                         aaaaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGa                               >  1:2314834/1‑93 (MQ=255)

                                                           aaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATc                             <  1:3508951/93‑1 (MQ=255)

                                                            aCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATcc                            <  2:2758807/93‑1 (MQ=255)

                                                            aCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATcc                            <  2:2319786/93‑1 (MQ=255)

                                                            aCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATcc                            <  1:2007094/93‑1 (MQ=255)

                                                             ccTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCa                           <  1:2413985/93‑1 (MQ=255)

                                                              cTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAg                          >  1:1083176/1‑93 (MQ=255)

                                                                gtgATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTc                        <  2:4584069/93‑1 (MQ=255)

                                                                  gATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCTAGTccc                      >  1:4055518/1‑93 (MQ=255)

                                                                   aTCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTcccc                     >  1:973109/1‑93 (MQ=255)

                                                                    tCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCg                    >  2:788224/1‑93 (MQ=255)

                                                                    tCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCg                    <  2:3907332/93‑1 (MQ=255)

                                                                      tgctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGc                  <  2:2769651/93‑1 (MQ=255)

                                                                       gctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCa                 >  1:2269168/1‑93 (MQ=255)

                                                                       gctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCa                 >  1:2877704/1‑93 (MQ=255)

                                                                         tgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgc               >  2:449370/1‑93 (MQ=255)

                                                                           ggCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTGCCGATCCAGTGCCCGGCAgctg             >  2:818034/1‑93 (MQ=255)

                                                                            gCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgctgc            >  2:1588407/1‑93 (MQ=255)

                                                                            gCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgctgc            <  2:1952500/93‑1 (MQ=255)

                                                                               gccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc         >  2:2913769/1‑93 (MQ=255)

                                                                               gccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc         >  2:2266502/1‑93 (MQ=255)

                                                                               gccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGCCCCCGGCAGCtgctgc         >  2:3687752/1‑93 (MQ=255)

                                                                               gccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGCCCCCGGCAGCtgctgc         >  2:3687750/1‑93 (MQ=255)

                                                                                ccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCGAGTCCCCGGCAGCTGCTGcc        <  1:387798/93‑1 (MQ=255)

                                                                                ccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGcc        <  2:3290463/93‑1 (MQ=255)

                                                                                 cagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCa       <  2:2617300/93‑1 (MQ=255)

                                                                                   gccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGa     >  1:3890490/1‑93 (MQ=255)

                                                                                     caACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAAc   <  1:4606398/93‑1 (MQ=255)

                                                                                      aaCTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAAcg  >  1:3087070/1‑92 (MQ=255)

                                                                                      aaCCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACt  >  1:3450571/1‑93 (MQ=255)

ACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACT  >  NC_000913/3817067‑3817245

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